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Experiments have been performed at Chicago with two Compton-Bennett cosmic-ray 
meters fo determine the frequency of occurrence of very large bursts (120 to more than 1000 
particles) as a function of lead shielding thickness. The thickness of lead for optimum burst rate 
was found on both meters to be in the neighborhood of 3 centimeters. This is larger than the 
optimum for showers and small bursts, a fact which is in qualitative agreement with the theory 
of Carlson and Oppenheimer. Quantitatively, the maximum determined experimentally is of the 
same order of magnitude as that predicted by the theory though lower. 


URING the past few years many experi- 

ments have been performed on the effect of 
shielding upon cosmic-ray showers and small 
bursts. Much less work has been done where very 
large bursts are involved. It has, therefore, 
seemed worth while to report here in brief some 
experiments carried out at intervals during the 
past two years upon very large cosmic-ray 
bursts. These results are in general agreement 
with the work on showers and small bursts and 
are principally of interest because the bursts 
measured are probably as large as any which 
have been thus systematically studied. 

The measurements were made at Chicago with 
two type C cosmic-ray meters of the Carnegie 
Institution of Washington which have been 
described in detail elsewhere.! The meters were 
both placed in a shop building directly under 


1A. H. Compton, E. O. Wollan and R. D. Bennett, 
Rev. Sci. Inst. 5, 415 (1934). 


the roof in such a position that no side wall was 
less than fifteen feet from the instrument and 
the roof was at no point closer than twenty feet. 
The roof was a thin one of slate supported by 
steel girders. The ionization chamber of each 
instrument consisted of a spherical chamber 35 
cm in diameter, the walls of which were } inch 
(1.25 cm) steel. These chambers were filled with 
argon at a pressure of approximately fifty atmos- 
pheres. The lower hemisphere of the ionization 
chamber was at all times submerged in lead shot : 
the upper was covered with a succession of 
hemispherical caps, each of } in. (0.32 cm) 
thickness. For any given shielding thickness, the 
number of bursts occurring during a given 
interval of time was read from the recording 
tape and hence the average burst rate could be 
computed. 

The results for meter No. 4 and meter No. 6 
are given in Table I. The smallest burst which 
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could be definitely distinguished from a random 
fluctuation on the trace was arbitrarily taken as 
one giving a sharp displacement of 1 millimeter. 
This corresponds to the formation of 23 X 10° ion 
pairs within the chamber. Above this minimum, 
bursts were observed ranging in magnitude up 
to a few scattered ones giving over 400 X 10° ion 
pairs. In Table I, the number of bursts observed 
in each size range is indicated and below it in 
parenthesis the rate of occurrence in bursts per 
hour. The small number of bursts counted 
within each size range precludes any great 
accuracy for the distribution curve of the number 
of bursts with size of burst. This distribution is, 
however, in agreement with previous work?: * on 


TABLE |. Summary of data. 
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THICK- 
NESS OF 
LEAD DISTRIBUTION OF BURST SIZE IN ION PAIRS 
cm Hours PRODUCED 
23-45 | 45-68 | 68-90 | 90-113 over 
x 10) | 10° | X10° | 11310 Total 
Meter No. 6 
0 422 21 14 1 1 0 36 
(0.050) | (0.033) | (0.002) | (0.002) | (0.00) | (0.085) 
0.95 210 15 10 4 1 0 30 
(0.07 1) | (0.048) | (0.019) | (0.005) | (0.00) (0.143) 
1.90 165 21 11 1 1 3 37 
(0.127) | (0.067) | (0.006) | (0.006) | (0.02) | (0.224) 
2.5 216 30 14 3 4 3 54 
(0.139) | (0.065) | (0.014) | (0.018) | (0.014) | (0.25) 
2.9 219 32 14 7 5 6 64 
(0.146) | (0.065) | (0.032) | (0.023) | (0.028) | (0.292) 
3.2 240 28 23 2 3 3 59 
(0.117) | (0.096) | (0.008) | (0.012) | (0.012) | (0.246) 
3.8 66 7 4 1 1 3 16 
(0.106)§| (0.06) | (0.015) | (0.015) | (0.045) | (0.242) 
10.7 752 70 27 7 2 3 109 
(0.093) | (0.036) | (0.009) | (0.003) | (0.004) | (0.145) 
Meter No. 4 
0 314 3 2 4 1 0 10 
(0.01) | (0.006) | (0.013) | (0.003) (0.032) 
0.95 381 14 9 4 3 4 34 
(0.037) | (0.024) | (0.105) | (0.008) | (0.01) | (0.089) 
1,90 364 21 10 2 2 3 38 
(0.058) | (0.027) | (0.005) | (0.005) | (0.008) | (0.104) 
2.9 488 60 28 12 4 5 109 
(0.123) | (0.057) | (0.025) | (0.008) | (0.01) | (0.223) 
3.8 504 41 18 8 2 4 73 
(0.081) | (0.036) | (0.016) | (0.004) | (0.008) | (0.145) 
10.7 290 18 13 3 3 3 40 
(0.062) | (0.045) | (0.01) | (0.01) | (0.01) | (0.138) 


2? R. L. Doan, Phys. Rev. 49, 107 (1935). 
3D. Heyworth and Ralph D. Bennett, Phys. Rev. 50, 
589 (1936). 
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Fic. 1. Curve showing relation between burst rate and 
thickness of lead shielding for bursts of all sizes observed. 


the same type meter. The total burst count and 
burst rate, including bursts of all sizes, are given 
in the last column. The burst rates for this 
column are plotted in Fig. 1 against the shielding 
thickness of lead. The statistical uncertainty, 
indicated by the length of the vertical lines 
through the points, is rather large. In addition, 
because of the scattered points, there is some 
uncertainty in the drawing of the curve for meter 
No. 4. Both curves are observed to have a 
maximum in the neighborhood of 3 cm of lead. 
The curve for meter No. 6, however, lies through- 
out somewhat above that for meter No. 4. 
The explanation for this is not apparent. The 
effect, however, has been observed before with 
two similar meters.” 

If the number of ions produced by a particle 
per centimeter path in air at atmospheric pres- 
sure is taken as 60 ion pairs per centimeter,‘ 
then from the geometry of the system the 
ionization produced by each shower particle 
traversing the chamber is computed to be in 
the neighborhood of 100,000 ion pairs. This is 
at best only an approximation. Hence for our 
minimum burst size we have about 230 particles 
traversing the chamber. For the largest bursts 
observed the number of particles is of the order 
of 4000. 


DISCUSSION 


In comparing these results with the results 
obtained by others, we meet with the usual 
difficulty that the geometry of the systems used 
by various experimenters is quite different. In 
addition to this, the fact that the wall of our 


*W. F. G. Swann, Phys. Rev. 44, 961 (1933). 
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jonization chamber was quite thick (1.25 cm of 
iron) presents an added complication in the com- 
parison. The effect of this wall thickness on the 
burst production observed might be twofold. 
First, we add an additional thickness of shielding 
material. This would probably amount to an 
equivalent thickness of only a few millimeters of 
lead, the exact value depending upon the form 
of the law® used to represent the relative burst 
frequency in iron and in lead. Second, an addi- 
tional uncertainty would also be introduced by 
the possibility of lead-iron transition phenomena. 
These effects have been discussed by others® 
but the exact correction in the present experi- 
ments cannot be accurately evaluated. It might 
be pointed out that such phenomena may 
possibly be present as a disturbing factor, to a 
greater or less degree, in any experiments where 
the wall of the ionization chamber is of a dif- 
ferent material from the shielding material. No 
attempt has been made to correct the data for 
either of these two effects. 

The results of the present experiment are in 
good agreement with earlier? * less quantitative 
experiments on the type C meter. The curve in 
Fig. 1 is also in agreement with the results of 
Carmichael’ for bursts of comparable size (160- 
720 particles). Nie,’ however, for about the same 
range of burst sizes finds an optimum shielding 
thickness of between 4 and 5 centimeters of lead. 


1937} M. Nielsen and J. E. Morgan, Phys. Rev. 52, 568 
( 

6C. G. Montgomery, D. D. Montgomery, W. F. G. 
Swann, Phys. Rev. 47, 512 (1935). 

7H. Carmichael, Proc. Roy. Soc. 154, 223 (1936). 

8H. Nie, Zeits. f. Physik 99, 453 (1936). 


Finally, the present maximum at about 3 
centimeters of lead is definitely larger than 
the maximum below 2 centimeters found by 
others’: for showers of only two or three 
particles and of about 2 centimeters found for 
small bursts" consisting of 30 particles or less. 
This increase in the optimum thickness of 
shielding with the size of the burst has already 
been noted experimentally by the authors cited 
within the region of showers and small bursts, 
and is qualitatively in agreement with the 
multiplicative theory of showers of Carlson and 
Oppenheimer.” Because of the large uncertain- 
ties both in theory and experiment, it is im- 
possible to make an exact quantitative com- 
parison between the theoretical prediction and 
the present experimental value for the optimum 
lead thickness. Suffice it to say, that for a burst 
of 450 particles, about the mean size of the 
bursts observed above, a thickness in the neigh- 
borhood of 4 centimeters would be expected. 
This is of the same order of magnitude though 
larger than the present experimental result. 
The authors wish to express their thanks to 
Professor A. H. Compton for placing at their 
disposal the facilities for this work and for his 
continued interest in it. They also wish to thank 
the Carnegie Institution of Washington for the 
use of the meters employed for the experiments. 


®D. K. Froman and J. C. Stearns, Phys. Rev. 52, 382 


(1937). 

10M. A. Starr, Phys. Rev. 53, 6 (1938). 

"R. T. Young, Phys. Rev. 52, 559 (1937). 

® J. F. Carlson and J. R. Oppenheimer, Phys. Rev. 51, 
217 (1937). 
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A representation of all known data on cosmic-ray intensi- 
ties in their dependence on geomagnetic latitude and on 
altitude, respectively, depth below sea level, has been 
obtained on the following assumptions: 1. A soft component 
consisting of electrons which follow the Bethe-Heitler 
theory up to the highest energies with a primary distribu- 
tion E~" (E=energy, n between 3 and 2). 2. A hard com- 
ponent of different nature which in addition to normal 
ionization is absorbed by a discontinuous process leading 


also to the production of numerous secondaries with a 
cross section inversely proportional to E for high energies, 
The primary distribution in energy is about the same form 
as for the soft component, and also the absolute numbers 
of hard and soft primaries are of the same order of magni- 
tude. This picture is furthermore well compatible with the 
observed secondary and transition effects, especially the 
Rossi transition curves for showers, and their occurrence 
below thick layers of heavy material and underground. 


I. THe DuAcistic NATURE OF COSMIC 
RADIATION 


URING the last years it has become more 
and more clear that there are two distinct 
components in the cosmic radiation, the hard and 
the soft one. They both have electronic charge 
(positive and negative) but are distinguished by 
their penetrating power through heavy material, 
for instance, lead. The soft component is ab- 
sorbed in a few cm of this material and produces 
secondaries (showers) in profusion, whilst the 
hard component can pass through one meter 
without being reduced by more than about 50 
percent. The soft component can be well ac- 
counted for by the assumption that it consists of 
electrons (and photons) which behave according 
to the theory of radiation. 

It is indicated by the following facts that the 
hard rays cannot possibly consist of electrons of a 
different, higher energy range for which the 
theory does not hold any more: 

1. The ionization-curvature measurements of 
Anderson and Neddermeyer' and others which 
seem to show (though it can hardly be taken as 
definitely proven) the existence of radiating and 
nonradiating particles of the same curvature. 

2. The good representation of many shower 
phenomena by the ‘multiplication theory” as 
worked out by Carlson and Oppenheimer? and 


* A report on this paper has already been given at the 


Indianapolis Meeting of the American Physical Society, 
December 30, 1937. 

1 Anderson and Neddermeyer, Phys. Rev. 50, 263 (1936); 
51, 884 (1937). 

? Carlson and Oppenheimer, Phys. Rev. 51, 220 (1937). 


Bhabha and Heitler* which demand the validity 
of the theoretical formula for radiative effects to 
very high energies (even in lead) together with 
the observation that only a small fraction of all 
ionizing rays produce showers.* 

3. The absence of equilibrium between the 
two components as found by Rossi> and Auger 
and others® and discussed by Heitler.” 

4. The absorption curves at high altitudes and 
different latitudes, as will be discussed in 
Section IT. 

In view of this situation the question arises 
whether it is possible to arrive at a complete 
description of all the cosmic-ray phenomena by 
ascribing to the hard component a definite set of 
properties. The following attempt seems, as a 
whole, to give a simple picture of the majority of 
the observed facts without being in contradiction 
to any particular detail. 

It should be realized, however, that such a 
phenomenological approach can only be of a 
preliminary nature until a definite theory has 
been found which connects the high energy 
phenomena with our other knowledge of the laws 
of physics. Until this is accomplished it is not 
possible to prove the uniqueness of a particular 
picture as at least two functions, viz. the primary 
distribution and the absorption and secondary 
effects in dependence of energy, can be adjusted 


3 Bhabha and Heitler, Proc. Roy. Soc. A159, 432 (1937). 

4 This point is made particularly. clear by the extensive 
investigation of Starr, Phys. Rev. 53, 6 (1938). 

5 Rossi, La radiation cosmique (Hermann, Paris, 1935). 

6 Auger, Ehrenfest and Leprince Ringuet, J. de phys. et 
rad. 7, 58 (1936). 

7 Heitler, Proc. Roy. Soc. A161, 261 (1937). 
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to give the desired results. It is, therefore, not 
surprising that previous attempts have failed so 
far to converge, and also the present one should 
not be taken as more than a working hypothesis 
which might possibly lead to a better under- 
standing of the connections between the observed 
phenomena. 

In the next section we give firstly a discussion 
of the high altitude measurements in order to be 
able to separate off the contribution of the soft 
component. 


II. DiscussIon oF THE H1iGH ALTITUDE DATA 


The hypothesis of two different components in 
the cosmic radiation, one of which consists of 
electrons and photons behaving according to 
theory, is strongly supported by the intensity 
versus altitude data at different geomagnetic 
latitudes. 

Figure 6 (full lines) contains the observational 
results as available to the author. Curves I and II 
have been taken from the ionization curves of 
Bowen, Millikan and Neher® (B.M.N.) at 
geomagnetic latitudes \= 3° and 38° converted to 
vertical incidence by the well-known Gross® 
transformation 


V(x) =I(x) —x(dI(x)/dx), (1) 


where x is the depth below the top of the atmos- 
phere (in mass units), J the intensity for all- 
around incidence (as obtained by an ionization 
chamber) and V(x) the corresponding vertical 
intensity. Curve III is the one measured by 
Pfotzer'’ with a counter telescope at \=49°. The 
units are those used by B.M.N.° Pfotzer’s data 
have been expressed in approximately the same 
units by multiplying his values by a factor 2.7 
which gives the correct latitude effect at a 
pressure of 6 to 7 m H,O." 


§ Millikan, Neher and Haynes, Phys. Rev. 50, 992 
(1936); Bowen, Millikan and Neher (B.M.N.) Phys. 
Rev. 52, 80 (1937). 

* Gross, Zeits. f. Physik 83, 217 (1933). 

10 Pfotzer, Zeits. f. Physik 102, 23, 41 (1936). 

" This factor gives a latitude effect of ~10 percent at 
6.7 m H.O between \=49° (Pfotzer curve) and \=38° 
(transformed B.M.N. curve at Fort Sam Houston). It 
should be emphasized that the data contained in Fig. 6 
cannot be considered astoo accurate. The transformation 
(1) magnifies any error of the original measurements, 
especially when the second term involving a differentiation 
becomes larger than the first one; this is the case for all 
x>4m H;,0 whilst for smaller x the measurements them- 
selves are necessarily less certain. The assumption of 


The minimum energy for vertical incidence at 
a latitude \ is given by Stoermer's formula 


Emin = 18 cos* 10° ev. (2) 


Hence curve I represents the effect of un- 
charged primary radiation and of charged pri- 
maries above about 18 billion electron volts. 
II includes the charged primaries down to 6.7 and 
III to 3.1 billion electron volts. 

As the percentage of ‘the soft component 
increases strongly with altitude® it is suggestive 
to compare the observations first with the results 
which would follow from a suitable primary 
electron distribution alone. Such a comparison” 
is given in Fig. 1. The full line gives the difference 
between the curves II and III of Fig. 6 corre- 
sponding to the energy range “‘C”’ from 3 to 
6.7X10° ev, and the barbed line the difference 
between I and II, the energy range ““B”’ from 6.7 
to 18X10° ev. The dotted curves are calculated 
from a primary electronic distribution which be- 
haves essentially like E~?* (E=energy) in the 
range from 3 to 18X10° ev. The details of the 
calculation are given in the Appendix. The calcu- 
lated values are the same as for the soft compo- 
nent in Table III for the corresponding energy 
ranges multiplied by 1.35. For the upper part of 
the curves the agreement is quite good and cer- 
tainly within the limits of accuracy of observa- 
tion and theory. The relative heights of the peaks 
can, of course, be adjusted by the choice of a 
corresponding primary distribution so that the 
only check on theory is given by the shape of the 
curves which is not very sensitive against a 
change of the primary distribution within the 
intervals. One sees that the position of the maxima 
is given correctly in both energy ranges and that 
the curves for interval ‘‘C’’ are steeper and nar- 
rower than for interval “B.’’ Below about 4 m 
H,O both the theoretical curves fall off much 
uniform all around incidence on which (1) is based is 
furthermore not quite correct owing to the complicated 
nature of the geomagnetic effect, but it seems not to be 
worth while at the present stage to attempt a more elab- 
orate discussion. Though the points of B.M.N. follow a 
much smoother curve than those of Pfotzer, curves | and 
II are not inherently better than III. The difference 
III-II, of course, is still less certain owing to the empirical 
conversion of the units. 

#2 A similar comparison with nearly the same results has 
already been made by B.M.N. (reference 8). The argument 
is repeated here as our conclusions are somewhat different 


from theirs and in order to obtain data which are needed 
later. 
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Fic. 1. Field sensitive vertical intensities (dashed, and 
dot-dashed curves, calculated from primary electron 
distribution ; solid and barbed curves, observed). 


more rapidly than the observed ones. According 
to theory there is no perceptible influence farther 
down than 6 m H,O whilst a marked geomagnetic 
effect is known to reach down to sea level 
(10 m H,O) and even somewhat below. This 
discrepancy cannot be removed, at least not in a 
simple way, by the assumption of a failure or 
break-down of the theory. The fit at small 
depths is equally good and the disagreement for 
larger depths is equally pronounced for both 
intervals so that the long tail of the empirical 
data cannot be ascribed to a break-down at high 
energies. A ‘‘stretching”’ of the curves along the 
abscissa (corresponding to a change in the 
constants of the formulae of radiation theory) to 
obtain more intensity at lower depths would 
make the upper regions entirely too broad and 
shift the maxima too far down. Fig. 1 in con- 
junction with the observations mentioned in 
Section I suggests thus strongly a practically 
quantitative validity of radiation theory up to at 
least 18X10° ev plus the existence of a separate 
component of higher penetrating power. 

If we then adopt this interpretation we have, 
of course, to reduce somewhat the amount of soft 
radiation in the upper reaches of the atmosphere. 


We can further conclude that the hard com- 
ponent alone is responsible for the geomagnetic 
effect at sea level, that it also Mcreases strongly 
with elevation and that it is still responsible for 
the greater part of the geomagnetic effect up to 
about 4 m H,O. 

A few remarks should be added about the third 
interval “‘A’’ pertaining to curve I, Fig. 6, i.e., 
primary energies over 18X10° ev and contri- 
butions from uncharged radiation (neutral 
particles and y-rays). For this practically every 
law of absorption at low depths can be obtained 
by a suitable choice of the primary distribution. 
The observed absorption coefficient of 0.5 m-! 
H.O for the soft component which near sea level 
is nonfield-sensitive would be given by a primary 
electronic distribution proportional to E-** at 
high energies, as shown by the author, and a 
smooth combination of this with an initial E-** 
law for lower energies accounts for practically 
every property of the soft component." 7 It 
should be mentioned, however, that according to 
this view the rest of the soft component near sea 
level is due to electrons of primary energies of 
about 10" ev and over. It is quite possible, and 
even probable that the radiation theory which 
seems to hold well up to some 10"° ev will have to 
be modified in this region. For instance, the finite 
mass of the atomic nuclei and their recoil 
momenta in radiative processes would have to be 
taken into account. The same final absorption 
law might therefore be due to a more rapidly 
decreasing primary distribution together with 
smaller radiation probabilities at these extremely 
high energies. 


III. GENERAL PROPERTIES OF THE HARD 
COMPONENT 


We consider now a number of facts pertaining 
to the hard component, which seem to allow some 
significant conclusions. 

1. The apparent absorption of the hard rays is 
nearly proportional to the mass of the traversed 
layers, certainly not to Z’. 

2. At sea level about 70 percent of the total 
radiation belongs to the hard group. The fraction 


* caused by it may even be somewhat higher as it 


also produces some soft radiation (showers) by 


13 Nordheim, Phys. Rev. 51, 1110 (1937). 
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secondary processes. The large majority of the 
more energetic rays over 10° ev is certainly hard 
(i.e., they do not radiate appreciably). 

3. The energy distribution at sea level has 
been measured by Blackett." He finds an about 
equal number of positives and negatives. Leaving 
out the particles below 5.10% ev the number of 
which is not very accurately determined and 
which seem to be mostly electrons, the center of 
gravity of Blackett’s distribution lies around 
3X 10° ev with a slow tailing off at high energies. 
The particles with energies above 10!° ev form 
about 10 percent, the ones above 20 X 10° ev (not 
deflected by his magnetic field) about 5 percent 
of the total. According to 2, we can take these 
figures as representative for the hard component 
alone. 

4. The geomagnetic effect at sea level shows 
that 80 to 85 percent of the total, and therefore, 
(when allowance is made for a field-insensitive 
contribution due to the soft component) that 
around 70 percent of the hard rays are not 
influenced by the magnetic field of the earth. 
This means, that so much of it must be due to 
primaries over 18 X 10° ev (or neutral radiation). 
On the other hand, the geomagnetic effect at sea 
level extends up to 50° latitude which means that 
primaries of 3X10° ev must already be able, 
with some probability, to penetrate through the 
atmosphere. 

From 3 and 4 it seems necessary to conclude 
that also the majority of the hard rays at sea 
level are not the primaries themselves but of 
secondary nature.'® Though the majority of the 
rays is due to primaries over 18X10° ev their 
average energy at sea level is of the order of 
3X10° ev. This difference of about 15X10° ev 
cannot be due to a continuous slowing down 
process as this would be in contradiction to the 
extension of the geomagnetic effect to 50°; and it 

4 Blackett, Proc. Roy. Soc. A159, 1 (1937). 

1% Also in case that part of the hard rays were due to an 
uncharged primary radiation they had to be classified, of 
course, as secondaries. There is, however, no other evidence 
in favor of such an assumption which would only make the 
mn more complicated. On the other hand, a considerable 
raction of the hard radiation, larger than indicated by the 
geomagnetic effect at sea level alone, must be due to 
charged primaries, as the east-west effect persists even at 
the equator, and also the simple geomagnetic effect in- 
creases considerably at moderate altitudes, where the soft 
component should still not be appreciably field sensitive. 


For data on the east-west effect compare Johnson, Phys. 
Rev. 48, 287 (1935). 


would furthermore be difficult to reconcile such a 
large energy loss with the great penetrating power 
of the hard component after having come down to 
sea level. The degeneration of the high primary 
energies must be effected by a discontinuous 
process, leading to very extended straggling. 

Further indications for the largely secondary 
nature of the hard rays at sea level are given by 
the following observations. The east-west effect 
shows that the majority of the field sensitive 
rays are due to positive primaries whilst Blackett™ 
has found no significant deviation from a half-and- 
half division. The absorption curves of the hard 
component in lead interposed between counters 
seem to be practically the same at all latitudes 
and also below a thick layer of rock'® (of about 
60 m water equivalent) which means that their 
energy distribution must be nearly the same 
under these different conditions. Lastly the 
existence of a long range secondary radiation has 
directly been proved by the transition effects for 
showers under thick layers of heavy material 
(second maximum of the Rossi curve) and the 
transition effect air to water or rock as. will be 
discussed in Section IV. 

As a last significant fact we take the con- 
tinuation of the absorption curves to depths 
below sea level, investigated extensively by 
Ehmert!’? and Wilson.'* It is found that the 
intensity as a function of depth is very well given 
by simple power laws over very large ranges. 
For the first 20 m water equivalent below sea 
level a transition effect'® apparently takes place. 
From about 20 m to 250 m (Ehmert'’s end point) 
both Ehmert’s and Wilson’s data follow very 
closely the function x~" where x is the total mass 
of the layer above including the atmosphere and 
the exponent lies between 1.8 and 2 (Ehmert 
gives 1.87). From 250 m to about 1500 m 


1 Measurements by Morgan and Nielsen, now in 


progress. 

17 Ehmert, Zeits. f. Physik 106, 751 (1937). 

18 Wilson, Phys. Rev. 53, 337 (1938). The author is 

eatly indebted to Professor Compton and Dr. Wilson for 

ind communication of their data. 

19 This transition effect could formally be made to dis- 
ed by ascribing to the atmosphere a water equivalent 
of about 15 m instead of the 10 m according to its mass; 
i.e. with this assumption the ~x~? law would hold im- 
mediately from sea level on. Such an interpretation seems, 
however, to be unlikely as there is no significant deviation 
between Ehmert’s measurements in water and Wilson's 
in rock, and it would be difficult to understand, why there 
should be a considerable difference between air and water, 
but none between water and rock. 
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(Wilson’s end point) the exponent is somewhat 
higher (about 2.4) corresponding to a more rapid 
decrease. Some single points by Barnothy and 
Forré”’ fall in between the extrapolated 
curve and Wilson’s observations. 

These results are the more surprising in view of 
the discontinuous absorption mechanism dis- 
cussed above. From this one would have expected 
an exponential type law but the observations do 
not show the least resemblance to the latter as the 
apparent absorption coefficient 


k= —1/V(dV/dx) 


=intensity) changes from 0.1 m~' H,O near 
sea level to about 0.003 m~ at a depth of about 
1000 m H,O. Such a slow absorption and ap- 
parent hardening of the radiation can, however, 
result also from a discontinuous absorption 
mechanism, provided the absorption cross section 
decreases with energy. The combination of such a 
behavior with a suitable primary distribution 
tailing off slowly to high energies might, of 
course, produce any kind of absorption curve. 
When we admit, however, that the validity of 
simple power laws over such large ranges is more 
than a mere coincidence, we can limit the 
possibilities so far, as to give at least a model for 
the observed behavior without too much arbi- 
trariness. This shall be attempted in the next 
section where various schemes will be discussed. 
The absorption law described above consti- 
tutesarather strong argument against a hypothe- 
sis put forward by Heisenberg”! that the cosmic 
radiation is brought down to the greatest depths 
by neutrinos. According to his ideas, the 
neutrinos, generated themselves as secondaries 
in the atmosphere and higher strata of the 
surface of the earth, would be distributed over a 
comparatively narrow energy band, and their 
absorption by conversion into ionizing tertiaries 
should consequently follow an exponential law. 


IV. Tue Apsoreprion LAW FOR THE HARD 
COMPONENT 


We try now to bring the qualitative con- 
siderations of Section II11—which seem to be well 
founded—into a quantitative form, which neces- 
sarily will be much more uncertain. 


2° Barnothy and Forré, Zeits. f. Physik 104, 744 (1937). 
* Heisenberg, Zeits. f. Physik 101, 533 (1936). 


The hard rays firstly lose energy by normal 
ionization, an effect which certainly must be 
taken into account, when treating the penetra- 
tion to great depths. At a distance x from the top 
of the atmosphere the energy of a ray of initial 
energy E will be 


E(x) =Eo—Bx, (3) 


where £8 is the specific energy loss. 8 can be 
considered as a constant since the density of 
ionization seems to be practically independent of 
energy. 

We furthermore were led to assume an ab- 
sorption coefficient «(#) dependent on energy. 
The chance for a primary of energy Eo to survive 
down to a depth x will then be given by an 
attenuation factor 


A(Eo, x) =exp | 


0 


Eo 
f K(E)dE, 3] (4) 
E 


o—Br 


and the number of rays of energy FE from a 
primary distribution N(£o, 0) will be 


N(E, x) = N(E+8x, 0) 


xexp| f (5) 


The energy lost in the absorption processes 
will reappear through the production of second- 
ary radiation which, as it seems, is very im- 
portant. We can obtain an estimate of this effect 
on the basis of the following assumptions. 

1. The secondaries have a shorter range than 
the primaries so that they can be considered to be 
in equilibrium with the primary beam. This 
appears to be well justified in our case where we 
have an absorption coefficient decreasing with 
increasing energy. (It is, however, absolutely 
wrong for electrons and photons.) 2. The dissi- 
pation of the primary energies is ultimately 
always effected by ionizing particles all of the 
same ionizing power 8, regardless of their nature. 
This too seems to be very nearly correct as the 
fraction of heavy ionizing tracks and of actual 
nuclear disintegrations in the cosmic radiation 
seems to be quite small. 
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With these two assumptions the number S(£) 
of secondaries in equilibrium with a primary 
beam of energy E will be simply the energy lost 
by absorption from the primary bundle over the 
energy loss of a single ray by ionization, i.e., 


S(E) = Ex(E) 8. (6) 


This formula will hold also in case the secondaries 
are of different nature (hard and soft rays). It 
is irrelevant whether the secondaries are emitted 
singly or in bundles. Eq. (6) will also remain 
correct in case of further multiplication of the 
secondaries (even with nonionizing links) if only 
the total conversion of energy takes place in a 
short enough distance. Of course, it is not 
possible to give a more detailed discussion of 
transition effects between different materials 
without more specific statements regarding the 
production of secondaries. 

With these assumptions the total number of 
rays penetrating to a depth x becomes 


Z(x) = N(E+8x) exp] f 
0 


E=0 


X(1+Ex(E)/B)dE. (7) 


Z(x) is empirically found to be close to x~* (see 
Section III). Such a behavior can be obtained 
from (7) in various ways, as there are two 
adjustable functions and «(£)). 

The following cases may be noted. 


(a) If the absorption is neglected (x=0), i.e., 
if only ionization is present, a primary distribu- 
tion N(E,))dEy~ E~"dE will give a Z(x) 

(b) If the energy loss through ionization is 
neglected (8=0), a power law for the absorption 
coefficient together with a power 
law N(E,) = will give asymptotically 


Z(x)~x- 


(c) For the combination of the two effects a 
simple solution is possible only for 


xk=«,/E, (8) 


N(Eo, 0)((Eo— Bx) /Eo)*(1+0) 


Eo—Bx 
I(Eox) =+ N(Eo, exp| 
4 Bo 


LN (Eo, 0) exp /Bo} 


i.e., an absorption coefficient inversely pro- 
portional to energy. This gives for the attenua- 
tion factor (4) 


0 
(4a) 

and for N(E») = CE 

C (o+1)'(n—2)! 


(ax)"" 


Z(x) = 
(otn—1)! 


(7a) 


For n=3 this would be 
Z(x) = (7b) 


The law (8) is practically the only one which 
combines the previously deduced qualitative 
features without being too complicated or 
artificial and seems, therefore, to be very suitable 
for a discussion of the hard component. But (8)* 
would give rise to very large cross sections at low 
energies. To avoid this singularity we make the 
further assumption that the absorption cross 
section becomes constant below a certain energy 
Bo, i.e. 

ki/Bo for E<Bo, 


E)= 


This law should not be taken too literally. It 
means simply that below a certain energy 8» the 
behavior of the hard rays can be represented by 
a sort of average absorption coefficient which 
might well conceal a more complicated depend- 
ence on energy while above 8, a definite tendency 
to decrease with energy is present. The change in 
the low energy region does not influence the 
asymptotic behavior at large depths. 

For the practical application of the above 
hypothesis we define an influence function 
I(Eo, x) giving the number of rays (primaries and 
secondaries) at a depth x, due to primaries of 
initial energy Eo. Using (7) and (8a) we obtain 


for By 


for Eyg—Bx<Bo<Eo, (9) 
0 


for Bx < Ey <Bo. 
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The total intensity due to primary energies above 
a given E,,in becomes 


Z(x) = ° I(E, x)dE. (10) 


Emin 


By proper choice of E,,in (compare (2)) we can 
evaluate also the geomagnetic effect. In (9) 8 
represents the specific energy loss through ioniza- 
tion, Bo the critical energy in (8a) where the law 
of absorption changes, and o=«x,;/8 the average 
number of secondaries in equilibrium with a high 


energy primary. 


V. NUMERICAL REPRESENTATION OF CosMIC-RAY 
INTENSITIES AT ALL LATITUDES 
AND ALTITUDES 


With the hypotheses about the hard component 
discussed in the preceding section their primary 
distribution for medium and high energies is 
ealready determined by the absorption curve 
below sea level as CE-" with n between” 2.8 and 
3. We extrapolate this distribution also to low 
energies as the simplest possible assumption. To 
obtain numerical agreement with the data and 
units given in Section II we take 


N(Eo, 0) = (10/E)*** (11) 


(E> in billion electron volts). The absorption 
curve below sea level is then taken care of 
automatically.* To adjust the two available 
constants in the absorption law (8a) we use the 
geomagnetic effect at sea level. The following 
values have finally been assumed 


22 According to the observations of Wilson (reference 18) 
the absorption exponent eng a to slightly higher values 
below 250 m H,0. From this follows a somewhat stronger 
falling off of the energy distribution at extremely high 
—. It is not necessary, to take this into consideration 
for the discussion of the behavior of the cosmic radiation in 
the atmosphere, as about 60X10° ev are required from 
the ionization losses alone before a ray can penetrate to 
these depths. 

23 The transition effect air underground (water or rock) 
mentioned in Section III might be due to any of the follow- 
ing causes or a combination of them: (a) deviation of the 
primary distribution from (11) at lower energies. A change 
of the exponent to 2.5 below ~15 X 10° ev would be suff- 
cient. (b) The values of the constants in (8a) in different 
media might be different. A change in the product 
«,8 (compare (7b)) by about 30 percent would be sufficient. 
(c) Inadequacy of the assumption of equilibrium between 
primaries and secondaries, i.e., the secondaries at a given 
depth would already be created in a higher stratum with 
greater primary intensity. (d) Change of counting efficiency 
due to the surroundings by dense material. It seems not 


. possible, at present, to decide between these possibilities. 


energy loss through ionization 
Bo=0. 3X10°ev/m 

critical energy at which absorption law 
changes By) =9X10°% ev 

absorption coefficient below 8; (12) 
Ko=0, 5/m H.O 

absorption coefficient above By; 
=x,/E=4. 5/E (Ein 10° ev) 

from this the average number o of sec- 
ondaries for E>8o;¢=«,/B=15. 


Fig. 2 shows the influence curves (9) calculated 
for the values (11) and (12) for the depths 
x=10 (sea level), 5 and 2 m HO. The insert 
gives the geomagnetic effect at sea level obtained 
by integrating the influence curve from the 
corresponding minimum energies (2) on. It has 
been made* to be 20 percent between the 
equator and the plateau reached at \=50°. (A 
contribution of 2 in our units for the soft com- 
ponent which is independent of latitude at sea 
level has been added to the total.) 

The calculated curve gives very well the type 
of the observed effect. The smallness of the effect 
is due to the greater penetrating power and 
larger ratio of secondaries to primaries of high 
energy, the pronounced beginning at 50° to a 
large number of primaries of low energy together 
with a comparatively high absorption coefficient. 
The beginning of the plateau at 50° is, of course, 
due to our choice of 8 with which just 3X 10° ev 
are required for a vertical traversal of the 
atmosphere. The irregularity at \= 30° (from the 
dip in the influence curve) comes from the break 
in the functional behavior (8a) of the absorption 
coefficient. It could be removed, of course, by 
rounding off (8a). As a matter of fact, such an 
irregularity in the geomagnetic effect seems 
actually to exist, though, of course, the interpre- 
tation by the change in the absorption coefficient 
is highly doubtful. We have given no experi- 
mental points for comparison as the geomagnetic 
effect has not been measured very exactly with 
counters (i.e., vertical incidence). Also the 
observations with ionization chambers (which 
give a somewhat smaller value) show large 
fluctuations due to seasonal and local influences*® 
(irregularities in the earth magnetic field). 

24 This corresponds to the He given by Johnson and 


Read, Phys. Rev. 51, 557 (193 
2% Compton and Turner, Phys. Rev. 52, 799 (1937). 
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Fic. 2. Influence curves for the hard component. Insert, 
geomagnetic effect at sea level. 


With the constants thus determined the 
intensity of the hard component can be calculated 
for all latitudes and altitudes. This has been done 
by drawing the influence curves (9) for various 
depths (some of which are shown in Fig. 2) and 
numerical integration. For the long drawn out 
high energy tail a suitable transformation has 
first been made. The obtained values are just of 
the kind required by the high altitude analysis in 
Section IT. 

For the soft component the distribution previ- 
ously determined from Pfotzer’s data’: * were 
used as described in the appendix. Its form is 
nearly the same as for the hard component, i.e., 
an initial E~** law, which only falls off more 
slowly at high energies, the exponent being —2.3 
in the region over 10" ev. 

The relative numbers of hard and soft pri- 
maries can be seen from the first column (x =0) of 
Table III in the appendix. They are of the same 
order of magnitude with the electrons slightly 
predominant. 

The results of these calculations are shown in 
Figs. 3 to 6. The first three give the energy 
ranges “C” of Section II with the 
contribution of hard component alone, the sum of 
the hard and soft components as calculated (dot- 


ted lines), and the observational data (full lines). 
It is seen that the hard component too shows a 
strong increase with altitude especially in the 
low energy range. A synopsis of the results is given 
in Fig. 6. The fit is extremely good and certainly 
within the accuracy of theory and experiment. At 
very high altitudes the calculated values fall 
somewhat below the observed ones. An additional 
soft primary y-radiation (of energy below 3X 10° 
ev) of small intensity would just account for the 
differences shown in all three curves but as all the 
data, the experimental as well as the theoretical 
ones, are least certain in this region, its reality 
can hardly be taken as proven. The numbers of 
particles outside the atmosphere before the 
secondary effects set in are given in Table III. 
They are much smaller than at the intensity 
maxima, in qualitative agreement with the 
recent observations of Curtiss and others.*® 

The data observed in the balloon flights of 
B.M.N. and Pfotzer are, of course, less accurate 
at low altitudes as their apparatus are adapted to 
the high intensities farther up. They can be used 
only to about 6 to 7 m H,O. From 5 to 6 m down 
to sea-level ionization chamber measurements 
have been given by B.M.N.?’ and Regener and 
Hoerlin.?* From a Gross transformation of their 
curves (which does not give very accurate results 
in this region) the author finds a ratio of between 
3.5 and 4 for the vertical intensities between 
10 m and 6.7 m H,O in good agreement with the 
calculated figures. (Compare Table III.) On the 
other hand, Auger and others® report a ratio of 
not more than 2 as measured by counter tele- 
scopes. The reason for this disagreement is not 
clear. It might be due to a change in efficiency of 
the counter set used. Further precision measure- 
ments both with counters and unshielded ioniza- 
tion chambers at small and moderate altitudes 
would, therefore, be very desirable. 

* Curtiss, Astin, Stockmann, Brown and Korff, Phys. 
Rev. 53, 23 (1938). The value of about 1/150 of the maxi- 
mal ionization reported by these authors seems, however, 
to be far too low. It is to be noticed that the maximum 
intensity at our latitudes is, at least, twice as high as near 
the equator and this difference has to be accounted for 
by charged primaries. This would require the production 
of approximately 75 secondaries per primary in about 2 m 
H,0 which can neither be obtained from radiative multi- 
plication nor any other known effect. 

—-* Millikan and Neher, Int. Conf. Phys. London 


28 Regener and Hoerlin, Int. Conf. Physics, London 
1934; Hoerlin, Zeits. f. Physik 102, 652 (1936). 
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Fic, 3. Vertical intensity due to primaries over 18 X 10° ev. 
(Range 


A discussion of the measurements obtained 
with shielded ionization chambers would be 
possible only with a more detailed knowledge of 
the secondary processes of the hard component 
which seem to be of quite an involved nature (see 
following section). The general results for a 
shielded chamber in comparison to an unshielded 
one, i.e., a much less rapid increase of ionization 
with altitude (due to stronger absorption of the 
soft component in the shield) and an at least as 
high latitude effect at moderate altitudes are in 
qualitative agreement with our picture.”® 

Our formulae give thus a very good represen- 
tation of the cosmic-ray intensities at all lati- 
tudes, altitudes and depths and also of the 
relative influence of the hard and soft com- 
ponents. They allow us to retain fully the 
multiplication theory of showers and explain 
therefore the general parallelism*® between the 
intensities of showers and the soft component as 
already pointed out by Heitler’ and the author.® 

Some remarks should be added on how far the 
assumed constants could be varied without 
destroying the agreement with experience. With 

29 The filtering of part of the soft component will in- 
crease the geomagnetic effect below 5 m H,O. The absorp- 
tion of the hard component in the shield will tend to 
decrease the geomagnetic effect. The balance of these two 
effects might produce an apparent insensitivity of the 


geomagnetic effect against thickness of shielding as re- 
rted by some observers. 

30 The showers produced by the hard component cer- 
tainly play an insignificant role in the upper atmosphere. 
At sea level, however, we estimate that between 1/4 and 
1/2 of the soft rays and showers are actually generated by 
the hard component (which are included in our figures for 
the latter). The small latitude effect for showers as re- 
ported recently by Neher and Pickering, Phys. Rev. 
52, 111 (1938) might be well due to this portion. The 
absence of an exact parallelism of their curve to the one 
for single rays (if it really exists) cannot be taken as an 
argument against our view as the shower production of the 
hard rays might well depend on energy. 


wrensity 
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Fic. 4. Vertical intensity due to primaries between 6.7 and 
18 10° ev. (Range “*B.’’) 


the power law (11) for the primary distribution 
and the geomagnetic effect at sea level the 
constants are actually fixed within rather narrow 
limits. The comparatively high value for the 
absorption constant «x; in (8a) (and with that the 
high value for the average number of secondaries 
a) is necessary to account for the smallness of the 
geomagnetic effect near sea level. The high value 
for By has then to be taken in order to keep the 
absorption coefficient for primaries of low and 
medium energy sufficiently down to give them a 
chance to penetrate through the atmosphere. 
With this the amount of hard radiation in the 
atmosphere is completely determined. Only the 
constant 8 for the linear energy loss could be 
changed in rather wide limits. It has been found 
that all calculated intensities remain the same 
within a few percent for all values of 8 between 
0.2 and 0.4X10° ev/m H,O (the assumed value 
was 0.3) provided the other constants (x1, 8») are 
kept fixed. (¢ = x:/8 has, of course, to be changed 
accordingly.) 

However, one has to remember that the pri- 
mary distribution (11) is made necessary by the 
absorption curve underground only for energies* 
above 15 to 20X10° ev. A less pronounced 
increase in the number of primaries below this 
energy would lead to smaller values for «;, ¢ and 
By and a somewhat smaller geomagnetic effect 
from the hard component at moderate altitudes. 
It would furthermore offer an explanation for the 
transition effect air-underground (compare 

31 The maximum of the influence curve for a depth x oc- 


curs at the energy Emax = 8xo/n (n=exponent in primary 
distribution) and the greatest contribution comes from 


energies> Emax. 
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reference 22). These features would facilitate the 
representation of the observational data.* We 
have, however, refrained from introducing such 
modifications in order not to make too many 
hypotheses ad hoc. 


VI. THe SECONDARY EFFECTS OF THE HARD 
COMPONENT 


We have, so far, not attempted any detailed 
description of the absorption mechanism of the 
hard rays. From our general discussion we 
deduced a number o=15 of secondaries® in 
equilibrium with energetic hard primaries, and 
it remains to be seen whether such a high value is 
compatible with our other knowledge. It allows 
firstly a natural explanation of Blackett’s energy 
spectrum." Practically all the rays below 10 to 
1510° ev will be secondaries according to our 
view and the ratio of rays below and above this 


INTENSITY 
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Fic. 5. Vertical intensity due to primaries between 3 and 
6.7 10° ev. (Range 


® B.M.N. (Abstract, Phys. Rev. 53, 329, (1938)) have 
reported that their newly measured absorption curves at 
\=50° and 60° show only a very small difference even at 
high altitudes. This would indicate definitely a limitation 
of the primary distribution at lower energies. 

% This figure has to be taken as entirely provisional. 
It could well be lowered for the following reasons. (a) 
change in the primary distribution as discussed at the end 
of last section. (b) not all absorbed energy goes into pro- 
duction of normally ionizing rays. (c) the average linear 
energy loss is higher than assumed. 


is of the order 10 to 1. The greater part of the 
total intensities at sea level and underground is 
then carried down by the rays of high energy 
which are accompanied by a secondary spectrum 
which does not depend much on the primary 
energies and is therefore nearly the same every- 
where. The low value of the ionization near the 
top of the atmosphere is also an indication for a 
high number of secondaries, though only of a 
qualitative nature. 

As for the chance of a direct observation of the 
secondary processes we note first that our 
maximum absorption coefficient (holding for the 
majority of observable rays) is of the order of 
0.5 m~ H,O. This means that such an event 
should happen in the average after traversing a 
layer of about 200 g/cm’, i.e., for instance, about 
20 cm of lead. These processes will be in most 
cases indistinguishable from electronic processes 
(forks or showers) and the above figure seems to 


_be low enough to account for the small ratio of 


showers to single rays observed in cloud chambers. 

The hard rays have to be responsible also for 
the showers observed below thick layers of heavy 
material, i.e. the tail end of the Rossi curves, and 
for showers underground. These showers seem to 
be in general also of electronic nature as shown by 
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Fic. 6. Vertical intensity in the atmosphere at different 
latitudes; full lines observed, dotted lines calculated. 


‘ 
‘ 
\ 800 
200 
4 
‘ 
‘ 
‘ 
‘ 
\ ‘ 
100 / 00 H 
/ HARO (CAL) 
U 
400 \ 
\ 
> 8 \ 
/ 
200 ‘ 
365 
1 


PO 


704 Ws 


the transition effects between different materials.™ 
To explain this it is only necessary to assume 
that in some of the absorption processes at least 
one electron or photon can be produced which 
automatically will appear as a shower if it is 
energetic enough. Since the number of soft 
particles in equilibrium with the hard component 
(as found at some meters underground when the 
rest of the effect of soft primaries has apparently 
already been filtered out) seems to be of the 
order 10 to 20 percent, about the same fraction*® 
of the absorbed energy should go into production 
of soft radiation following from the argument in 
Section IV. 

This picture of a composite effect of a hard 
secondary radiation of average range of about 
200 g/cm? and a soft one corresponding to elec- 
tronic nature, finds considerable support in the 
recent experiments of Bothe and Schmeiser.*® 
They find the ‘‘second maximum” of the Rossi 
transition curve to come out much more mark- 
edly*”? for showers of small angular divergence 
so that its existence is now beyond doubt. Such 
a maximum can only be due to a transition effect 
of rays of corresponding range which agrees 
remarkably well with our absorption coefficient 
for the less energetic hard rays. On the other 
hand, the final shower rays, also emerging from 
the region of this second maximum are highly 
absorbable, i.e., they are of “‘soft’’ nature. 


* Morgan and Nielsen, Phys. Rev. 52, 564 (1937) have 
shown that the effects of thin plates of Pb under thick 
plates of Fe and vice versa conform very well to the pre- 
dictions of Carlson and Oppenheimer (reference 2). The 
absence of the first Rossi maximum for showers under- 

ound or below thick layers of light material (concrete, 
ice) is simply an effect of geometry. The multiplication 
distance in these materials (40 cm in H,Q) is still small 
enough so that all soft radiation (electrons and photons) will 
already appear in bundles in contrast to air with a free 
path of several hundred meters. This difference, therefore, 
constitutes an argument for and not against the multiplica- 
tion theory of showers. 

% According to Ehmert and Wilson (references 17 and 
18) the proportion of showers to single rays increases 
somewhat when going to very great depths. This might 
be due to a dependence of shower production on energy 
as the average energy of the rays will increase with depth. 

% Bothe and Schmeiser, Naturwiss. 25, 669, 834 (1937). 
For the proof of existence of a hard secondary radiation 
“a also Maass, Ann. d. Physik 27, 507 (1936). 

37 That this second maximum is, in general, so little 
pronounced must be caused by the fact that the absorp- 
tion coefficient is apparently very nearly proportional to 
the mass of the material so that the distribution of the 
hard rays does not change very much when going to 
another medium. 


NORDHEIM 


The absorption mechanism of the hard com- 
ponent is, thus, certainly of a rather complex 
nature. It seems to be premature to attempt a 
detailed analysis of these secondary effects. We 
therefore leave it open whether the absorption 
mechanism consists in the emission of single 
rays, bundles or even explosion phenomena, 
whether the emission of hard and soft radiation 
occurs in the same or different events, and 
whether the nuclear disintegration processes 
sometimes observed in the cloud chamber are 
representative of its nature or belong to a dif- 
ferent category. 

The chief point of the foregoing considerations 
was to show the compatibility of our analysis 
with the observed effects and it seems even that 
it offers an explanation of many of them in a 
very natural way. 


VII. ConcLusions 


The picture at which we finally arrive might 
be summarized as follows : There are two different 
components, electrons (always accompanied by 
photons) which behave according to the theory 
of radiation, and hard rays. The latter are 
absorbed by a discontinuous mechanism in 
which both hard rays and soft ones are emitted. 
Their absorption coefficient is more or less 
constant up to energies of some billion electron 
volts and decreases thereafter with energy. The 
primary distributions of both kinds are very 
similar** and nearly a power law E~" with n 
between 2 and 3. Also the absolute numbers of 
hard and soft primaries are nearly the same. As 
a matter of fact, with only slight deviations from 
the assumed simple form of the absorption laws 
or with small changes in the observed intensity 
curves, even an exact coincidence between the 


38 This similarity of the two distributions is somewhat 
surprising as they have been derived from quite different 
arguments, for the electrons from the high altitude data 
(shape of the absorption curve and geomagnetic effect) 
and for the hard component from the absorption measure- 
ments underground. The two distributions might even be 
exactly the same within the limits of error of our determina- 
tion. This likeness is suggestive perhaps of a common origin. 
It also seems rather certain that the distribution must fall 
off rather slowly with energy, i.e. not exponentially but 
with a power law. E~* would still be possible according 
to our analysis though we obtain a somewhat better fit 
with £~:§, As from ordinary processes one would expect 
rather a Maxwellian kind of distribution, a power law might 
point to a cosmological origin. 
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two primary distributions might be effected.*® 

The given figures and tables show that it is 
possible to obtain on these lines a more or less 
quantitative description of the intensities of both 
components at all depths, altitudes and latitudes, 
and at least a qualitative understanding of the 
Rossi transition curves. 

In judging the success of this attempt, it 
should be borne in mind that throughout the 
simplest possible assumptions have been made, 
and complications such as irregularities in the 
primary distributions and in the behavior of the 
absorption coefficient have been avoided. Though 
the qualitative features seem in most cases to be 
well founded, the numerical details will in course 
of time certainly need considerable readjustment, 
and it might well be that some of the facts which 
we accepted as fundamental will turn out to be 
a premature generalization. 

With all due reserve as to the final meaning of 
our speculations a few general remarks might be 
added. This analysis does not say anything about 
the nature of the hard particles. We only have to 
conclude that their mass must be large compared 
to the mass of an electron as otherwise they 
would radiate strongly*® and the energy loss 
through radiation would increase with energy, in 
contrast to the behavior of our absorption coef- 


ficient. A mass of about a hundred times the elec- 
tronic mass, as suggested by some _ recent 
observations,“ would be amply sufficient while 
a protonic nature seems only to be excluded by 
the ionization curvature measurements of Ander- 
son,' Street and Stevenson, and the failure to 
find slow protons in sufficient amount in the 
cosmic radiation.” 

In view of the experimental uncertainty about 
the hard component two other possible inter- 
pretations should be mentioned. It is also com- 
patible with our present knowledge that the 
primaries belonging to the hard rays are of a 
different nature than the secondary rays observed 
near sea level. To be more specific they might be 
protons to which the absorption coefficient «,/E 
might belong. They, in turn, would generate a 
new particle to which belongs our average ab- 
sorption coefficient at lower energies. Or else the 
hard rays might be created by electrons of high 
energy in addition to ordinary radiative processes 
in which case most of them would come from the 
upper region of the atmosphere where the elec- 
tronic intensity is highest.“4 From the point of 
simplicity, however, it seemed to be desirable, 
to discuss the hard rays as they appear, i.e., as a 
definite separate primary group. 


APPENDIX: NUMERICAL DATA 


The calculation of the contribution of the soft com- 
ponent has been carried out in the same way as by Heitler.” 
The total number of electrons (negatives and positives) 
with an energy larger than E, due to a primary electron of 
energy Eo after traversing a layer of thickness /, can be 
represented by a function Z(/, y), where y=log Eo/E, 
tabulated in Table I. 


39 The number of primaries (column x=0 in Table II1) 
is not very exactly determined by our procedure and might 
easily be in error by a factor of order 2. In using the Bhabha- 
Heitler influence function Z(l, y) (compare Table 1) the 
low energy electrons with E<0.15X10° ev have been 
neglected. As the distribution in energy of the soft rays 
does not change much with altitude (with the exception of 
the very highest layers) the inclusion of these rays of low 
energy will approximately give a constant factor to the 
soft total and the number of primaries should be reduced 
accordingly. The number of hard rays, on the other hand, 
depends essentially on the average number of secondaries 
o which is also rather uncertain (comp. footnote 33) and a 
lower value of ¢ would require a higher number of pri- 

. maries. Both these effects would work in the direction of 
equalizing the two primary distributions. 

‘© Though the energy loss through radiation cannot con- 
tribute much to the absorption coefficient it might pos- 
sibly be responsible for some of the shower production. 


Some new values have been added to those given by 
Heitler, part of which have been calculated, the others 
having been obtained by logarithmic interpolation. For the 
atmosphere the unit for / is 280 m standard air=0.35 m 
water equivalent. The final energy E has to be chosen so 
that there is no appreciable multiplication below it, i.e., 
for air E=0.15 X 10° ev. 

Our adopted primary distribution for the soft component 
is listed in Table II, also in a logarithmic energy scale. 
The total number of rays at depth /, coming from a certain 


Street and Stevenson, Phys. Rev. 52, 1003 (1937); 
on Takeuchi and Ichimiya, Phys. Rev. 52, 1198 
(1937). 

# Street and Stevenson, Phys. Rev. 51, 1005 (1937) (A). 

48 Montgomery, Montgomery, Ramsay and Swann, 
Phys. Rev. 50, 403 (1936). 

“ This point of view has recently been strongly advo- 
cated by B.M.N., Phys. Rev. 53, 217 (1938). The argu- 
ments used by these authors are, however, very similar to 
some in this paper and entirely compatible with our an- 
alysis. This shows again the impossibility to prove the 
uniqueness of a phenomenological picture as long as the 
physical laws regarding the hard component are unknown 
(note added after completion of the manuscript). 


L. W. NORDHEIM 


TABLE I. Influence function for electrons Z(I, y). 


\ 
: 3 3.5 4 4.5 5 6 7 8 9 10 11 12 13 
2 2.1 2.6 3.4 4.4 5.8 10 17 21 26 31 33 34 
3 2.7 3.6 5.0 6.9 8.9 15 26 43 63 76 89 96 
4 2.4 3.6 5.4 7.9 12 22 38 59 84 140 190 250 
5 2.2 3.2 4.9 7.4 13 25 51 98 160 250 360 550 
6 1.3 2.2 4.0 7.2 13 27 66 120 220 370 620 | 1000 
7 0.9 1.7 3.2 6.2 12 32 68 140 300 540 830 | 1400 
8.5 | 0.43 1.0 2.2 4.7 9.3 29 73 170 400 830 1600 2800 
10 0.22 0.63 1.5 3.5 6.6 24 76 200 500 1150 2200 4100 
11.5 | 0.1 0.28 0.78 1.9 4.1 18 61 180 540 1250 3000 6000 
13 0.05 0.15 0.40 1.0 2.6 12 46 160 480 1300 3500 8300 
15 0.02 0.06 0.17 0.49 Lz 6.4 27 100 360 1150 3200 9300 
17 0.008 | 0.02 0.07 0.20 0.53 3.1 17 71 260 910 2800 7800 
19 0.003 | 0.008 | 0.03 0.08 0.21 15 8.7 40 170 650 2100 6000 |17000 
22 0.002 | 0.006 | 0.019 0.056 | 0.45 3.0 16 73 310 1050 3000 | 8100 
25 0.001 | 0.004 0.015 | 0.14 0.95 5.6 30 140 500 1600 | 4900 
29 0.002 | 0.023} 0.21 1.5 2.5 42 180 700 | 2300 
TABLE II. Primary electronic distribution F(y)dy = F(Eo)dEo/Eo. 
y 3 {35/4 | 45] 5 6 7 8 9 10 11 12 13 
Eo in 10° ev 3 5 18211335) 22 60 165 450 1200 3300 9000 2.5 X10! | 6.7104 
F(y) 250 | 89 | 33 | 12.8) 5.0 | 0.87 | 0.22 | 0.062 | 4.7 1.3 X 107-3 | 3.6 10-4| 110-5 


TABLE III. Cosmic-ray intensities for the different latitude intervals at various depths. 


DeptH BeELow Top oF ATMOSPHERE IN m H2O 0 2 3.3 5 6.7 10 
Interval A (A=0°) soft 2.2 90 97 49 18 e 
E o> 18 10° ev hard 1.9 22.5 18 14 11 7.4 
total 4.1 112.5 115 63 29 9.1 
: Interval B soft 24 125 62 10 0.5 an 
6.7>Eo>18 10° ev hard 9.5 53 29 13 7 1.6 
total 33.5 178 91 23 75 1.6 
Interval C soft 96 93 48 5 0.3 — 
3>Eo>6.7 10° ev hard 37 95 42 16 6.4 0.7 
total 133 188 90 21 6.7 0.7 
Intensity at \=49° soft 122.2 308 207 64 12.8 2 
hard 48.4 170.5 89 43 24.4 9.4 
total 170.6 478.5 296 107 43.2 11.4 


interval of primary energies, is then 
NO= »ddy. (13) 
To the limits of our intervals A, B, C correspond the values 


A: E=18X10°ev_ corresponding to y=4.75, 
B: E=6.7X10°ev corresponding to y=3.75, 
C: E=3xX10°ev corresponding to y=3. 


The integrals (13) have been evaluated graphically. Table 
III shows some of our final results from which the Figs. 3 
to 6 have been drawn. It gives the calculated intensities for 
the soft and hard component separately for the three 
intervals at various depths below the atmosphere. The 
values for x=0 give, of course, the number of primaries. 
A small fraction (perhaps 15 percent) of the figures for the 


hard group will actually be soft rays, though produced by 
hard primaries and in equilibrium with them. For instance, 
1.4 of the 9.4 for the “hard”’ component at sea level (x = 10) 
might be electrons, which brings the percentage of highly 
absorbable rays up to about 30 percent as observed. At 6.7 
m H;0 corresponding to an elevation of 3500 m, perhaps 
3.5 of the hard intensity of 24.4 might be soft and the per- 
centage of the soft rays is then about 50 percent as it should 
be. The increase of the total intensity between sea level 
and 3500 m is from 11.4 to 43.2 that is 1 to 3.8 in agreement 
with the values deduced from ionization chamber measure- 
ments. Thus our formulae give also the correct behavior for . 
the intensities of the hard and soft component in the range 
between sea level and 3500 m altitude, which is not covered 
by the data summarized in Figs. 3-6. 
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The distribution of fast neutrons about the Li+D 
source has been investigated by means of radioactivity in- 
duced in silver (24.5 min.) and aluminum (14.8 hr.). The 
observations are roughly in agreement with a distribution 
calculated from statistical considerations. The maximum 
neutron energy is 20.8 Mev. An increase in the deuteron 
energy results in an increase in the neutron energy as meas- 
ured by the »—3n reaction of scandium and the n—2n 
reactions of silver and copper. The intermediate nucleus, 
Be®, formed by the deuteron bombardment of lithium must 
have a life short compared to 10~" sec. A number of ma- 


terials when bombarded with high energy detuerons vield 
neutrons of sufficient energy to produce the 7 — 2” reaction 
in silver. Only the neutrons from lithium can produce the 
n—3n reaction in scandium. Neutrons of somewhat less 
energy are required for the » — 2” copper reaction and of ap- 
preciably less for the » — 2” silver reaction. When scandium 
is bombarded with fast neutrons, reactions of the type 
n—3n, n—2n, n—a, n—p and n—y take place in the fol- 
lowing percentages, 24.0, 27.8, 20.6, <0.5, and 27.0, 
respectively. 


INTRODUCTION 


HE fast neutrons from the Li+D reaction 
have been used largely to discover new 
induced radioactive periods in the various ele- 
ments. This use has been very successful and is of 
special interest because these neutrons are very 
proficient in the production of the n—2n type of 
disintegration.!: 

In this paper, however, known radioactive 
periods are used to indicate the distribution of 
these high energy neutrons about the lithium 
source, the dependence, if any, of the neutron 
energies on the bombarding deuteron energy and 
finally the potentialities of other fast neutron 
sources. Branching ratios for the disintegration 
products from scandium are also presented. 


Fast NEUTRON YIELD 


The metallic lithium to be bombarded with 
deuterons was securely fastened upon a water 
cooled copper plate. Directly behind the copper 
plate was placed weighed samples of Cu, Sc, Ag, 
Al and Dy. These samples were all irradiated at 
the same time with neutrons obtained from 
bombarding the lithium with 5.7 Mev deuterons 
produced by the cyclotron at the University of 
Michigan. Other sets of samples were given 
similar treatments but with deuterons of lower 
energies. Aluminum foils were used to decrease 
the energy of the deuterons. 


1M. L. Pool, J. M. Cork and R. L. Thornton, Phys. Rev. 
52, 239 (1937). 
*F. A. Heyn, Physica 4, 160 (1937). 


The intensities of radioactivity induced in four 
sets of samples are shown in Fig. 1. All the curves 
were arbitrarily put through the point 4 on the 
intensity scale. The rise in intensity of the 
dysprosium activity (strong 2.5-hr. slow neutron 
period) with increase of deuteron energy probably 
indicates the rise in the slow neutron background 
which, in turn, is a measure of the increase in the 
total number of Li+D disintegrations taking 
place. The rise in intensity in excess of the 
dysprosium rise would indicate the influence of 
the increased number of faster neutrons. 

From the curves it is inferred that faster 
neutrons are required to produce the 10 min. 
period in copper than the 24.5 min. period in 
silver. Both periods are the result of an n—2n 
type of reaction. The energy of neutrons required 
to produce the 4.0 hr. period in scandium, an 
n—3n type of reaction,’ is probably a little more 
than that required to produce the 10 min. copper 
period. Sagane reports that neutrons of more than 
12 Mev are required to produce the 10 min. 
copper activity. The 12.8-hr. copper activity 
which is produced by both fast and slow neutrons 
is shown and was measured from the same copper 
samples that gave the 10 min. period. The n—a 
type of reaction responsible for the 14.8-hr. 
period in aluminum apparently takes place as 
readily with neutrons of quite low energy as with 
those of much higher energy. 


3M. L. Pool, J. M. Cork and R. L. Thornton, Phys. Rev. 
52, 41 (1937). 
4R. Sagane, Phys. Rev. 53, 212 (1938). 
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INTENSITY OF RaADIORCTIVITY 


1 
4.3 4.8 3.3 5.7 


Deuteron Enency in Mev 


Fic. 1. Intensity of radioactivity induced by neutrons 
from the Li+D reaction versus energy of bombarding 
deuterons. 


The extra energy acquired by the neutrons 
must be derived from the kinetic energy of the 
impinging deuterons. For this to happen it is 
necessary that the intermediate nucleus Be® at 
the time of its disintegration be in motion 
through the lithium metal. A 5.7 Mev deuteron 
gives a velocity of about 5107 cm/sec. to Be’. 
To bring this nucleus to rest in the metal a time 
of 10-" to 10-” sec. is required according to a 
simple calculation based upon estimated range- 
energy relations.® The life of Be® must therefore 
be short compared to 10-" sec. 


DISTRIBUTION OF FAst NEUTRONS 


Strips of silver and aluminum 1 by 13 cm were 
bent in an arc of about 2.9 cm radius. The 
centers of the circular strips coincided with the 
region of the metallic lithium most strongly 
bombarded by deuterons of about 6 Mev. If the 
neutrons ejected from the lithium source had 
spherical symmetry in energy and number, then 
the radioactivity induced in each centimeter 


5 F. Kirchner, H. Neuert and O. Laaff, Ann. d. Physik 
30, 527 (1937). 


POOL 


length of the 13 cm strips would be of the same 
intensity. The intensities however were found 
not to be the same. The asymmetry is shown by a 
polar plot in Fig. 2. 

As may be seen the induced radioactivity in 
the direction of motion of the deuterons is 
definitely greater than in the direction at right 
angles and suggests therefore a marked forward 
throw of the neutrons resulting from the trans- 
mutation of the lithium. The 24.5-min. period of 
silver (Ag!*) and the 14.8-hr. period of aluminum 
(Na**) are known to be produced by only fast 
neutrons and are not induced by slow neutrons or 
stray deuterons. Since the asymmetry of the 
silver radioactivity is greater than that of the 
aluminum it is concluded that the fast neutrons 
responsible for the production of the silver 
activity have a greater asymmetrical distribution 
than those responsible for the aluminum activity. 

Inquiry may be made concerning what energy 
distribution might be expected of the trans- 
mutation products when the intermediate nucleus 
transforms according to 


Li’+ D*Be*—Be$+ n. (1) 


Stephens using a 0.91 Mev deuteron beam ob- 
tained neutrons coming off at 90° which produced 
alpha-particle recoils in a helium filled cloud 
chamber.‘ His results are shown in Fig. 3 and may 
be interpreted as consistent with the above 
reaction. For a 6 Mev deuteron beam the 


@Al Activity (14.8 hes) 


Fic. 2. Polar plot of intensities of radioactivity induced 
in samples of silver (Ag!®, 24.5° min.) and aluminum 
(Na*, 14.8 hr.) placed symmetrically about the Li+D 
source of fast neutrons. 


E.Stephens, Phys. Rev. 53, 223 (1938). 
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FAST NEUTRONS 


neutrons are calculated to be ejected in about 
equal numbers at right angles and in the forward 
direction with energies of 17.3 and 20.8 Mev 
respectively. Since fast neutrdéns over a wide 
range of energies seem about equally efficient in 
producing the m—a reaction in aluminum, the 
above-mentioned two groups of neutrons should 
induce in the aluminum approximately equal 
radioactivity, and give, therefore, a symmetrical 
intensity distribution. However, asymmetry is 
observed. This suggests that reaction (1) does 
not play the dominant role in this experiment. 
The asymmetry is easily accounted for pro- 
vided the transmutation takes place according to 


D’—Be*—He'!+ He*+2!. (2) 


For this type of reaction a statistical treatment 
of the neutron energy distribution may be made 
by assuming for the two alpha-particles and 
neutron, equal elements in phase space are 
equally probable. The treatment of this three- 
body problem is similar to that given by Fokker, 
Kloosterman and Belinfante on the energy distri- 
bution between the products of the transmutation 
of boron atoms into three alpha-particles.’ In our 
problem it is of particular interest to calculate the 
number of neutrons, having energies between E 
and E+dE, that are ejected parallel Noe and at 
right angles No» to the 6 Mev deuteron beam. 
The results are expressed in the following 

equations where £& is the energy of the neutron 
in Mev: 

No? = (17.3+0.77E!— E)(E£)!}, 

Noe = (17.3—£)(E)!. 


These equations are plotted in Fig. 3. The 
maximum energy expected for the forward 
neutrons is 20.8 Mev. The neutron distribution 
curve for zero energy deuteron bombardment is 
also shown and may be compared with the 
distribution observed by Bonner and Brubaker 
by means of proton recoils in a cloud chamber 
placed at right angles to an 0.85 Mev deuteron 
beam.* They observed that about 95 percent of 
the reactions were of type (2). The agreement 
between theory and experiment and between the 
experiments themselves is not good but some of 


7A. D. Fokker, H. D. Kloosterman and F. J. Belinfante, 
Physica 1, 705 (1934). 
(1935) W. Bonner and W. M. Brubaker, Phys. Rev. 48, 742 
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Fic. 3. Calculated energy distribution of neutrons in the 
directions parallel and at right angles to the impinging 
deuteron beam. Observations of Bonner and Brubaker 
(B and B) and Stephens (S) are also shown. 


the discrepancy may be attributed to the uncer- 
tain cross section for high energy neutron 
interactions. It may be necessary also to intro- 
duce into the theoretical treatment, factors other 
than those of pure chance as has already been 
suggested by Cockcroft and Lewis in connection 
with the break up of B"’. 

The ratio of the calculated total number of 
forward neutrons to the total number of right 
angle neutrons is 1.8. The observed ratios of the 
intensity of the radioactivity in silver and 
aluminum in these two directions are 1.43 and 
1.35, respectively. The observed asymmetry in 
the intensity of the radioactivity is therefore 
easily accounted for on the hypothesis of type 
(2) reaction. 


OTHER FAst NEUTRON SOURCES 


A number of materials were bombarded with 
6.3 Mev deuterons and the intensities of the fast 
neutron induced radioactive periods in Ag, Cu, 
Al and Sc were measured. Table I summarizes 
the results. 

As is seen the Li+D reaction furnishes the 
best source of fast neutrons. Again the n—3n 
reaction in scandium appears to require more 
energetic neutrons than does the n—2n reaction 
in copper. The m—2n reaction in silver is next in 


9 J. D. Cockcroft and W. B. Lewis, Proc. Roy. Soc. 154, 


254 (1936). 
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TABLE I. Intensities of radioactivities induced by neutrons 
from various sources. 


NEUTRONS [24.5 min. Ag! 10 min. Cu) 5 min. Cu | 14.8 HR. Al |4 HR. Sc 
FROM n—2n n—2n n-y n-—a n—3n 
Li+D 100 100 —- 57.5 100 
B+D 11.5 12 66 100 —_ 
CaF.+D 3.15 _- 71 59.5 
Ag+D 2.15 _ 20 38 — 
Fe+D 1.9 12 44 
Al+D 75 22 73 
Be+D 100 57.5 


order and the »—a reaction in aluminum requires 
the least energy. 

The maximum energy of the neutrons in the 
forward direction from the Be+D reaction is 
calculated to be about 10 Mev. Since the n—2n 
reaction in silver does not take place with these 
neutrons, it is inferred that more than 10 Mev 
neutrons are required to remove two neutrons 
from silver. There must, however, be neutrons 
with energies greater than this value released 
from Al, Fe, Ag, CaF, and B when they are 
bombarded by 6.3 Mev deuterons. It is inter- 
esting to note that neutrons from the B+D 
reaction seem very efficient in producing the n—a 
reaction in aluminum. 

It is also interesting that deuteron bombard- 
ment of aluminum produces neutrons which are 
energetic enough to then produce quite strongly 
the m—a reaction in aluminum. This tandem 
reaction probably accounts for at least some of 
the long radioactive periods observed in alumi- 
num when given extensive bombardment with 
high energy deuterons.'® Likewise, neutrons from 
the Ag+D reaction are capable of producing, al- 
though weakly, the m — 2n type of reaction in silver. 


BRANCHING RATIOS 


Since an appreciable quantity of scandium was 
used in the above experiments, six samples that 
showed strong radioactivity were chosen to be 
observed over a period of a few months for the 
purpose of calculating the branching ratios or 
relative rates of formation of the disintegration 
products. Two of the samples were followed in 
their chemical separations (K, Ca, Sc). The 
samples had been irradiated from one to four 
hours with neutrons from a lithium target 
bombarded with about 6.0 wA of 6.3 Mev 
deuterons. 


10M. L. Pool and J. M. Cork, Phys. Rev. 51, 383 (1937). 
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TABLE II. Percent of various radioactive nuclei formed in the 
transmutation of Sc® and V5! when bombarded 
by high energy neutrons. 


OF Sucu 
NUCLEUS TYPE OF NucLeus | RADIOACTIVE NUCLEI 
BOMBARDED | REACTION FORMED PERIOD FORMED 
n—3n Sc 4.0 hr. 24.0 
n—2n Sc# 52. hr. 27.8 
12.4 hr. 20.6 
n—p Ca* 2.4 hr. < 0.5 
Sc* 85 da. 27.0 
n—3n min. 0 
n—2n 3.7 hr. 3 
n—a Sc*8 41 hr. 32 
n—p Ti 2.8 min. 65 


Table II shows a summary of the results and 
expresses in the last column the relative rates of 
formation of the various radioactive trans- 
mutation products as the percent of such nuclei 
formed. It is noticed that scandium is transmuted 
most of the time according to the n—2n type of 
reaction. The n— p reaction is very weak which is 
consistent with the findings of Walke."' However, 
his observations that the number of n—3n 
reactions is nearly twice that of the n—2n 
reactions is not confirmed. As observed from the 
six samples in this research the ratio of the 
number of the former reactions to the latter is 
more nearly 6/7. 

It was thought that possibly the branching 
ratios for the neighboring odd nucleus »;V*! 
would be similar to those for :Sc**. Tentative 
values were obtained and are shown also in 
Table II. As may be seen, no similarity is 
apparent. 
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1 Harold Walke, Phys. Rev. 52, 669 (1937). 
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Formation of an Excited He’ in the Disintegration of Deuterium by Deuterons 


T. W. Bonner 


Rice Institute, Houston, Texas 
(Received March 14, 1938) 


The energy distribution of the neutrons from the reaction H?+H*—~He'*+n'+(Q, has been 
studied. Two homogeneous neutron groups with energies of 1.08 and 2.50 Mev have been ob- 
served at 90° to the 0.11 Mev deuterons. The emission of 1.08 Mev neutrons has been found 
to be approximately 1/10 as likely as the emission of 2.50 Mev neutrons. The corresponding 
values of Q; are Q)°=3.29+0.08 Mev and Q;'=1.40+0.11 Mev. The low energy group results 
when the He’ is left excited to a level of 1.89+0.11 Mev. The mass of He’ has been calculated 
from the value of Q;° and the mass-spectrographic value of H?=2.01473. The result is He® 
3.01700. This result indicates that H* may spontaneously disintegrate into He’ with the emis- 


sion of an electron. 


HE neutrons produced in the disintegration 

of deuterium by deuterons were at first 
thought to be monoenergetic. This was concluded 
on the basis of the experiments of Dee! and those 
of Bonner and Brubaker.? These experiments 
indicated that at least most of the neutrons had 
an energy of about 2.5 Mev. Later experiments 
by the writer,* undertaken primarily to study 
neutron-proton scattering, gave a first indication 
that there probably was a second lower energy 
group with considerably less intensity than the 
high energy group. The purpose of the present 
experiments was to investigate carefully the 
neutrons from the reaction: 


n'+ Qi. (1) 


A search was made for low energy neutrons and 
particular care was taken to determine Q; as 
accurately as possible. A precise value of Q, is 
important because it is from this reaction that 
the mass of He® is best determined. The exact 
masses of He*® and H® are particularly important 
as they each are made up of three fundamental 
particles (neutrons and protons). From their 
masses the relative magnitude of neutron-neutron 
forces and proton-proton forces may be deduced. 


EXPERIMENTAL PROCEDURE 


The accelerating tube and d.c. source of po- 
tential were the same as that previously de- 
scribed.’ The neutron energies were determined 


from the ranges of the recoil protons which were 


'P. I. Dee, Proc. Roy. Soc. A148, 623 (1935). 
? Bonner and Brubaker, Phys. Rev. 49, 19 (1936). 
°T. W. Bonner, Phys. Rev. 52, 685 (1937). 


observed in a cloud chamber. The cloud chamber 
was 17 cm in diameter and 4.8 cm deep. It was 
constructed from as small an amount of material 
as possible in order to reduce neutron scattering 
to a minimum. The chamber was of the rubber 
diaphragm type, thus eliminating the material 
otherwise used in a piston. The chamber proper 
was made of a Pyrex glass ring which was 3 mm 
thick. The plate glass top had a thickness of 3 
mm and the brass bottom plate a thickness of 1 
mm. The cloud chamber was connected to the 
operating mechanism, an electromagnetic air 
valve, by a brass tube which was 2.5 cm in 
diameter and 90 cm long. With this arrangement 
only about 5 percent of the neutrons which pass 
through the cloud chamber make collisions before 
entering the chamber. This reduction in the 
number of scattered neutrons reduces the low 
energy tail or background which is ordinarily 
found on neutron energy distribution curves.’ 
Consequently the experimental arrangement de- 
scribed above is well adapted for observing neu- 
tron groups of low intensity which are present 
along with higher energy neutrons. 

The nearest edge of the cloud chamber was 
placed at a distance of 13 cm from the ion spot on 
the target. The neutrons which entered the cloud 
chamber made an angle of 90+11° to the direc- 
tion of the deuterons on the target. The target 
was a 1 mm layer of D;PQ, at the bottom of a 
thin walled brass cup. Under these experimental 
conditions very few neutrons were scattered by 
the target or target holder. 

The cloud chamber was filled with a mixture 
of 85 percent CH, and 15 percent C.H¢ with 


711 


24.0 
27.8 
20.6 
27.0 
0 
3 
32 
65 
ts and 
ates of 
trans- 
nuclei 
muted 
J. M. 
sity of 
| 


712 

60 

2° 

a 

| 

S 


o 


ENERGY OF RECOIL PROTONS IN Mev 


Fic. 1, The differential energy distribution curve of the 
recoil protons in the forward direction. 


alcohol vapor. The expanded pressure in the 
chamber was slightly larger than atmospheric. 
The stopping power of the gas was found by 
measuring the mean length of polonium alpha- 
particles in the chamber. The stopping power for 
alpha-particles of range 3.805 cm was 1.07. This 
gas mixture has an appreciably different stopping 
power for polonium alpha-particles than for the 
highest energy protons which were observed. The 
stopping power for 10 cm protons was calculated 
from the data of Livingston and Bethe‘ and 
is 1.03. 


EXPERIMENTAL RESULTS 


Sixty-eight hundred pairs of stereoscopic pic- 
tures were taken of the cloud chamber when the 
deuterium target was bombarded with 0.11 Mev 
deuterons. On these photographs roughly 4000 
recoil protons were observed and 205 of these 
were projected within the angular cone of 0—-10° 
to the neutron’s direction. The track lengths of 
these 205 acceptable tracks were measured and 
their ranges computed from the calculated 
stopping power of the gas in the chamber. Their 
energies were then computed from the range 
energy curve of Bethe and Livingston.® Fig. 1 
shows the energy distribution of the recoil 
protons. There appear to be two definite proton 
groups with energies of about 1.1 and 2.6 Mev. 
These correspond to two neutron groups of very 
nearly the proton energies. The lower energy 
group has roughly 1/10 the intensity of the 
higher energy group. This ratio is complicated by 
the differing probability of observing tracks of 


‘ Livingston and Bethe, Rev. Mod. Phys. 9, 275 (1937). 
5 Bethe and Livingston, private communication. 
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Fic. 2, The integral range number curve for long range 
recoil protons in the forward direction. 


different lengths in the cloud chamber and the 
varying cross section for neutron-proton collisions. 

A more accurate maximum energy of the 
neutrons was obtained from the integral range 
number curve of Fig. 2. The extrapolated range 
was 10.60+0.30 cm. The calculation of the 
neutron disintegration energies was carried out 
by the method of Livingston and Bethe.‘ The 
experimental conditions correspond to their case 
of “good geometry.” In their symbols, the range 
straggling s for 10.5 cm protons is 2.2 percent. 
The angular straggling y/R was 4.55 percent and 
the total straggling s’’ was 5.09 percent. For the 
thick target correction, 8 was calculated to be 
11.9 and hence xXextr=0.95. Consequently the 
difference between extrapolated and mean range 
is 0.95-5.09-10.60/100=0.51 cm. It follows that 
the mean range is 10.09 cm. This corresponds to 
a proton energy of 2.46 Mev. The correction for 
using recoils which made angles of 0-10° to the 
direction of the neutron is }E2x0o? = 2.5/2 -0.00306 
= 0.04 Mev. Therefore the energy of the neutrons 
coming from the top layer of the target at 90° 
is 2.50+0.05 Mev. The disintegration energy is 
then calculated from the relation Q,;=(4/3)En 
—(1/3)Ey. The value of Q,° is 3.29+0.08 Mev. 
The corresponding value of Q,', calculated from 
the low energy proton group, is 1.40+0.11 Mev. 


DISCUSSION OF RESULTS 


The low energy neutron group corresponds to 


-an excitation level in He*® at 1.89+0.11 Mev 
- above the ground level. One should expect to 


find a gamma-ray with this energy with an 
intensity of about 1/10 that of the neutrons. 
In the experiments described above some high 
energy electron tracks were observed in the 
cloud chamber which probably were produced 
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by these gamma-rays. However their quantum 
energy must be determined before the correlation 
is made certain. 

In view of the excitation level in He*® at 
1.89 Mev it seems reasonable to expect a similar 
level in H* of approximately the same energy of 
excitation. This follows since the level in He* 
should differ from that in H* only by the per- 
turbation caused by the Coulomb field of the 
extra proton. This level in H* should be observed 
as a second lower energy group of protons from 
the reaction 


H?+H?—H'*+H!+(Q2. (2) 


One should expect a second proton group with a 
range of perhaps 5 cm. A careful search for such 
a group is planned, particularly since the results 
of Oliphant, Harteck and Rutherford® give some 
indication of such low energy protons. 

The mass of He* can be most accurately found 
from the disintegration energy Q,° together with 
Bainbridge and Jordan’s spectrographic mass of 


6 Oliphant, Harteck and Rutherford, Proc. Roy. Soc. 
A144 695 (1934). 


*=2.01473 and a neutron mass’ of 1.00893. 
The calculated mass of is 3.01700 +0.00010. 
Since the mass of H® is 3.01705, this indicates 
that H* may be unstable, disintegrating with a 
long half-life according to the reaction: 


(3) 


This predicted instability of H* is consistent with 
the failure of recent attempts to concentrate H® 
from large samples of heavy water.* 

The value of Q2°—Q,° is just the difference 
between the binding energies of He* and HH’. 
For the value of Q,° found above, this difference 
becomes 3.98 —3.29=0.69 Mev. Since the calcu- 
lated value® of the Coulomb repulsion between 
the two protons in He* is very nearly equal to 
the experimental difference of 0.69 Mev, this 
shows that neutron-neutron forces must be very 
nearly equal to proton-proton forces. 


7H. A. Bethe, Phys. Rev. 53, 313 (1938). 

8’ Lord Rutherford, Nature 140, 303 (1937). 

® Bethe and Bacher, Rev. Mod. Phys. 8, 82 (1936); 
Share, Phys. Rev. 50, 488 (1936); Rarita and Present, 
Phys. Rev. 51, 788 (1937). 
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The Energy-Range Relations for Deuterons, Protons and Alpha-Particles 


F. T. RoGers, JR. AND MARGUERITE M. RoGeErs 
The Rice Institute, Houston, Texas 


(Received March 16, 1938) 


The complete deuteron, proton and a-particle energy-range relations as given by Livingston 
and Bethe and as subsequently revised by Bethe to include the latest data, are presented as a 
set of eleven equations. Given the range of a known particle, the corresponding kinetic energy 
of the particle can be computed easily and quickly from these equations. 


N 1937 Livingston and Bethe! published a 

very extensive set of kinetic energy vs. range 
relations for deuterons, protons and a-particles. 
These relations were in the form of graphs, in 
which the particle energies were plotted as 
functions of particle ranges. Later in 1937 these 
graphical energy-range relations were revised? to 
include the recent data of Parkinson, Herb, 
Bellamy and Hudson’ on proton ranges and of 


1 Livingston and Bethe, Rev. Mod. Phys. 9, 266-269 


(1937). 


* See Bethe, Phys. Rev. 53, 313 (1938). 
’ Parkinson, Herb, Bellamy and Hudson, Phys. Rev. 
52, 75 (1937). 


Blewett and Blewett? on a-particle ranges. Now, 
to provide these latest revised energy-range 
relations in a complete and compact form suitable 
for the easy and rapid computation of energies, 
we present the following. 

Let r be the mean range in cm of a particle in 
air at a temperature of 15°C and at a pressure 
of 760 mm Hg; let V be the corresponding 
energy of the particle expressed in millions of 
electron volts. Let the subscripts D, P, and a 
refer the symbols to which they are affixed to 
deuterons, protons and a-particles respectively. 
Then, given a value of rp, rp, or ra (by observa- 
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TABLE I. Precistons of Eqs. (1)-(11). Average deviations. 
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0.003 
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tion or otherwise), the corresponding value of 
Vp, Vp, or V. can be obtained from the following 


equations 
for deuterons having rp < 12, 
logiy Vp = —0.0835+ (0.595 +0.45¢-2" 

+0.0025 sin® (rr, 9) cos (rr/6)) login r; (1) 
for protons having rp< 6, 
logio Vp = —0.3845+ (0.595 +0.45e-" 

+0.0025 sin? (2277/9) cos (rr/3)) logig 27: (2) 
for protons having 30, 
Vp=1.819+0.16758(r—6) 

—0.003(r—6)?+0.00004(r—6)*; (3) 
for protons having rp < 90, 
Vp=4.661+0.08715(r — 30) 

— 0.000454(r — 30)?+0.00000217(r—30)*; (4) 
for protons having 90< rp 140, 
Vp=8.714+0.054964(r—90) 

—0.00014(r (5) 
for protons having 140< 236, 
Vp=11.174+0.044137(r— 140) 

— 0.0000429(r — 140)? 
—0.000000108(r—140)?; (6) 
for a-particles having r.< 0.9, 
V.=0.98517+3.4958r? — 2.46817r° 
+0.03r sin (27r/0.9); (7) 
for a-particles having la 2.4, 
Ve=1.919+1.3167(r—0.9) 
+0.076 sin (27(r—0.9)/3); (8) 
TABLE II. Estimated accuracies of Eqs. (1)-(11). 


RANGE PossIBLE Error | EQUATIONS 

PARTICLES (cm) IN INVOLVED 
Deuterons Oto 0.6 10% 1 
Deuterons 0.6 to 12 0.06 X 10° ev 1 
Protons Oto 0.3 10% 2 
Protons 0.3 to 236 0.03 < 10° ev 2-6 
a-particles Oto 0.35 10% 7 
a-particles | 0.35to 1.05) ev 7-8 
a-particles | 1.05to 3.5 | 0.0310° ev 8-9 
a-particles | 3.5 to 10.5 | 0.01X10° ev 9-10 
a-particles {10.5 to 21 0.03 X 10° ev 10-11 


_‘In these equations the subscripts of the r’s are, for 
simplicity, omitted. 


ROGERS, JR. AND M. M. 
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for a-particles having 7.2, 


V4=3.894+ 1.0906(r — 2.4) 
— 0.07025 (r —2.4)?+0.00357(r—2.4)°; (9) 


for a-particles having 12, 


Va=7.897 +0.67375(r—7.2) 
(10) 


for a-particles having 12 20.8, 


10.799 +0.50865(r — 12) 
—0.000786(r — 12)?—0.000549(r—12)*. (11) 


Each of Eqs. (1)—(6) and (9)—(11) is associated 
with only one of the graphs of Livingston and 
Bethe. Eqs. (7) and (8) are together associated 
with only one of the graphs. 

The majority of the Eqs. (1)—(11) were ob- 
tained by fitting to the revised energy-range 
graphs, functions of suitable forms by a method 
of zero sums ;° Eqs. (1) and (2) were obtained by 
a “‘trial and error’? method, while Eq. (8) was 
obtained simply by inspection. As measures of 
the precisions with which Eqs. (1)—(11) represent 
these revised energy-range graphs of Livingston 
and Bethe, Table I shows for each equation the 
average deviation A (in millions of electron 
volts) of the computed energy-range graph from 
the corresponding revised energy-range graph. 
Obviously the representations afforded by Eqs. 
(1)—(11) are quite good. 

As estimates of the accuracies with which 
Eqs. (1)—(11) give the appropriate energies, we 
may legitimately adopt the estimates by Liv- 
ingston and Bethe! of the accuracies of their 
graphical energy-range relations. Accordingly in 
Table II we have tabulated as functions of 
ranges, the possible errors which (from Liv- 
ingston’s and Bethe’s estimates) might be 
present in the energies computed from Eggs. 
(1)—(11). Table II should always be taken into 
consideration when Eqs. (1)-(11) are used, 
because it acts as a set of necessary restrictions 
which must be applied to computations made 
by Eqs. (1)—(11). 

Finally, it is noteworthy that Eqs. (1)—(11) 
(of course restricted by Table II) are well suited 
for the simple and rapid numerical computation 
of V’s when the corresponding r’s are given, 
especially if tables of x’, x*, e~7, logio x, and sin x 
are available. 

5 See, e.g., Campbell, Phil. Mag. 39, 177 (1920). 
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Magnetic Field Corrections in the Cyclotron 
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It is pointed out that in order to obtain magnetic fields which will allow the ions in the cyclo- 
tron to attain high energies it is first necessary to render the field homogeneous within an ac- 
curacy of about one part in a thousand. The correction (insertion of ‘‘shims"’) for the most 
important source of inhomogeneity, namely the decrease in the field near the edge of the pole 
pieces, is considered. A general method is given for determining the dimensions for any shape 
of these shims. Results are given for the dimensions of shims in the form of a ring of rectangular 
cross section placed at the edge of the cyclotron chamber. The magnetic field resulting from the 
use of such shims is determined and compared with the field of plane parallel pole faces. For the 
special case of a ring whose thickness is 1/16 the magnetic gap, the field is homogeneous within 
the required accuracy up to a distance from the edge of the poles of 0.4 of the magnetic gap. 
The uncorrected field is homogeneous only up to a distance from the edge equal to 0.9 of the 


magnetic gap. 


I. INTRODUCTION 


T has been realized from the beginning of 
cyclotron technique that great homogeneity 
of the magnetic field is absolutely essential for 
the maintenance of resonance between the 
Larmor rotation of the ions and the oscillation 
of the accelerating electric field. A deviation of 
one percent, e.g., will throw the ions completely 
out of resonance after 100 accelerations so that 
they will actually be decelerated instead of 
accelerated. Now the field produced by plane 
parallel magnetic poles is always inhomogeneous 
near the pole edges (cf. Fig. 4 below). Unless a 
great part of the area of the pole faces is to be 
sacrificed it is necessary to introduce ‘“‘shims” 
near the edge which will have the effect of in- 
creasing the field in their neighborhood and thus 
render the field more homogeneous. 

This shimming must be done rather accurately 
because over-correction must definitely be 
avoided. An over-correction, i.e., a magnetic 
field increasing with the distance r from the 
center of the cyclotron chamber, leads to strong 
magnetic defocusing of the ions, as was shown in 
a previous paper,? and therefore to enormous 
losses of intensity. The magnetic field must not 
increase radially at any point except possibly 
near the center where the electric field provides 
some focusing. 

Even more exact design of the field is necessary 

1 Of course, this is apart from the considerations of ob- 
taining very high energies. 


M. E. Rose, Phys. Rev. 53, 392 (1938); hereafter re- 
ferred to as I. 


if very high energy ions are to be obtained from 
the cyclotron. It was shown in I that the impos- 
sibility of simultaneously maintaining exact 
resonance and focusing of the ions imposes a 
limit on the energy obtainable. At the same time 
it was shown that magnetic fields can be devised 
which make this energy limit rather large.’ In 
order to attain the highest energies, a field such 
as is given in Fig. 5 of I should be used. However, 
all such fields are somewhat critical in the sense 
that they necessarily do not differ by very much 
from a defocusing field. Hence the deviations 
from the prescribed field of the field finally ob- 
tained must be kept very small, of the order or 
less than one part in a thousand. 

The success in attaining high energies will 
depend very much on the elimination of unde- 
sirable inhomogeneities in the field. The most 
important source of such inhomogeneity is the 
decreasing field due to the finite size of the pole 
pieces. In fact, since the desired field has been 
defined only inside the cyclotron chamber, the 
shims required to produce this field can be deter- 
mined only if the effect of the pole edges is made 
negligible in the region where the ions move. 
The following section contains the calculation of 
the size of Shims required to compensate for this 
edge effect as far out from the center of the 
chamber as is possible.‘ 

3In the best case, field given in Fig. 5 of I, the energy 
is E=2.12 (VoAZ)* Mev where V4 is the accelerating dee 
voltage in kilovolts, A and Z the atomic weight and 
atomic number of the ion. 


* Of course, there will in general be some inhomogeneity 
due to the nonuniformity usually present in the large iron 
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II. Shims FOR HOMOGENEOUS FIELD 


In the following we shall make two assump- 
tions: The first is that the reluctance of the iron 
can be neglected so that the faces of the magnets 
(or lids of the cyclotron chamber) form equipo- 
tential surfaces.> The second assumption is that 
the shims may be placed on the inner surfaces 
of the lids of the cyclotron chamber rather than 
in the air-gap below and above the lids. 

Since the magnets are supposed to be axially 
symmetrical we should introduce cylindrical 
coordinates with the origin at the center of the 
cyclotron chamber. Then the magnetic potential 
V is a solution of the Laplace equation written 
thus: 

(1) 


where z is the vertical coordinate. However for 
our present purpose, only the field near the edge 
where r is large is of importance. In this case the 
second term in (1) will be smaller than the other 
terms by a factor of order magnetic gap divided 
by diameter of magnetic poles* and hence may 
be neglected. 

In this case we have a two-dimensional Car- 
tesian problem to which we may apply the 
Schwarz transformation method.’ We introduce 
a Cartesian coordinate system in a plane through 
the z axis, with origin at the edge and in the 
median plane. Denoting (radial) distances meas- 
ured inward by x we have 


(2) 
We introduce the complex variable 
f=x+iz (3) 


and the potential U conjugate to V (U=con- 
stant giving the magnetic lines of force). For our 
Cartesian problem we have the Cauchy-Riemann 
relations between the potentials: 


0U/dx=—dV/dz, 0U/dz=dV/dx. (4) 


magnets, and possibly to irregularities arising from machin- 
ing or assembling the magnets, etc. Obviously corrections 
for these inhomogeneities must be carried out empirically 
and it will be assumed that this can be done to the re- 
quired degree of accuracy. 

5 No currents in the cyclotron chamber. 

6 Cf. the result (27) for the field obtained below and the 
remark preceding (15). The ratio of the second and first 
term of (1) is 1/2 times gap/diameter which for the usual 
dimensions is of the order of one percent. 

= e.g., F. Kottler, Handbuch der Physik, Vol. 12, 
p. 480. 
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Hence the general solution of (2) is the imaginary 
part of the complex function W defined by 


= U(x, 2) +7V(x, 2). (5) 


The procedure now consists in the establish- 
ment of a connection by conformal mapping 
between an auxiliary variable ¢ and W (poten- 
tials) on the one hand and between ¢ and ¢ 
(coordinates) on the other. By elimination of ¢ 
the potentials or magnetic field can be obtained 
in terms of the coordinates. The variable ¢ is to 
assume real values on the equipotential surfaces 
and particular values are to be assigned at all 
places where the slope of these surfaces changes 
i.e., at all the corners of the pole faces. Only 
three of these particular values (say 0, 1 and ~) 
are independent, the remaining values being 
completely determined by the special geometry. 
If we regard the pole faces as forming the bound- 
ary of a polygon and if the particular value /, 
is assigned to the mth corner where the internal 
angle of the polygon is a,, then the connection 
between ¢ and ¢ is® 


d¢/dt= (t (6) 


where C; is a constant to be determined from the 
boundary conditions. 

As an example we may consider the case of plane parallel 
poles. In the finite part of the ¢ plane there is, of course, 
one corner on each pole to which the values ¢, = +1 may be 
assigned. At this corner the internal angle a,=3/2. At 
x=0,7=+ wecan set /,=+- and at x= ~, =0, and 
a, =0. Then we have 


dg /dt = —1)4/t. (6a) 


Since in all cases there are only two equipoten- 
tial surfaces we have, regardless of geometry, for 
the connection between W and ¢ (cf. Fig. 1, also 
Fig. 73 reference 7). 


dW/dt=C2/t, (7) 


where again C2 is a constant to be determined 
from the boundary conditions. 

It is convenient to introduce as unit of poten- 
tial the magnetic potential on the (lower) pole 
face so that we have as boundary conditions 


W(t=1)=-1, W(t=-—1)=+42. (8) 
8 The index m can assume a series of discrete or a con- 


tinuous range of values. The latter case is pertinent when 
the slope of the shims varies continuously. 
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Fic. 1. Conformal mapping of the ¢ plane on the W 
plane. The real ¢ axis coincides with the equipotentials 


Then we find from (7) and (8) 


W= —(2/7) log te™*?, (9) 
or by the use of (5) 


From the symmetry with respect to the median 
plane, the values ¢, at corresponding corners of 
upper and lower pole faces are opposite in sign. 
At these corresponding corners the angles a, 
will be the same. In addition, taking ¢,=0 
between the poles at x= ~ at which point a,=0 
we may write (6) in the form 


df /dt=(C1/t)(P tn?) (11) 


where the index m runs over the values on one 
pole face only. 
Introducing the components of the magnetic 
field 
H,=—dV/dx, H,=—dV/dz (12) 


we find from (4) and (10) 
H,—ill,= (13) 
or from (11) 
H,—iH,= — Ci) tn?) (14) 
We use as unit of length half the magnetic gap 
so that for x= «, =0 the magnetic field becomes 
H.=1, H.=0. (15) 


We may then determine the constant C,; and 


obtain 
H,—iH,=M(1 (16) 


Thus it is seen that the magnetic field may be 
obtained as a function of ¢ or from (10) as a 
function of the potentials without integration. 
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From this result we see that the inhomogeneity 
of the field can be expressed as a power series in 
f. Further, since the field is in any case almost 
homogeneous inside the cyclotron chamber, U 
will be proportional to x within very small cor- 
rection terms (cf. 24, 25 below). Hence from 
(10) the inhomogeneity is 


1—J],=k,e-** cos cos 2nz+---, (17) 


where k;, ky: -- are coefficients depending on the 
parameters /, and hence on the geometry. In 
general the coefficient k, will not vanish and for 
plane parallel pole face k;~} so that the in- 
homogeneity at x=1 in this case is about one 
percent. To make the field more homogeneous 
the dimensions of the shims may be made such 
that k; vanishes. In this case the inhomogeneity 
cos decreases very rapidly and 
will be negligible for x larger than unity. 
Application to ring shims 

A great variety of shapes of shims would be 
equally suitable for accomplishing our purpose. 
In many cases it would be possible to make the 
main term in the inhomogeneity (~e~**) vanish. 
Of course, in the various cases, the size of the 
remaining inhomogeneity, i.e., the coefficient of 
e~*** would depend on the geometry in different 
ways. However, since nothing is gained by com- 
plication, we consider here only the simplest case 
of a single ring shim of rectangular cross section 
on each pole face. The outer diameter of the ring 
is to be equal to the pole radius, the width of the 
ring a and its thickness }, cf. Fig. 2. The parame- 
ters ¢, and internal angles a, are: 


tn Qn 
+1(1—6) +1 3/2 
a+1(1—d) +h, (18) 
att n/2. 


In addition at t=+ and at 
+iz, t,=0. From (16) we find 


(19) 


The coefficient of # in the expansion of this 
expression vanishes if 


te? = t,?/(1+4,) (20) 
and with this result we have 
(21) 
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Fic. 2. Conformal mapping of the ¢ plane on the ¢ plane. 
The real ¢ axis coincides with the surfaces of the magnets. 


To make the inhomogeneity small, it is desirable 
that ¢; be not too small. This will mean as may 
be seen below that a should be small; e.g. a~0.1 
would be quite satisfactory. 

‘For we have from (13) 


dt /dt= —(2/nt)(1—P)! 
(22) 


Upon integrating (22) from t=1 to t=t, we 
obtain a as a function of ¢;. Again integrating 
(22) from to we obtain as a function 
of t;, using (20). Eliminating ¢; we have the 
thickness 6 in terms of the width a. This final 
result is given in Fig. 3. 

The magnetic field on the median plane is 
obtained from (22), (13) and (10) with V=0. 


H,=(1+e-*4 
(23) 


40 


Fic. 3. Dimensions of ring shim for homogeneous 
magnetic field. Thickness of ring 6, width a. The unit of 
length is half the magnetic gap. 


ROSE 


To obtain the field in terms of x we have 
/u 


where f(U) is the reciprocal of the expression 
given on the right-hand side of (23). The constant 
is given by 


log 


+ { (25) 


the integration being taken along the surface of 
the magnet where ¢ is real. This constant is only 
slightly geometry dependent and its value is in 
general about 0.2. 

The field obtained in this way is given in Fig. 4 
(upper curve) for a special case. In the same figure 
the field on the median plane for the case of plane 
parallel pole faces is given. It is seen that the 
homogeneity of the field is indeed very much 
improved. The curves also show that the ap- 
proach to the asymptotic behavior (cf. (17), (24)) 
of the field is rather rapid. By retention of only 
the first term in the inhomogeneity (~e~**") the 
error made at x= 1.0 is only 3 percent of the total 
inhomogeneity. Thus for larger x we may write 


(26) 
(cf. (24)) and V=—z. (26a) 
— 


) 10 15 


5 


Fic. 4. Magnetic field on the median plane near the 
edge of the magnets. The unit of field strength is the 
difference of potential between the magnets divided by the 
magnetic gap. The field labeled ‘‘with shims” applies in the 
case of a ring shim of width a=0.095 and thickness 
6=0.125 in terms of half the magnetic gap as the unit of 
length. The field labeled ‘‘without shims’ pertains to the 
case of plane parallel pole faces. Note the change in scale 
at H=0.90 and at x=0. 
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Then we obtain for the magnetic field in this 
region 
H, =1— cos 2rz, 


(27) 


IT, = sin 


For the uncorrected field of plane parallel pole 
faces the field inhomogeneity is one-tenth percent 
at x=1.8. For the field with ring shims in the 
case b=0.125, a=0.095 we find from Fig. 4 that 
this same inhomogeneity occurs at x = 0.80. Thus 
if the exit slit is placed at this distance (0.4 the 
magnetic gap) from the edge, the magnetic field 
over the entire region of motion of the ions will 
be homogeneous within the required degree of 
accuracy. 


Finally, we may return to a consideration of 
the assumptions made at the beginning of this 
section. First of all the assumption of low re- 
luctance of the iron will in general be fulfilled 
rather well. Of course, it is sufficient if only the 
cyclotron lids and not the large magnets them- 
selves be of low reluctance, high permeability 
iron. Secondly, the assumption that the shims be 
placed inside the chamber against the lids 
rather than in the air gap need not impose any 
restriction on the applicability of the results 
obtained here. Since the shims used to make the 
field homogeneous may be inserted at the time 
of construction it should perhaps be not incon- 
venient actually to place them inside the 
chamber. 
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The behavior of partly polarized beams of slow neutrons as regards their precession on passing 
through homogeneous magnetic fields has been investigated. From the experiments it is con- 
cluded that the neutron has a magnetic moment not far from 21/1840 Bohr magneton and 
that the sign is negative. Further, the precession of neutrons inside magnetized iron was 
investigated; it was found that the field accounting for the observed rate of precession is 
more than 100 times the actual field strength H and actually of the order of magnitude of 


the magnetic induction B. 


1. INTRODUCTION 


HAT a neutron should have a magnetic 

moment at all, seems somewhat surprising, 
on account of its being electrically neutral. On 
the other hand, from the magnetic moments of 
the proton!:? (2.5 to 2.8 n.m., 1 n.m.=1 nuclear 
magneton =1/1840 Bohr magneton); and the 
deuteron’: * (0.85 n.m.), a magnetic moment yp, 
of the neutron, of about 2 n.m., can be deduced ;4 


* Now at the Institut de Chimie Nucléaire, Collége de 

rance, Paris. 

10. R. Frisch, O. Stern, Zeits. f. Physik 85, 4 (1933); 
I. Estermann, O. Stern, Zeits. f. Physik 85, 17 (1933); 
I. Estermann, O. C. Simpson, O. Stern, Phys. Rev. 52, 
535 (1937). 

*1. 1. Rabi, J. M. B. Kellogg, J. R. Zacharias, Phys. Rev. 
46, 157, 163 (1934); 50, 472 (1936). 

31. Estermann, O. Stern, Phys. Rev. 45, 761 (1934). 
(1936) A. Bethe, R. F. Bacher, Rev. Mod. Phys. 8, 91, 205 


the sign of yu, should be negative, that is, the 
relative position of spin and magnetic moment 
should be the same as in the (negative) electron. 
A tentative explanation of this moment, based on 
the Fermi theory of beta-decay, has been offered 
by Wick.® 

A way of measuring, at least roughly, the mag- 
netic moment of free neutrons has been pointed 
out by Bloch.* He showed that the magnetic 
interaction between neutrons and electrons must 
have a measurable influence on the scattering of 
slow neutrons by magnetic atoms or ions (pro- 
vided the neutron has a magnetic moment of the 
order of 2 n.m.). Of special interest is the scatter- 
ing of neutrons from a ferromagnetic substance in 


5G. C. Wick, Att. Acad. Lincei 21, 170 (1935); see also 
reference 4. 
°F. Bloch, Phys. Rev. 50, 259 (1936). 
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which, by means of a magnetic field, the magnetic 
momenta of a large number of electrons are ad- 
justed parallel among themselves; then the scat- 
tering should depend on the spin orientation of 
the neutrons relative to that of the electrons. 
One should therefore be able to produce partly 
“polarized” beams of neutrons by sending them 
through a piece of magnetized iron; a second, 
similar piece might be used as an “‘analyzer.” 

The existence of such an effect was demon- 
strated by Hoffman, Livingston, and Bethe;’ 
they found that the total transmission, for slow 
neutrons, of a system of two iron bars (of 1 cm 
thickness each) was 2.3 percent smaller when the 
bars were magnetized antiparallel instead of 
parallel. The amount of the difference was in fair 
agreement with a calculation’ based on the as- 
sumption w,=2n.m.; the agreement may be 
regarded as an argument for the correctness of 
this assumption. Similar results have been re- 
ported by Beyer, Carroll, Dunning and Powers.*—" 

It was soon realized that the more or less ir- 
regular magnetic stray fields prevailing along the 
path of the neutrons would complicate the course 
of the phenomenon by causing frequent reorien- 
tation of the polarized neutrons. On the other 
hand, the systematic study of this reorientation 
might permit one to determine the magnitude 
and also the sign of the magnetic moment uy, of 
the neutron. This was pointed out by Rabi," who 
calculated the reorientation caused by a special 
type of an inhomogeneous magnetic field. Inde- 
pendently the author’ also demonstrated the 
existence of this reorientation effect by using a 
homogeneous field perpendicular to the polarizing 
and analyzing fields, and showed that it behaves 
roughly as one would expect from the assump- 
tion u,=2 n.m. 

Using a similar arrangement, the authors also 
proved" that the sign of u, is actually negative, 


7J. G. Hoffman, M. Stanley Livingston, H. A. Bethe, 
Phys. Rev. 51, 214 (1937). 

5 J. R. Dunning, P. N. Powers, H. G. Beyer, Phys. Rev. 
$1, 51 (1937). 

®P. N. Powers, H. G. Beyer, J. R. Dunning, Phys. Rev. 
51, 371 (1937). 

10P. N. Powers, H. Carroll, J. R. Dunning, Phys. Rev. 
$1, 1112 (1937). 

uT. I, Rabi, Phys. Rev. 51, 652 (1937). 

20. R. Frisch, H. von Halban, J. Koch, Nature 139, 
756 (1937). 

#0. R. Frisch, H. von Halban, J. Koch, Nature 139, 
1021 (1937). 


as expected from the moments of proton and 
deuteron. The same result was found by Powers, 
Carroll, Beyer and Dunning" who used an ar- 
rangement similar to the one suggested by Rabi." 

A third paper of the authors’ was devoted to 
the study of the precession of neutrons inside 
magnetized iron, where the question of some 
theoretical interest is, whether the magnetic 
field strength H or the magnetic induction B (or 
perhaps something in between) is responsible for 
the rate of precession. 

In the following we shall give a more detailed 
account of the experiments the results of which 
have been published in the three preliminary 
notes,!?: 3. 5 and some considerations and calcu- 
lations which have to do with the problem. 


2. GENERAL REMARKS 


In the classical theory a particle possessing a 
magnetic moment uw and an angular momentum J 
precesses in a magnetic field 17 like a gyroscope, 
with an angular velocity w=//p/J. According to 
quantum mechanics, however, this picture can- 
not in general be used, and in particular fails 
entirely in cases like the Stern-Gerlach effect 


Fic. 1, Experimental arrangement for the demonstration 
of the precession of neutrons. The path of the neutrons 
from the hole in the paraffin block to the ionization chamber 
was almost entirely surrounded by boron and cadmium 
sheets (not shown in the figure) to suppress scattered 
neutrons. 


4 P, N. Powers, H. Carroll, H. Beyer, J. R. Dunning, 
Phys. Rev. 52, 38 (1937). 
a bs “ R. Frisch, H. von Halban, J. Koch, Nature 140, 360 
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Fic. 2. Schematic diagram, illustrating principle of pre- 
cession experiment. 


where we have to do with a different motion of 
the particles in the separate stationary states 
corresponding to different orientations of the 
spin axes relative to the magnetic force. Still, as 
can be simply verified by a direct quantum- 
mechanical analysis,'* the precession picture re- 
tains, in accordance with the correspondence 
argument, its validity in all such cases where, as 
in the phenomena considered below, the transla- 
tive motion of the particles is not essentially in- 
fluenced by their angular orientation in the field. 
A particle with angular momentum 3/ and mag- 
netic moment yu in a magnetic field /7 precesses 
with an angular velocity w= 2y///h. If the veloc- 
ity of the particle is v, then it moves a distance 
lp=v/w=vh/2uH, while rotating through one ra- 
dian. Assuming »=2 n.m.=1.10~* c.g.s. units, 
we find w=2-10* H; for an average ‘‘thermal” 
neutron with v=2- 10° cm/sec. we get /p- 7=10, 
where /z is measured in cm and H in gauss. 
This value of /z-H is very much larger (about 
10‘ times) than the values encountered in experi- 
ments on the reorientation of atoms.'? In molec- 
ular beam experiments, unless special arrange- 
ments are made, the direction of the magnetic 
field changes only very slightly over a distance Jp 
along the path of the particles, and the spins 
therefore follow the magnetic field lines adiabat- 
ically, precessing within a very narrow cone 
around them. Only if the field is made very weak 


wn e.g., C. G. Darwin, Proc. Roy. Soc. 117, 258 

"OQ. R. Frisch, T. E. Phipps, E. Segré, O. Stern, Nature 
130, 892 (1932); O. R. Frisch, E. Segré, Zeits. f. Physik 80, 
610 (1933); J. M. B. Kellogg, I. I. Rabi, J. R. Zacharias, 
Phys. Rev. 50. 472 (1936). 


and at the same time strongly inhomogeneous at 
some point of the beam, are the spins ‘‘shaken off 
the field lines’’ and actual precessing occurs. Of 
course, all this may be expressed, as is usually 
done, in terms of space quantization and non- 
adiabatic transitions. 

In analogous experiments with neutrons the 
field need be neither so weak nor so strongly in- 
homogeneous, because is so much larger. 
Furthermore, it is possible, in the case of neu- 
trons, to study the especially simple case of pre- 
cession in a sharply limited homogeneous “‘pre- 
cession field’’ which is e.g. perpendicular to the 
polarizing and analyzing fields; in order to 
separate the different fields, a layer of copper 
winding carrying a suitable current may be ap- 
plied, since a thin layer of copper is almost trans- 
parent for neutrons. 


3. DEMONSTRATION OF THE PRECESSION OF 
NEUTRONS 


A. Experimental arrangement 


Our first experiment” was designed to demon- 
strate the existence of the precession and to 
measure its rate in a rough way. The arrange- 
ment (Fig. 1) consisted of a paraffin block con- 
taining a neutron source (400 mg Ra+Be), the 
polarizer, the precession field solenoid, the 
analyzer, and the detecting boron-lined ioniza- 
tion chamber connected to an amplifier and 
mechanical counter. The principle of the method 
is shown in the schematic Fig. 2. 

The polarizer was a flat iron ring (20 cm outer, 
11 cm inner diameter, 0.8 cm thickness) which 
was wound all over with one layer of 1.5 mm 
copper wire (enameled). In such a symmetrical 
arrangement the magnetizing current should 
produce practically no field outside the iron (the 
fields in the neighborhood of the single wires can 
be neglected). Actually, when a current of 10 
amp. was passed through the windings, irregular 
fields of several gauss were observed near the 
surface of the iron rings, presumably on account 
of slight local variations of the magnetic proper- 
ties. We decided therefore to switch on the cur- 
rent of 10 amp. only for a few seconds every time 
the magnetization was reversed, and then to 
switch it off again, and use the remanent field 
only. In this case the stray fields were consider- 
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© POLARIZATION EFFECT PER ceNT 


1 4 GAUSS 


44 


Fic. 3. Result of precession experiment. Polarization effect 
plotted as a function of precession field. 


ably smaller and could not have caused appre- 
ciable precession. o 

The precession field coil was made with rec- 
tangular cross section, 3.5 X5 cm, and 35 cm long, 
and consisted of a single layer of 0.65 mm copper 
wire, wound on a brass frame. The field was calcu- 
lated from the formula for an infinitely long 
solenoid. 

The analyzer was made identical with the 
polarizer. 

The boron-lined ionization chamber was of the 
“‘double-decker”’ design, in order to get a large 
useful boron surface. It was connected to a trans- 
former-coupled high gain amplifier acting on a 
mechanical counter; the resolving power of this 
arrangement was sufficient since only about 100 
neutrons per minute were counted. 

The experiment consisted in taking alternative 
counts, of 7 minutes each, with the magnetization 
of the polarizer reversed after each count; the 
analyzer was kept with the same magnetization 
all the time. Runs were made with three different 
currents through the precession field solenoid. 
Each run had to last about 50 hours in order to 
have a small statistical error. The whole experi- 
ment was carried through entirely automatically 
by means of an arrangement of relays and rotat- 


ing switches, controlled by a clock. Every 7.5 
minutes the counting was stopped during half a 
minute; within this time, a picture of the 
mechanical counter (and some control instri- 
ments) was taken on a motion-picture film, and 
the polarizer current was reversed and switched 
off again. Because of the short cycle of counting, 
slow variations of the sensitivity of the amplifier 
(which actually occurred) were practically 
eliminated. 


B. Results and discussion 


In the first run, with no current in the preces- 
sion field solenoid, the number of neutrons re- 
corded was 0.65+0.28 percent larger with the 
fields in the polarizer and analyzer parallel than 
with antiparallel fields. In the second run, witha 
precession field of 2 gauss, a difference of 0.29 
+0.36 percent was recorded, while in the third 
run, with 4 gauss precession field, the difference 
was —0.36+0.31 percent ; the negative sign indi- 
cates that the recorded number of neutrons was 
smaller with parallel than with antiparallel 
fields, in the third run. The errors have been cal- 
culated in the standard way from the square 
roots of the numbers of neutrons counted. It was 
found repeatedly that the fluctuations were not 
larger than those to be expected as purely statis- 
tical fluctuations. In this experiment, 1.2-10° 
neutrons were counted, on the whole. 

Supposing, for the moment, that all the neu- 
trons in the beam had the same velocity vo, the 
polarization effect (difference in intensity, with 
parallel and antiparallel fields) would be propor- 
tional to the cosine of the angle g=2y/I1 hi 
through which the neutrons precess on their way 
(of length 7) through the field 77. On account of 
the velocity distribution of the neutrons, the 
polarization effect as a function of // will be not 
a cosine but a superposition of cosines with 
different periods, which is roughly equivalent to 
a damping of the cosine function which corre- 
sponds to the mean velocity. We have calculated 
this function (see curve C, Fig. 5) making a 
number of assumptions which are discussed 
in §7. 

In Fig. 3 the observed differences are plotted 
together with the theoretical curve, and it is seen 
that they agree within the very large experi- 
mental errors. Obviously it would be hopeless to 
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discuss these results any further and to try and 
enclose the magnetic moment of the neutron be- 
tween definite limits. We started a new experi- 
ment with an improved arrangement, taking 
counts with five different values of the precession 
field, but from a fortnight’s continuous counting 
we found differences showing the same general 
trend but still slightly smaller than those re- 
ported here, perhaps merely because of an ad- 
verse statistical fluctuation. Substantial im- 
provement in the accuracy would have required 
such unreasonably long periods of counting that 
we decided to postpone a repetition of the ex- 
periment until a time when stronger sources 
should be available. 

The maximum polarization effect found by us 
(0.65 percent) must be multiplied by 2 in order 
to take account of the fact that about half of 
the recorded neutrons were faster than one volt 
(went through 0.5 mm of cadmium) and, there- 
fore, should not contribute to the polarization 
effect. The resulting figure, 1.3 percent, is smaller 
than the one found by Hoffman, Livingston, and 
Bethe,’ (2.3 percent). This difference is under- 
standable since we used a field of 10,000 gauss in 
polarizer and analyzer (the remanent field) while 
Hoffman, Livingston and Bethe used a field of 
15,000 gauss. 


4. SIGN OF THE MAGNETIC MOMENT OF THE 
NEUTRON 


A. Experimental arrangement 


In order to determine the sign of the magnetic 
moment of the neutron it is necessary to find the 
direction of precession. This is obviously not 
possible with the symmetrical arrangement of 
our first experiment, and the arrangement was 
therefore changed in the following way (see Fig. 
4): (1) the analyzing iron ring was rotated 
through an angle of 90°, with the beam as an 
axis, so that the polarizing and analyzing fields 
were perpendicular to each other (and to the 
neutron beam); (2) the precession field solenoid 
was removed and a short (2.2 cm long, 5 cm 
diameter) coil was inserted producing a field 
parallel to the beam; this field was found (by 
means of a small flip coil) to be fairly homogene- 
ous, presumably on account of the iron below and 


above. 


The experimental procedure this time was to 
reverse the precession field alternately while the 
rings were kept magnetized in the same direction 
all the time. The direction of precession changes 
with the precession field, and the neutron in- 
tensity recorded must be larger when the neu- 
trons precess in the same direction as one would 
have to rotate the polarizer in order to make the 
polarizing and the analyzing fields parallel (see 
Fig. 4). From the direction of the precession field 
for which the transmitted neutron intensity is 
larger than for the opposite direction, one can 
deduce the sign of the magnetic moment of the 
neutron. 

In order to get a large effect, the magnetizing 
current on the rings was left on (2 amp.) all the 
time. Consequently the stray fields were quite 
strong. By rotating the rings independently in 
their own plane, however, a position was found 
where the stray field was nearly homogeneous and 
parallel to the beam, and of the right magnitude 
(3 gauss) to turn the neutrons by 90°, on the 
average. This field could be reversed by passing 
a suitable current through the solenoid. A large 
number of counts (of 7 minutes each) were then 
taken with the solenoid current alternately on 
and off. 


Fic. 4. Determination of the sign of the magnetic mo- 
ment of the neutron, schematic diagram of the arrange- 
ment. If the neutrons precess in the way indicated, a 
larger fraction of them goes through the analyzer than 
when the current in the solenoid is reversed and the neu- 
trons correspondingly precess in the opposite way. For the 
sake of clearness, the solenoid is made much longer in 
Fig. 4 than it was in the actual experiment. 


| 


724 FRISCH, VON HALBAN, JR. AND KOCH 


B. Results 


A difference of 1.05+0.24 percent was re- 
corded between the intensities observed with the 
two directions of the precession field. The amount 
of the difference is in accordance with the ‘as- 
sumption yu, =2 n.m., although the arrangement 
is not well suited for measuring the magnitude 
of wn. From the sign of the difference it followed 
that the neutrons precess in the direction of the 
(positive) current in the solenoid; from this 
again it follows that the sign of the magnetic 
moment of the neutron is negative, that is, the 
relative direction of spin and magnetic moment 
is the same as in the (negative) electron. 

Shortly after our first publication, the same 
result was published by Powers, Carroll, Beyer 
and Dunning,’ who had used a type of an un- 
symmetrical inhomogeneous field suggested by 
Rabi."" The negative sign had been generally 
expected on account of the moments of the proton 
and the deuteron (see §1). 


5. PRECESSION OF NEUTRONS INSIDE 
MAGNETIZED IRON 


A. The problem 


Rossi!’ and Mott-Smith'® have tried to detect 
the deflection which cosmic-ray electrons undergo 
when passing through magnetized iron. For the 
interpretation of such experiments it is essential 
to know, whether the amount of deflection corre- 
sponds to the magnetic induction B (about 
20,000 gauss at saturation) or the magnetic field 
strength J/7 (in general only a few gauss). 
Weizsicker®® has examined the question with 
the help of Dirac’s theory and found that the 
induction B is effective. 

We have put to ourselves the analogous 
question: Is the rate of precession of neutrons 
inside magnetized iron determined by // or B? 
So far, no theoretical investigation of the prob- 
lem as such has been published ; but Bloch*! made 
some interesting remarks from which it appears 
that the theoretical treatment of this and some 


8B. Rossi, Rend. Lincei 11, 478 (1930); Nature 128, 
300 (1931). 

1%L. M. Mott-Smith, Phys. Rev. 37, 1001 (1931); 39, 
403 (1932). 

20 C.F. v. Weizsiicker, Ann. d. Physik 17, 869 (1933). 

21F, Bloch, Phys. Rev. 51, 994 (1937). 


similar problems would require a much more 
intimate knowledge of the interaction between 
electrons and neutrons than is available at pres- 
ent. In this state of affairs we thought that even a 
crude experiment might be of interest. 


B. Experimental arrangement 


The arrangement was in principle similar to 
our first one, but yet different in several respects. 
As polarizer and analyzer we used straight iron 
bars of cross section 1X5 cm, and 80 cm length, 
carrying a single layer of 0.65 mm copper wire; 
by means of iron pieces at the ends they were 
connected to form a single magnetic circuit. The 
magnetic induction in the bars was about 14,000 
gauss. The space between them was investigated 
with a small flip coil, and the neutron beam was 
sent through at a point where the stray field was 
weak enough to cause no disturbance. 

A long, flat solenoid (7 X0.7 cm cross section, 
15 cm long) was placed across the beam, similar 
to the first arrangement, and a strip of thin (0.15 
mm) iron sheet inside the coil; the ends of the 
strip were connected by a strong iron yoke (made 
of round iron of 1 cm diameter) to close the 
magnetic circuit. The current through the 
solenoid was alternately set at two different 
values corresponding to magnetic fields of 2.8 
gauss and 35 gauss, respectively ; the polarizing 
and analyzing fields remained unchanged (they 
were bound to be antiparallel since they formed 
a magnetic circuit). 

In order to enhance the polarization effect we 
used a low temperature source of neutrons,'° 
consisting of an ice block at liquid-air tempera- 
ture, of similar shape as the paraffin block used 
before. 


C. Results 


To check the whole arrangement, a run was 
made with the iron strip removed. With a preces- 
sion field of 2.8 gauss, along a path of only 7 mm, 
the neutrons should hardly precess at all and the 
transmission of the system should have the low 
value corresponding to the polarizing and analyz- 
ing fields being antiparallel. With the precession 
field at 35 gauss, however, most of the neutrons 
precess by more than 360° and the beam should 
be practically depolarized, and consequently the 
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Fic. 5. Curves representing the precession of 
thermal neutrons. 


transmission of the system should lie halfway be- 
tween the values for parallel and antiparallel 
fields. Such a difference was actually observed ; 
the intensity was 0.83+0.26 percent larger with 
35 gauss than with 2.8 gauss. This figure of 0.83 
is to be compared with one-half of 0.65, the result 
in our first experiment (see §3B). It should be 
remembered that in this first experiment the 
polarizing field was alternately reversed. The 
improvement by a factor of 2 or 3 is to be at- 
tributed partly to the cooling, and partly to the 
larger magnetic induction in polarizer and ana- 
lyzer (14,000 instead of 10,000 gauss). 

The iron strip was now inserted (and closed by 
the yoke) and a second run, quite similar to the 
first one, was carried out. If the rate of precession 
in the iron was determined by the field strength 
H, then the iron strip should make no difference 
at all, because the field strength is the same within 
the iron as below and above it. If it is B that 
matters, however, then the neutrons should 
precess many times on their way through the iron 
sheet and be completely depolarized; since the 
iron is practically saturated even with the weak 
field of 2.8 gauss, the neutrons reaching the 
analyzer should be depolarized and the intensities 
transmitted the same both for 2.8 and 35 gauss. 
This is actually what we found: the difference 
was —0.05+0.24 percent, that is, zero within 
the limits of error. 

It is concluded that the field accounting for the 
rate of precession of neutrons inside magnetized 
iron is certainly much larger than // ; a numerical 
discussion shows that it is probably larger than 
500 gauss, or, in other words, larger than 0.03- B. 


6. CALCULATION OF THE PRECESSION OF 
NEUTRONS WITH MAXWELLIAN 
VELOCITY DISTRIBUTION 


We assume the velocity distribution by the 
neutron beam and the velocity dependence of the 
indicating system to be such that f(v)dv is the 
number of neutrons with velocities between v and 
v+dv counted per unit time, with polarizer and 
analyzer removed. With the polarizer (of thick- 
nesg x) in place (and magnetized) the number of 
neutrons is reduced to 


Jo(v)dv = 3 


where yu; and ye are the attenuation coefficients 
for parallel and antiparallel neutrons, re- 
spectively. 

If now the analyzer (also of thickness x) is 
brought in place, the transmitted intensity de- 
pends on its position : 


parallel J\(v)dv= bf | (1) 
antiparallel (2) 
Introducing and po(v)=~—B(v) 


(we may assume that uw does not depend on v) we ° 
find 


[Ji(v) —J2(v) Jdv 
J. (3) 


If the analyzer forms an angle ¢ with the 
direction of polarization, (3) must be multiplied 
with cos ¢. We assume that the analyzer forms 
an angle a with the polarizer, and that the 
neutrons pass through a magnetic field // (per- 
pendicular to both polarizer and analyzer) and 
precess through an angle ¢,=2y,///] hv. Then we 
have ¢ = ¢,—a, and hence, if we integrate over v 


f 
-[cos a cos ¢,(v)+sin @ sin ¢,(v) ] 


=const. cos a f s(o)do-sinh? [x-B(v) ]-cos ¢,(v) 


+sin a f s(o)do-sinh? [x-B(v) ]-sin 


=const. [cos a:-C+sin a: S]. 
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Fic. 6. Curve representing the precession of thermal 
neutrons. The distance OA indicates the polarization effect 
obtained for a given value (e.g., 0.4) of ger and an angle a 
between polarizer and analyzer. 


The functions C and S have been calculated 
(by numerical integration) and are plotted in 
Fig. 5. The abscissa is 
the angle through which neutrons with an energy 
of kT precess on their way / through the field 77. 
For neutrons with a magnetic moment of 2 n.m. 
and thermal energy distribution at room tem- 
perature we have ¢xr=0.088/H. The scale of 
ordinates is chosen so that C(O) =1. 

In calculating the functions C and S, we have 
made the following assumptions: 

(1) f(v)dv=const. which means 
Maxwellian distribution of the neutrons and 1/v 
sensitivity of the boron chamber. The integra- 
tions were carried from zero to infinity. 

(2) sinh [x8(v) ]~x8(v), which is correct for the 
larger part of the velocity spectrum where 
xB(v)<1 (with x~1 cm). 

(3) B(v) =const./v. This assumption’ should be 
correct for the mean part of the spectrum, but 
tends to overemphasize the importance of the 
slowest neutrons. We have therefore recalculated 
a part of the function C for neutrons of room 
temperature, using for 8(v) the more accurate ex- 
pression’ 


const. [K+1—(2K+1)!]/K?(K 


It is seen (Fig. 5, curve C’) that the initial slope 
is less steep, but the intersection with the axis 
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of abscissae is not considerably shifted. It should 
be remembered, however, that the expression for 
B(v) depends on the density distribution of those 
electrons which are responsible for the ferro- 
magnetism of the polarizer material. Since this 
density distribution is not too well known, any 
determination of u, from the precession of neu- 
trons with a wide velocity band would be some- 
what uncertain. For high accuracy, experiments 
with neutrons of homogeneous velocity will be 
required. 

In Fig. 6 a somewhat different representation 
of the functions C and S is given. Corresponding 
values of C and S are plotted as abscissae and 
ordinates, respectively, and the parameter ¢gyr is 
noted along the curve. The polarization effect for 
any angle a between polarizer and analyzer and 
for any value of g;7 may then be obtained in the 
way indicated in the figure. For a given value of 
gir the largest polarization effect is obtained 
when the analyzer is placed in the direction of 
the vector R (a=¢*). The vector R represents 
therefore the “‘mean spin direction”’ in the beam, 
and its length is proportional to the polarization 
effect. If all the neutrons in the beam had the 
same velocity v=v,.7=(2kT/m)!, then ¢* would 
be equal to g:r and the corresponding vector R 
would rotate at a uniform rate without changing 
in length, and would represent directly the pre- 
cession of the neutrons. The velocity distribution 
of the neutrons has the effect that, firstly, the 


amount of R decreases with increasing 


(“‘depolarization”’) and, secondly, ¢* is not pro- 
portional to g,r, but increases more slowly at 
large values of gr; this is because the slowest 
neutrons are the first to be depolarized, and the 
remaining faster ones require stronger fields for 
precessing through the same angle. 

In conclusion we wish to thank Professor Dr. 
N. Bohr for his kind interest and for putting the 
experimental resources of the Institute of 
Theoretical Physics at our disposal. For many 
stimulating discussions we wish to thank Drs. 
C. Mgller, G. Placzek and V. Weisskopf. 
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The Isotopic Constitution of Gadolinium, Dysprosium, Erbium and Ytterbium 


A. J. DEMPSTER 
Ryerson Physical Laboratory, University of Chicago, Chicago, Illinois 
(Received February, 24, 1938) 


The mass spectra of four rare earth elements have been photographed, by the use of a 
spark source for the ions. New faint isotopes have been found in all four. The mass numbers 
of the new isotopes and rough estimates of their abundance in percent are: Gd, 154 (1.5), 
152 (0.2); Dy, 160 (1.5), 158 (0.1); Er, 164 (2), 162 (0.25); Yb, 170 (2), 168 (0.06). 


HE isotopic constitution of all the rare 

earth elements has been examined by Dr. 
Aston! with anode rays.2 The author has also 
photographed the mass spectra of cerium*® and 
neodymium and found in each of these elements 
two new faint isotopes that had not been ob- 
served by Aston. In samarium‘ no new isotopes 
were found. 

The source of the neodymium ions was a 
spark between electrodes of the metal itself, and 
for cerium ions the spark was passed between 
electrodes made by packing nickel tubes with a 
mixture of the oxide and magnesium and 
aluminium powder. It was suggested by Hénig- 
schmid® that the existence of two new isotopes 
at 148 and 150 might remove the discrepancy 
between his value of the atomic weight of 
neodymium 144.27 and Aston’s mean value 
143.5, but a measurement of their intensity by 
Mattauch and Hauk® showed they could account 
for only a part of the discrepancy. The remainder 
is removed if we take new values determined by 
them for the abundance of the strong components 
in place of those obtained by Dr. Aston. 

Recently the attention of the author has been 
called to the possible existence of faint isotopes 
in other rare earth elements by Dr. Pool. His 
experiments on the radioactivity induced in these 
elements by fast neutrons’ show that at least 
three new stable rare earth nuclei should exist, 
and it was thought desirable to examine the even 
numbered rare earth elements for possible faint 


1 F, W. Aston, Phil. Mag. 47, 385 (1924); 49, 1191 (1925). 
“” W. Aston, Mass Spectra and Isotopes, (1933) p. 65, 

8A, J. Dempster, Phys. Rev. 49, 947 (1936). 

*A. J. Dempster, Phys. Rev. 51, 289 (1937). 

°O. Hénigschmid, Naturwiss. 25, 701 (1937). 

6 J. Mattauch and V. Hauk, Naturwiss. 25, 781 (1937). 

7M. L. Pool and L. L. Quill, Phys. Rev. 53, 437 (1938). 


isotopes. The stable atomic nuclei now appear 
to be bounded by radioactive forms, and it is 
thus of greater interest to complete the list of 
stable isotopes. 

The mass spectrograph has been previously 
described.*: * Dr. Pool kindly supplied a few 
milligrams of the oxides of the three elements 
Gd, Er, Yb from the samples used by him in 
his experiments. The electrodes were made by 
packing nickel tubes about one millimeter in 
diameter with the oxide mixed with a reducing 
agent. In these experiments, the reducing agent 
was either calcium, lanthanum or neodymium 
metal in the form of small filings. The heats of 
combination’® are, for LagO3, 457 calories; for 
Nd.O;, 435 calories; for 3CaO, 456 calories. 
For the other rare earths, the values are ap- 
parently unknown. It was found that with 
lanthanum or neodymium as reducers, the 
spectra were several times as intense as with 
calcium. It was found that the oxides of the 
reducing agents appeared faintly, and with Gd, 
Dy, and Er, the oxides of neodymium could 
have been confused with their isotopes. In these 
cases, lanthanum or calcium was also used to 
reduce the oxides. 


GADOLINIUM 


Gadolinium oxide was first reduced with 
neodymium metal. A new isotope at mass 154 
was observed on four photographs with exposures 
of ten to seventy minutes, and an isotope at 152 
on two plates with seventy minutes exposure. 
A second series was made in which the oxide 
was reduced by lanthanum metal. Here the mass 


8 A. J. Dempster, Proc. Am. Phil. Soc. 75, 762 (1935). 
®A. J. Dempster, Rev. Sci. Inst. 7, 46 (1936). 
1° Bichowsky and Rossini, Thermochemistry (1936). 


727 


728 A. J. DEMPSTER 


155 


162 164 166 168 170 
te 
44 


Fic. 1. Mass spectrum of erbium showing two new isotopes at 164 and 162. The mass at 169 is ascribed to thulium 
present as an impurity (about } percent). The mass at 155 is due to lanthanum oxide. The lines due to the strong com- 


ponents are broadened because of excessive over-exposure. 
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Fic. 2. Mass spectrum of ytterbium showing two new isotopes at 170 and 168. The lighter masses 166, 164, 162, etc. 
are due to neodymium oxide. 


at 154 was observed on thirteen photographs 
with exposures ranging from half a minute to 
twenty minutes and the mass at 152 was ob- 
served on four with exposures of six to twenty 
minutes. From the series of photographs with 
different exposure times, it was estimated that 
the abundance of the new isotope at 154 was 
approximately 1.5 percent of the total and that 
the isotope at 152 was 0.2 percent of the total. 
The photographs also showed two europium 
lines at 151 and 153 with about the same in- 
tensity as the weakest gadolinium isotope at 
152. 
DyspROsIUM 

Two different samples of dysprosium oxide 
were used and reduced with lanthanum. In 
addition, one was reduced with calcium. Both 
samples showed two new isotopes at masses 160 
and 158. One sample also showed strong lines at 
165 and 159 due to holmium and terbium. With 
the second sample, these lines were much weaker. 
The isotope at 160 was found on 12 photographs 
with exposures from one-half to forty minutes, 
the one at 158 on 4 photographs. From the 
exposure times, it was estimated that the isotope 
at 160 makes up about 1.5 percent of the total, 
and the one at 158, 0.1 percent. In the second 
sample which was kindly supplied by Professor 
Quill, the holmium was estimated at 0.4 percent 
and the terbium at 0.2 percent. No line was 
observed at 166. 


ERBIUM 


Two new isotopes were also observed in erbium 
reduced with lanthanum at masses 164 and 162, 
the first on eleven photographs with exposures 
of ten seconds to twenty minutes and the second 
on four photographs with seven to twenty 
minutes’ exposure. An example of the mass 
spectrum is given in Fig. 1. The abundances 
were estimated as approximately 2 percent for 
the mass at 164 and 0.25 percent for the mass at 
162. The spectra also show a faint line at 169 
(0.5 percent) which was ascribed to thulium. 
No other isotopes were observed as strong as 0.03 
percent. 

YTTERBIUM 


The ytterbium oxide was reduced by neo- 
dymium and also showed two new isotopes at 
masses 170 and 168. The heaviest of the faint 
neodymium oxides came at 166. A spectrum is 
shown in Fig. 2. The mass at 170 was found on 
eight exposures of 20 seconds to 80 minutes, 
and the isotope at mass 168 on two exposures of 
60 and 80 minutes. The abundances were esti- 
mated at 2.0 percent and 0.06 percent. Faint 
lines at 175 and 169 were ascribed to lutecium 
and thulium (<0.03 percent). 

The thanks of the author are due to Professor 
M. L. Pool and Professor L. L. Quill for kindly 
supplying samples of the rare earth elements for 
use in the investigation, and to Dr. A. E. Shaw 
for assistance in preparing electrodes. 
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On the Sensitivity of the Balanced Space Charge Method for Detecting Ionization of 
Gases by Collision of Ions and Atoms 


ANTONIO ROSTAGNI 
Istituto di Fisica, R. Universita di Messina, Italy 


(Received November 30, 1937) 


Factors which influence the balanced space charge method of detecting ionization are dis- 
cussed. These factors are (1) the velocities of the ions, (2) the distances traveled by the ions 
before they are neutralized, (3) the effect of secondary electrons on the space charge. It is 
concluded that the reduced sensitivity in some experiments may account for the failure to 


observe ionization. 


NSET potentials for the ionization of gases 
by positive ions and atoms have been 
investigated in several papers by R. N. Varney,' 
who used a balanced space charge method. Last 
year he reported experiments on noble gases 
bombarded by their own accelerated atoms. His 
results for A and Ne agree quite well with those 
obtained previously by me.? With He atoms in 
He he does not observe any ionization up to 400 
volts primary energy, while I have found a 
positive effect between 60 and 460 volts, and 
have measured the ionization efficiency, which is 
about ten percent that for A in A.’ 

Varney estimates, merely from the size of 
galvanometer readings, that his method would 
have revealed an ionization effect five percent 
that for A in A. For this reason and on account 
of the similitude of A and He curves, he is 
inclined to regard my results on He as having 
been due to some A impurities. Varney must 
implicitly have made some assumption about the 
relative effects of He and A on space charge, but 
on this very essential point no specification is 
given in his paper, and I shall show that, whatever 
the assumption may be, it must necessarily 
contain much that is arbitrary. 


'R. N. Varney, Phys. Rev. 47, 483 (1935); 50, 159 
(1936); 50, 1095 (1936). Varney, Gardner and Cole, Phys. 
Rev. 52, 526 (1937). 

2A. Rostagni, Nuovo Cimento 11, 621 (1934); 13, 389 
(1936); Ricerca Scient. 7, 511 (1936). 

* My lowest measured point for He corresponds in the 
published diagram to 62 volts on the abscissa, and 0.01 
cm?/cm* on the ordinate (cross section for ionization); 
which represents, as explicitly stated in my paper, the 
sensitivity limit of the method. Direct extrapolation from 
the traced curves would point to an onset potential ranging 
between 50 and 60 volts. I therefore question Varney’s 
statement that my results indicated ionization energies 
less than twice the electron ionization potential of He. 


Neither does Varney stop to consider the a 
priori rather surprising fact that helium will 
behave so differently from argon and neon as his 
experiments indicate. We may remark by the 
way that a similar anomaly is shown also by 
previous experiments with the same method ; for 
example Varney‘ does not observe any ionization 
of noble gases by the lightest alkali ions Lit and 
Na*, with the exception of Ne by Na*, while he 
has positive results with any other alkali ion 
noble gas combination. 

It appears a legitimate question to ask whether 
the sensitivity of the space charge method has not 
been overestimated, in respect to its adaptation 
to the problem of ionization by massive particles. 
Investigations of this matter have not been 
published. I think they are necessary, because, as 
we shall see, the mechanism involved is quite 
different from all previous applications of the 
method. 

The enormous efficiencies of positive ions in 
neutralizing electron space charge observed by 
Kingdon,*® Foote and Mohler® and others are the 
result of two factors: 

(1) The velocities of the ions are, for the same 
kinetic energy, i.e., for a given. potential differ- 
ence, some 10° times smaller than those of 
electrons (in inverse proportion to the square 
root of masses). 

(2) Positive ions produced inside the diode 
withasmall initial velocity component transverse 
to the filament must travel many times to and fro 


4R. N. Varney, Phys. Rev. 47, 483 (1935). This partially 
negative result contradicts also Beeck and Mouzon's ex- 
periments by another method: Beeck and Mouzon, Ann. 
d. Physik 11, 737, 858 (1931). 

5 Kingdon, Phys. Rev. 21, 408 (1923). 

® Foote and Mohler, Phys. Rev. 26, 195 (1925). 
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across the inner space before they discharge 
themselves. 

Under Kingdon’s experimental conditions the 
average distance which an ion travels before 
reaching the filament is of the same order of 
magnitude as the kinetic mean free path. This 
may be some 10? times the cathode-anode 
distance at sufficiently low pressure. We thus 
obtain amplification factors of 10* to 10° for 
electron current with regard to positive ions as 
observed by different authors. Now the experi- 
mental conditions of Varney seem to me not very 
different from those of Kingdon, as regards the 
discriminating factor R)*; but the pressure 
being of 10°? mm of Hg, we should expect 
amplification factors of 10* both for A+ and He* 
ions, with a ratio of about 3 in favor of A. 
Varney gives much greater values; it is possible 
that he has determined them directly, by some 
method equivalent to those used by Kingdon, 
although no mention of that is made in his papers. 

But there is an essential remark to be made. 
A necessary condition in order that factor (2) may 
be effective at all is that the initial energies of the 
ions should be less than the applied potential, so 
that if they will not hit the cathode they will be 
also unable to reach the anode. That is of course 
verified with ions produced by photoelectric 
effect or by electron collision as in the experi- 
ments of Foote and Mohler, Kingdon, etc ; but in 
the present case, where the masses of colliding 
particles are of the same order of magnitude as 
the gas molecules, we must expect most of the 
secondary ions to receive a finite fraction of 
primary velocities. That is shown directly by my 
own experiments. And as the applied potential 
ranges between 1 and 2 volts (it must be less than 
the resonance potential of the gas) most of the 
ions will therefore reach the anode. Factor 
(2) then becomes inefficient. As the energy of 
many of the ions is large compared with that of 
electrons, also factor (1) will be affected. 

Now the sensitivity of the space charge method 
is actually greater than it would appear from this 
remark, or the method would not have yielded 
any valuable result. The sensitivity may only 
arise from the partial neutralization of the 
secondary ions on the way to the anode, giving 
rise to slow tertiary ions for which both sensitivity 
factors may be efficient. This fact seems not to 


have been taken into consideration: the neu- 
tralization plays an essential part in the working 
of the ion detector and the mechanism is more 
complicated than in usual applications.’ The 
efficiency must depend on the following quanti- 
ties: (1) cross sections for neutralization; (2) 
cross sections for scattering ; (3) velocities of slow 
ions arising from neutralization ; all as functions 
of the velocity of secondary and eventually of 
tertiary ions. These quantities should be known 
for every ion velocity between primary and 
thermal velocity, but they are not completely 
known :* and in any case the problem will be 
much too complicated to allow a calculation 
either of absolute or of relative efficiency. To 
determine it by experiment, the ions under 
question must be produced under proper 
conditions. 

Another effect which has been totally neglected 
in the discussion because of the relatively low 
sensitivity of the ion detector to electrons, is the 
secondary electron emission from metal surfaces. 
As it has been shown that the amplification factor 
measuring the ratio of sensitivities may be less 
than expected it seems advisable to examine 
whether secondary electrons may not sometimes 
influence the measurements and perhaps in 
extreme cases counterbalance the effect of posi- 
tive ions on space charge, thus putting a practical 
limitation to the efficiency of the method. 

The principal source of secondary electrons in 
the space charge cylinder would be of course the 
metal surface hit by the primary beam in front of 
the entrance hole. As a discriminating factor to 
the relative efficiency of the method in the 
investigation of various colliding particles in 
differents gases, we may of course adopt the 
quotient qg; k of the cross section for ionization of 
the gas by the coefficient of electron liberation 
from the metal surface by the particles under 
question. These quotients have been measured by 
me®* for both systems A gas A atoms, He gas He 


7 The part of neutralization in the working of the ion 
detector also in previous applications is probably more 
important than hitherto supposed, and ought to be taken 
into account for a complete theory of it; however it is 
not so essential as here. 

8’ A. Rostagni, Nuovo Cimento 12, 134 (1935); Atti di 
Torino 70, 472 (1935). F. Wolf, Ann. d. Physik 29, 33 
(1935). For anomalies presented by both cross sections for 
ion velocities under 30 volts see also A. Rostagni, Atti del 
Congresso Galvani (Bologna, October, 1937) in press. 

® A. Rostagni, Nuovo Cimento 11, 621 (1934). 
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atoms and copper surfaces. For energies of the 
particles between 60 and 460 volts they are 
always about fifty times larger with A than 
with He. 

We are forced to conclude that if secondary 
electron emission really can exert any limiting 
effect on the sensitivity of the method, that will 
be much stronger in He than in A. 

In order to calculate the magnitude of the 
effect we may consider a particular case. With 
200 volt He atoms in He at a pressure of 10-? mm 
Hg qi k& is 0.006 ions per electron-cm. Since the 
beam length in the space charge cylinder was 2 
cm, and since the effective cross section for 
collision of He atoms of 200 volts energy in He 
gas at 10~* mm is less than 0.1 cm, it follows 
that more than 80 percent of the primary atoms 
will hit the top of the cylinder with full energy, 
and liberate there about 70 electrons per positive 
ion produced in the gas. As the wall is positive 
with respect to inner space, not all the electrons 
will be able to travel far from it, but, since the 
field is small and concentrated near the filament, 
there is no doubt that a large fraction of them 
will contribute to increase the negative space 
charge. 

Besides, in the special arrangement for in- 
vestigating ionization by atoms,'® we may look 
for another possible source of secondary electrons 
at the back of aperture 3, where practically all 
primary ions stopped by the field between 3 and 4 
fall with their full energy. Unless the cylinder is 
negative with respect to aperture 3 (no mention 
of that is made in Varney’s paper), a large 
fraction of these secondaries will be conveyed by 
the same field, through aperture 4, into the 
cylinder. 

The proportion of secondary electrons to the 


1” See Fig. 1, in Phys. Rev. 50, 159 (1935). 


ions inside the space charge cylinder should then 
be of the order of 100 to 1, in the case of helium. 
This result seems to justify our doubts on the 
legitimacy of neglect of secondary electrons. 
Moreover, the coefficient & for He atoms is 
rapidly increasing between 50 and 100 volt (from 
0.01 to 0.1), at a rate which is comparable with 
that of ionization ;" therefore it seems possible 
that the augmentation of secondary emission has 
masked the onset of ionization. Above 100 volts 
the quotient g; k for He remains practically 
constant. Any further increase of the velocities 
of the atoms therefore will be of no use in making 
ionization detectable. 

By way of conclusion: while the positive 
results of the balanced space charge method in 
determining onset potential are certainly the 
most reliable, in order to secure reliability of 
negative results the sensitivity of the method 
needs to be more thoroughly investigated. As 
regards the particular case of the ionization of 
helium I am of the opinion that Varney’s 
published results do not prove anything against 
the consistency of mine. These are confirmed on 
the one side by independent measurements with 
two quite different arrangements, on ions as well 
as on atoms, and by the perfect agreement of 
part of them with Varney’s positive results, on 
the other side. I do not see any specific reason why 
if my results on A and Ne are correct, those on 
He should be completely wrong, since I reject the 
too simplifying hypothesis that my He had been 
contaminated, and since the principal disturbing 
effect of collecting field methods, the secondary 
electron emission, has been directly accounted 
for. | have, however, pointed out two independent 
reasons why Varney’s result on He could be 
inconsistent. 


"A. Rostagni, Nuovo Cimento 11, 99 (1934); Zeits. f. 
Physik 88, 55 (1934). 
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Additional information is given on the reaction of the balanced space charge detector to 
fast positive ions in the presence of gas in the tube. The sensitivity of the detector is analyzed, 
primarily in terms of experimental observations. Further notes on the mechanism of the de- 


tector are added. 


N view of the questions raised by Rostagni! 
on the sensitivity of the balanced space charge 
method of observing ionization of gases, used by 
the writer in several recent experiments,” * it is 
necessary to describe more in detail the observa- 
tions which led to the estimated sensitivities 
reported for the method. 

The basic features of the detector may be 
summarized : 

1. A hot filament cathode emits a space charge 
limited electron current to an almost completely 
enclosing cylindrical anode. 

2. A positive ion in the space charge region has 
the effect of neutralizing part of the space charge 
and thus allowing an increased electron flow to 
the anode. 

3. This increased flow is measured in the 
balanced space charge method by a sensitive 
galvanometer connected in a balanced bridge 
circuit. 

4. The size of the increased electron current 
per positive ion depends on the mean length of 
time a positive ion remains in the space compared 
with the same time for an electron. 

5. The time a positive ion spends in the space 
depends on: 


(a) Slow movement of the ion. 
(b) A long free path for the ion. 


Rostagni points out very accurately that 
neither of the conditions 5 (a) or 5 (b) are 
satisfied in the writer’s experiments. (a) The 
positive ions all had high speeds. (b) The gas 
pressure was so high that free paths were short. 

On the other hand, the writer did observe !arge 
effects by positive ions on the electron current, 

1 A. Rostagni, Phys. Rev. 53; 729 (1938). _ 


2? R. Varney, Phys. Rev. 47, 483 (1935). 
3R. Varney, Phys. Rev. 50, 159 (1936). 


and this paper describes some experiments which 
indicate how Rostagni’s proposed mechanisms 
actually operate. 

The direct effect on the space charge electron 
current produced by a beam of positive alkali 
ions of various speeds was first observed. Gas 
pressures (A, Ne, No, He) ranging from less 
than 10-° up to 10-? mm of Hg were used. The 
results are summarized here: 

1. From 0 to 10 volts ion speed, the ordinary 
positive ion effect on the space charge was 
observed. As predicted, it fell off rapidly with 
increasing speed. 

2. Above 10 volts, the beam had a very nearly 
constant effect on the space charge current. 

3. The size of the effect in 2 depended on the 
gas pressure. In vacuum the effect was very small. 
With gas pressures from 5 X10-* to 10-? mm, the 
increase in electron current reached 250 to 300 
times the actual positive ion current. 

4. Helium pressures had to be higher than 
argon pressures to produce a given size of effect. 

5. No conspicuous differences were observed 
among Lit, Nat, K+, Rb*, and Cs*. 

It is concluded from these experiments that 
scattering of ions in the gas plays an important 
role. The exchange of charge between alkali ions 
and atoms of the gases used is very small so can 
hardly influence the results. That scattering of 
fast ions from the anode cylinder walls may under 
favorable conditions influence the space charge 
current considerably is indicated in recent results 
of R. Kienzle* who tried similar experiments in 
vacuum. 

Further analysis of the detector resulted from 
the failure to observe ionization in some cases, 
chiefly neon and helium, where Beeck and 


*R. Kienzle, Ann. d. Physik 30, 401 (1937). 
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Mouzon® and Rostagni® were more successful. 
The writer attempted to set an upper limit on the 
jonization efficiencies which might escape de- 
tection by the balanced space charge method. A 
discussion of the effect in helium as compared 
with argon is given, as an example, below. 

About 1/150 of the deflection produced by 
jonization of argon by 106 volt argon atoms could 
have been observed. The published figure, 
namely, that 1, 20 (or 5 percent) of the ionization 
in argon could have been observed was based on 
the following : 

(a) Helium ions of a given energy have over 3 
times the speed of argon ions of the same energy. 
Since slow ions are more effective in space charge 
neutralization, this factor of 3 goes against the 
sensitivity to helium. 

(b) It is more difficult to produce a beam of 
helium neutral atoms, as the charge exchange 
cross section for (He* — He) is much less than for 
(A+—A). It is the opinion of the writer that this 
factor far outweighs any considerations of space 
charge insensitivity. It is inapplicable to the 
results on ionization by positive ions, of course. 
To offset the difficulty, the initial intensity of the 
helium beam was enlarged by increasing the 
intensity of the preneutralization ion beam* by a 
factor of 2 of 3. The pressure of helium in the 
apparatus was varied over a range from 5X 10-* 
to 5X10-? mm of Hg. (This figure was unfortu- 
nately omitted in the original report.) No trace of 
ionization was ever observed. The combination of 
factors (a) and (b) was set at 7.5 as a result. The 
sensitivity of the detector to helium must be 
reduced by this factor in comparison with argon. 
1/150 of the deflection in argon could have been 
observed, so in helium 7.5 times this amount or 
1/20 of the deflection in argon should have been 
detectable. 

It may be pointed out here that Rostagni! has 
given an argument which the writer overlooked 
showing that helium ions will be more effective 
than argon ions in one respect. Since the collision 
cross section of argon is much greater than of 
helium (factor of at least 3), helium ions will 
travel much farther between collisions than 
argon ions. The sensitivity of the space charge 

5O. Beeck and C. Mouzon, Ann. d. Physik 11, 737, 858 


(1931). 
§ A. Rostagni, Nuovo Cimento 11, 34 (1934). 


detector depends on this distance between 
impacts, so a factor of 3 or more in favor of the 
sensitivity to helium results. 

One experimental result, however, is probably 
more significant than all the complex analysis. 
The result is the following: The relative effi- 
ciencies for ionization of Ne by Nat, A by K*, 
etc., as observed, uncorrected, by the space 
charge method were of the same order of magni- 
tude as Beeck’ had observed by an absolute 
method. In fact they agreed as closely as 10 
percent. No attempt at accurate comparisons of 
ionizations was made with the balanced space 
charge method, but the ion beam intensity and 
the gas pressures were kept practically constant 
for all ion sources. The galvanometer deflections 
in the balanced space charge circuit were thus 
roughly an indication of the combined efficiency 
of ionization and efficiency of detection. The 
close agreement with Beeck’s measurements of 
ionization efficiency led to the conclusion that the 
efficiency of detection must have been nearly 
constant for different gases all the way from neon 
to xenon. The decreased sensitivity factor of 7.5 
in going from argon to helium was thus believed 
to be by no means an underestimate. 

Rostagni! has suggested a masking of the 
ionization in helium by secondary electrons. The 
use of more than one gas pressure should auto- 
matically eliminate this possible error, for a 
surface and a gas phenomenon can only coincide 
at one particular gas pressure. It would further- 
more require that the coefficient of ionization and 
the coefficient of electron liberation be identical 
functions of voltage, and Rostagni’s® results do 
not indicate this to be true in the important 
region near the onset voltage. 

An error particularly to be feared was that 
secondary electrons ejected from the back of a 
slit might in some manner be accelerated into the 
detector and actually ionize the gas. We ruled 
this out quite certainly by intentionally placing 
an electron accelerating voltage between the last 
slit and the detector. It not only did not -alter 
previous results, but ionization of such gases as 
He, Ne, and Hg still did not appear. This same 
technique was used in all experiments including 
those on ionization by neutral atoms. Rostagni’s! 


suggestion that secondary electrons might enter 


70. Beeck, Ann. d. Physik 6, 1001 (1930). 
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from the back of the last slit seems to be definitely 
outlawed. The effect could at best be only of the 
second order of small quantities. 

Just what mechanisms are involved in the 
space charge detector containing gas and used for 
fast ions is still difficult to conclude. Scattering 
from the walls and from gas atoms has already 
been advanced. I am indebted however to 
Professor L. B. Loeb*® for an additional and 
important suggestion concerning the operation of 
the detector. If the filament of the detector is 
run at a very high temperature and the applied 
voltage is small, the space charge may be so 
heavy that an actual positive ion trap or potential 
trough of appreciable size exists between the 
filament and cylinder. If as the result of collisions 
or any other process, a positive ion finds itself in 
this trough with insufficient energy to escape, it 
may be extremely effective in neutralizing 
electron space charges. That this does occur is 
indicated by the frequent observations of the 
writer of an apparent increase in sensitivity of 
considerable magnitude with increase of filament 
temperature. Neither the space charge limited 
current nor the applied voltage was allowed to 
vary as the temperature was raised. The filament 
was already far into the space charge controlled 


8 Private correspondence. 
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region of emission at the start. The slight 
increase in potential drop down the filament was 
readily shown to be inadequate explanation of 
the sensitivity increase, for the plain variation of 
the applied potential gave no such effect. 

Rostagni! has challenged the writer’s negative 
results in some cases on the grounds that such 
results are surprising after all the positive ones 
with other noble gases. Without attempting to 
discuss what should or should not be expected, 
the writer feels obliged to point out that his 
statement that ionization efficiency in helium 
must be less than 5 percent of that in argon is 
hardly a denial of all ionization of helium by 
helium atoms. Rostagni’s challenge could possibly 
be directed at the smallness of the result but not 
at its absence. 

It has been pointed out previously, further- 
more, that Beeck’s data on the onset of ionization 
by positive ions fail to be even self-consistent in 
the cases (and only in those) where the writer has 
found no ionization by the balanced space charge 
method. Rostagni’s® brilliant work on ionization 
by neutral atoms is clouded by a similar difficulty ; 
the onset potential has been observed in every 
case to rise as the atomic weight decreases and as 
the ionization efficiency decreases. His result on 
He — He is the only one incompatible with thislaw. 
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Quadrupole-Quadrupole Interatomic Forces 
Jutian K. Knipp 
Laboratory of Physics, University of Wisconsin, Madison, Wisconsin 
(Received March 15, 1938) 


Atoms the normal states of which are not S states ex- 
perience quadrupole-quadrupole repulsive and attractive 
forces which vary as the inverse fifth power of the separa- 
tion distance. At large distances, exchange being negligible, 
the expression for the interaction energy contains as factors 
two atomic coefficients and a root of a secular equation 
determined by the molecular state of the system. For a 
given atom the coefficient is determined by the nature of 
its lowest level. If one uses the Hartree-Fock approxinia- 
tion, it is proportional to the average of the square of the 
radial atomic distance for those electrons with orbital 


HEN two atoms are far enough apart 
that overlapping and consequent exchange 


angular momentum not in complete shells. Atomic coeffi- 
cients have been calculated for most of the atoms of the 
periodic table having incomplete p and d shells and the 
secular equations have been solved for a number of cases. 
At distances of twice the sum of the atomic radii, diatomic 
molecules resulting from the combination of such atoms of 
the first row of the periodic table have quadrupole-quadru- 
pole energies of a few tenths of a volt, such energies be- 
coming rapidly larger for smaller distances because of the 
inverse fifth power dependence. 


effects are negligible, the interaction can be con- 
sidered as being composed of mixed terms of 
different pole strengths. The terms are obtained 
from the expansion of the potential in a Taylor's 
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series in inverse powers of the interatomic dis- 
tance. If the atoms are different, the dipole-dipole 
and dipole-quadrupole first-order interaction 
energies are zero because of the odd parity of the 
dipole portions of their potential terms. That is 
to say the expression for the energy is the sum of 
pairs of integrals, each referring to one atom 
alone. The integrand of each contains the product 
of two atomic functions, which product is even 
under inversion since the functions belong to the 
same atomic configuration. In addition each 
integrand contains a linear or quadratic function 
of the electron coordinates, referred to the nu- 
cleus of the atom in question. For first-order 
dipole-dipole or dipole-quadrupole interactions 
at least one of the any pair of such integrals 
contains a linear factor in its integrand, which 
being odd under inversion causes the integral to 
vanish. Moreover if either atom has spherical 
symmetry (S state), the first-order quadrupole- 
quadrupole energy is zero. Nonvanishing first- 
order dipole energies are obtained only if the 
atoms are the same and are in states with dif- 
ferent parity. Second-order interaction energies 
do not in general vanish. They are sometimes 
referred to as polarization or induced energies. 
Since for two neutral atoms the expanded poten- 
tial begins with the inverse third power of the 
interatomic distance, second-order energies vary 
with the inverse sixth and higher inverse powers 
of the distance. These interactions are those 
which produce van der Waals forces. They were 
first discussed extensively on the basis of quan- 
tum mechanics by London.' The most careful 
studies have been made for hydrogen and 
helium.” 

Many atoms are not spherically symmetrical 
in the ground state; for these the quadrupole 
moment does not vanish. If we have two such 
atoms, there will be, in general, a nonvanishing 
first-order interaction. The quadrupole-quadru- 
pole potential term varies with the inverse fifth 
power of the distance. London has noted in this 
connection that for two rotating molecules the 


1F, London, Zeits. f. Physik 63, 245 (1930). 

2 Slater and Kirkwood, Phys. Rev. 37, 682 (1931); H. R. 
Hassé, Proc. Camb. Phil. Soc. 27, 66 (1931); H. Margenau, 
Phys. Rev. 38, 747 (1931); Pauling and Beach, Phys. Rev. 
47, 686 (1935). 


quadrupole-quadrupole interaction energy will 
depend on the relative orientation of the two 
angular momenta.* The diagonal terms of the 
energy of interaction of two rotating neutral 
systems have in fact their dependence on orien- 
tation given by 


(1/R*){J(J+1) 


where J, \/; J*, M/* are the angular momentum 
and magnetic quantum numbers describing in 
part the states of the two systems and R is the 
separation distance. As London has pointed out, 
the interaction vanishes on the average over all 
orientations. However nondiagonal terms must 
be considered in an accurate treatment of the 
problem. 


2 


Let us introduce for the two atoms parallel 
sets of Cartesian axes, one with its origin in each 
nucleus, and let each have the z axis along the 
internuclear line. We shall use an asterisk to 
distinguish quantities referring specifically to 
the second atom. Thus coordinates, quantum 
numbers, wave functions (or states) are written, 
for example, x, J, ¥ when they refer to the first 
atom and x*, M*, ¥* when they refer to the 
second. If we have but two electrons, one near 
each nucleus, the expression for the quadrupole- 
quadrupole potential energy is 


(3e?/4R®) Pr? — — — *? 
+2(xx*+yy* —422*)?}. (1) 


It is convenient to write this in terms of the fol- 
lowing functions, 


f(—1) = (3) (x—ty)2, 
f(0) (2) 
=3(3) 


These functions have the property that they are 
partners belonging to the rows of the five- 
dimensional irreducible representation of the 
three-dimensional rotation group, which property 


3F. London, Zeits. f. physik Chemie B11, 222 (1931). 
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is helpful in evaluating matrix elements. Ex- 
pression (1) becomes 


| f(—2)f*(2) +4f(— 1) f*(1) 
+6f(0)f*(0) 1) +f(2) f*(—2)}. 


We define for a given atom the functions 
F(m)=DUf(m), (2a) 


where the summation is over all the electrons 
1---N of the atom. These atomic functions have 
the same rotational properties as the one electron 
functions. Summing (la) over the electrons 
1---N of the first atom and N+1---N* of the 
second atom, we obtain for the quadrupole- 
quadrupole interaction potential energy 


(e?/ R*)a(m) F(m) F*(—m), 
a(—2)=a(2)=1, (3) 
a(—1)=a(1)=4, a(0)=6. 


3 


The calculation of the energy of interaction of 
atoms at large separation distances is com- 
plicated by the fact that the atomic spin-orbit 
energy differences may be of the same order of 
magnitude. The doublet separation is only 
0.0018 ev in the lowest term of boron but in the 
lowest term of thallium it is 0.94 ev. The limiting 
cases in which the energy differences are very 
small or very large compared to the interaction 
energy are most easily treated. For some atoms 
the atomic splitting is so small that it is negligible 
over a rather wide range of R. For other atoms 
it is at large distances much larger than the 
interaction; then at smaller distances the inter- 
action becomes comparable to the splitting; at 
still smaller distances the molecular interaction 
is the dominant portion of the energy. We dis- 
tinguish therefore the cases that the splitting in 
both of the atoms is negligible, the splitting in 
the first atom is large and in the second is 
negligible, and the splitting in both is large. 
Whether one of these cases is applicable to a 
pair of atoms depends of course on the value of 
R. Transition cases involve very complicated 
treatment except in the very simplest cases. 

Let us suppose that Russell-Saunders coupling 
to exist in the states of the atoms. We restrict our 
attention to the states of an atom belonging to a 


definite term, distinguished by a symbol y 
indicating the configuration of which it is a part 
and definite Z and S indicating the total orbital 
and spin angular momenta.* Even or odd parity 
we distinguish by w=+1. If the spin-orbit 
splitting is neglected, all the levels of a term fall 
together and the atomic states are best described 
in the SLMsM,_ coupling scheme. On the other 
hand if the levels of a term are widely separated, 
we need consider only states of a given level, 
indicating the total angular momentum by J; 
under such circumstances the SLJM coupling 
scheme is most suitable. 

Let us now consider the case in which the 
atomic spin-orbit splitting is negligible. The 
characteristic values of the observables of the 
widely separated system are ywSLMsM_; 
y*w*S*L*M* Let ¥(.MsM_) be the Schri- 
dinger representation of the state of the first 
atom and W*(M*;M*,) be that of the second. 
We choose the relative phases so that the func- 
tions are partners belonging to the rows of the 
direct product of the Sth and Lth irreducible 
representations of the three-dimensional rotation- 
reflection group.® In neglecting exchange we 
we could very well use as our functions for the 
system simple products of the atomic functions. 
It is however almost as easy to use functions 
antisymmetrical in all the electrons. Such a 
presentation has the advantage that the cor- 
relation between atomic and molecular states is 
more apparent. 

We first note that the original interaction 
energy was symmetrical in all N+ N* electrons 
of the system while (3) is not. An unsymmetrical 
expression has been obtained from a symmetrical 
interaction by assuming particular electrons to 
belong to a given atom; a specific assignment is 
necessary if we are to expand the potential at 
this point of the treatment. If we use completely 
antisymmetrical wave functions, the expansion 
cannot be made until the functions have been 
substituted in the matrix integral and the latter 
reduced to expressions in which there are specific 
assignments of the electrons in each function. 
Those terms in which electrons are assigned to 


*We use for the most part the atomic terminology of 
Condon and Shortley, Theory of Atomic Spectra. 

5 More accurately of the “unitary group” or the group 
of the two-dimensional unitary matrices of determinant 1; 
see. for instance, Wigner, Gruppentheorie. 
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different atoms will be small under the assump- 
tion of small overlapping. We therefore retain 
only terms in which the assignment is the same 
for each atom. It is then possible at this final 
stage to substitute expressions like (3) for the 
interaction. 

It is simpler, however, to define from the 
start a function similar to (3) but having the 
desired property of symmetry and giving the 
same result for the energy as obtained by the 
process outlined above. We add to the definition 
of W the restriction that it vanish when the 
electrons 1---N are not associated with the first 
atom and then sum over all exchanges Q, that 
is all permutations which lead to an ordered 
assignment in each atom. We obtain the sym- 
metrical interaction 


(e?/ R®)T=(e/R®) SOW. (3a) 


As before, in computing energy matrices we 
finally obtain elements in which the electrons 
are assigned to specific atoms; we retain only 
those parts in which the assignment is the same 
for both functions for each atom. On introducing 
T there is only one term in the sum over Q 
which does not give zero, namely that which has 
the same assignment. 

The symmetries of the system are those of 
rotations around the internuclear line, reflections 
in planes through this line, and in the case of 
like nuclei of inversion, that is reflections in the 
midpoint between the nuclei holding the nuclei 
fixed. The states of the system will, correspond- 
ingly, have quantum numbers A and & for the 
total orbital angular momentum and spin along 
the axis, be positive or negative (of importance 
only in » states), and for like nuclei be g or u. 
The possible molecular states and their relation 
to the states of the atoms when widely separated 
have been considered in detail by Wigner and 
Witmer.® 

Let S be the total spin of the system. 
Its possible values are S+.S*, S+S*—1,---, 
| S—S*|. We define the functions 


MsM*,)(SS*Mg2— Ms|SS*S>), (4) 


® Wigner and Witmer, Zeits. f. Physik 51, 859 1928; also 
R. S. Mulliken, Phys. Rev. 36, 1440 (1930) ; F. Hund, 
Zeits. f. Physik ‘63, 723 (1930). For a good discussion see 
H. Sponer, Molekiilspektrum, II pp. 132-137. 


37 


which have been built up out of simple products 
by use of the vector model and belong to the Sth 
irreducible representation of the three-dimen- 
sional rotation group. The coefficients on the 
right are the elements of the usual transformation 
matrices for angular momenta. 

Take first unlike atoms and A+0. Let 
The two antisymmetrized func- 
tions 

A—M1), 
(5) 
c(—)4¥(Sz, —A, — Mz) 


are partners belonging to the two rows of the 
Ath representation of the two-dimensional rota- 
tion-reflection group. The operator A’ is defined 
as? 


A'=(NIN*!/(N+N*) 


where Q, as before, is a permutation which leads 
to an ordered assignment of electrons in each 
atom and g is the parity of Q. For a given A, 
M, takes on the following values (L 2 L*) 


L+L*2A2L—L*+1, 
M,=L,L—-1, ---,A—L*; 
L-L*2A>0, (6) 


M,=A+L*, A+L*—1, ---, A—L*. 


This can be seen either by constructing a diagram 
or by calculating the character of the reducible 
representation. In addition we have for A=0 


a 1 
¥ o(Sz, M,)}, 
M,=L*, L*—1,-+--->0, (Sa) 
¥(S> b0) =A’ 00). 


The functions are positive (2+) if c=—1, 
negative (L~-) if c=1. On the other hand the 
functions ¥(b) are positive if c=1, negative if 
c=-—1. For a given c these functions belong to 
different irreducible representations. The matrix 
T is diagonal in S, , and A (also a and b) and 
moreover elements in which both states are 
replaced by partners will have the same values. 
The diagonality in S results from the fact that 
T is independent of the spin and therefore is 


7 Condon and Shortley, reference 4, p. 215. 
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symmetrical under all rotations of the spin axes 


alone. 
Let us compute the matrix of W using the 


functions (4). It is no longer necessary to carry 
the quantum numbers SY. We find, because of the 
unitary nature of the transformation coefficients 
and the fact that the spin is not contained in W, 


(M1! 
= Ya(m)(LM,'| F(m)|LM1) 
(L*M* F*(—m)!L*M*,). (7) 


Now the functions F belong to the representation 
2 of the three-dimensional rotation group. The 
product F(m)¥(.W/,) occurring in the first of 
the elements on the right can be resolved by the 
vector model into a sum of functions x belonging 
to representations L’ according to the relation 


F(m) == x(L2L'Mi+m) 
(8) 


where again we have as coefficients the trans- 
formation functions for angular momenta. The 
summation is over the values L’=L+2, --., 
|L—2_). Except for the case L =0 (S state), there 
is a x in the sum belonging to the Lth repre- 
sentation. We have 


(LM,'| F(m)| LM) =A(L)6(Mz', M,+m) 
X(L2M (9) 


using the orthogonality theorems for the scalar 
product. A(L) is the scalar product of V(.1/,) 
and x(.V/_) and is therefore independent of 
and m. In this way we obtain for (7) 


(M ,!M*,!|W| M,M*,) =A(L)A*(L*) 
x M,+M*,)a(M_' M,) 
(L2M,M_' M, | L2LM,') 

(L*2M*,M*,'’— M*,| L*2L*M*;'). (7a) 


The non-zero elements of T for the states of 
(5) and (5a) (unlike atoms) are found to be 


Mz, A—M1), 


a a 


(6M — W 00), 


(b0| T'| 0) = W| 00). 


Take next like atoms. The additional sym- 
metry gives a reduction of the energy deter- 
minant, for g and u states do not mix. We define 
the functions 


(ss sats) A—M_z) 
M1)}, (11) 
M,.=L, L—-1, --->}A, 
and in addition if A is even 
=A’ (11a) 


¥(e) is g if S is odd (molecular triplets, etc.), 
u if S is even (molecular singlets, quintets, etc.). 


e 
( AM,'|T| A—M,'|W! M1, A—M,'|W|A—Mz, M1), 


On the other hand y(/) is g if S is even, u if § 
is odd. Moreover the states ¥(e0.1/;) are negative 
and the states ¥(f0./,) are positive. We have 
then the correlation (A=0) that negative states 
are gif Sis odd and u if S is even, while positive 
states are g if S is even and u if S is odd. These 
functions belong to different representations of 
the group of inversions in the midpoint between 
the atoms. When A#0, they belong to the 
second row of the Ath representation of the two- 
dimensional rotation-reflection group, their part- 
ners being the functions obtained by changing 
the signs of A and M, and multiplying by (—)*. 
The non-zero elements of 7 for the states of (11) 
and (11a) are 


12) 


(f{4M ;!|T|fA3A) A—M,'| 


=(3A3A| 


| 

| 

| 

| 

| 

| 


(10) 
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Let us next consider briefly the case of two 
different atoms in one of which the atomic spin- 
orbit splitting is large and in the other it is 
negligible. We describe the states of the first 
atom in the SLJ.\/ coupling scheme, keeping the 
second in the SLs, scheme, and consider 
only molecular states arising from a given 
J level. Introducing the new function 


¢o(JS*M* (13) 


we have at large distances the molecular states 
represented by the functions 


V(JS*M* (14) 


for ' #0. If J is an integer, [ can in addition 
have the value zero, for which value we define 
the functions 


a 1 


o(JS*M*s, —M,M)}, (14a) 
V(JS*M* b 0) =A’ e(JS*M* 00). 


The possible values of 7 in these expressions 


depend on whether L* is larger or equal to or 
less than J. We shall not write down the detailed 
expressions, for they are quite similar to the 
corresponding restrictions of the last section in 
any given case. There is one peculiarity which 
we note. The simultaneous rotation of spin and 
position axes does not give us a transformation 
group that is particularly useful. Instead we 
might introduce an artificial type of spin trans- 
formation depending on the space state of the 
particle in question. This is not necessary how- 
ever, for one readily sees that T is diagonal in 
and We define the expression 


(J.M'M* JMM*,_) 
= Ya(m)(JM’| F(m) | JM) 
—m) | L*M*,) 
= C(J)A*(L*)6(.M'+ M*,', M+M*,_) 
Xa(M’— M)(J2MM'— M 
(L*2M*,M*,' — M*,| L*2L*M*,'). (15) 
C(J) is a scalar product independent of AV and is 


zero if J is zero or one-half. In terms of this 
expression the non-zero elements of 7 are 


(JTM'|T| ITM) =(JM', JM, 


a a 
M’| u)= —M'|W\ JM, —M)#(JM’', —M'|\W|J, —M, M), 


(JoM’' T | Jb0)=v2(J.M', —M'| W) J00), 


(Jb0| T| Jb0) = (J00) W| 


5 


We come finally to the case in which the spin- 
orbit splitting of the atomic term energies is 
large in both atoms. All atomic states are repre- 
sented in the SLJ.M scheme. The molecule has 
a quantum number © for the total angular 
momentum along the internuclear line. The 
components of the orbital angular momentum 
and spin are not defined separately. We introduce 
the function 


¢o(JI*MM*) =¥(JM)v*(J*M*). (17) 
Consider unlike atoms and 20. Let 


d=ww*(—)/+", 


(16) 
The two antisymmetrized functions 
¥(JI*QM) =A'o(JJ*M, Q-—M), 
d(—)*¥(JJ*, —2, —M), (18) 


are partners belonging to the two rows of the Qth 
representation of the two-dimensional rotation- 
reflection group in which the rotation of the spin 
as well as the position axes is carried out. For a 
given 2, takes on the values (J 2 J*) 


J+ J*> 

J—J*> Q>0, (19) 

M=04J*, 24+J*-1, ---, Q-—J*. 


If J and J* are integers or are half-integers, 2 


trans- 
. The 
there 
repre- 
, (9) 
scalar 
(M1) 
(7a) 
tes of 
™ 
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can in addition have the value zero. We define 
in case 2=0 the two sets of functions 


a 1 
M)=—A'le(J PM, —M) 
o(JJ*, —M, M)}, 
M=J*, J*—1, --->0, 


(18a) 


including if J and J* are integers, 
¥(JJ*b0) = A’y(JJ*00). (18b) 


The functions ¥(a) are positive (0+) if d= —1, 
negative (0-) if d=1 while the functions y(d) 
are positive if d=1, negative if d=—1. The 


matrix T is diagonal in 2 (and a and 6). We can 
express the results in terms of the expression 


(JJ*M'M* | W| JJ*MM*) 
= Ya(m)(JM’| F(m)| JM) 


x (J*M* | F*(—m) | J*M*) 
M+ M*) 
x (J*2.M*M* — M*| J*27*M*). (20) 


The non-zero elements of T for unlike atoms are 
found to be 


(JJ*QM’ | T| JJ*QM) =(JJ*M', Q2-M), 


a a 
M’'|T|JJ* =(JJ*M’, —M’| W|JJ*M, —M)¥(JJ*M’, —M’'|W|JJ*, —M, M), 


(JJ*bM’ | T| JJ*b0) = v2(JJ*M’', —M'| W| JJ*00), 


(JJ*b0| T | J*b0) = (JI*00| W| JJ*00). 


For like atoms in states of the same configura- 
tion and term (y*=y, S*=S, L*=L, w*=w) and 
the same J values, we define the functions 


1 
Fo(JJQ—M, M)}, (22) 
M=J, J—1, >}2. 


If J is an integer and 2 is even or if J is an half- 
integer and @ is odd, we have in addition the 
function 


e e 
(42 QM’ | r\sy oat) =(JJM’', 2—M’'|W| JIM, 29—M)#(JIM’, M), 


(21) 


V(JIF23Q) =A’ o(JIF (22a) 


¥(e) is wif J is an integer, g if J is a half-integer. 
¥(f) is gif J is an integer, u if J is a half-integer. 
When 20, these functions belong to the second 
row of the Qth representation of the two- 
dimensional rotation-reflection group, their part- 
ners being the functions obtained by changing 
the signs of 2 and M and multiplying by 
(—)?/+®. For 2=0, g functions are positive and u 
functions are negative. The non-zero elements of 
T for like atoms with J*=J are found to be 


(23) 


(JIFQM’ | T | Q— M’| W| Q32), 


(JIFQZQ| T | JTFQZQ) = (FIZ W! 


The remaining possibility of like atoms in 
states belonging to the same configuration and 
term but with different J values leads to ex- 
pressions which are somewhat more complicated 
than those considered because the coefficients 
involve both J and J*. Since we are primarily 


interested in the lowest states of the atoms, we 
shall not discuss this case. 
6 


The elements of the matrix T in the cases 
discussed have as a common factor the product of 


| 
| 
T 


can 
on 


(20) 


5 are 


(21) 
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TABLE I. A(hh*5) for integral 5. If 5 is not zero there are 
two states for each number listed. If the atomic term splitting 
is negligible in both atoms (h=L, h*=L*, 5=A), the letters 
mean the following: Let c=ww*(—)"*4*. If the atoms are 
different, an a after a number means that the state is positive 
(=+) if c= —J1, negative (=~) if c=1, while b means that the 
state ts positive if c= 1, negative if c= —1. If the atoms are the 
same, an e means the state is g 1f the molecular spin ts odd, 
u if tt is even. An f means that the state is g if the spin is 
even, u tf it ts odd. Also for \=O states with e are negative 
and states with f are positive. For other cases see the text. 


h=1 2 3 4 Ls 2 
3 | h*=1 1 1 1 Ls 2 
ca | —1.014a| -1.225a| -1.305a| 2.77400 
0 ea | —1.060¢a 
3.6fb | 2.3246] 2.0966] 2.0136] 3.2/6 | 3.261fb 
0 fb | -0.2956| —0.3826| —0.4186| | 1.50876 
~0.484 fb 
Oe 2.456 | 2.061 1.972 | 1.202 ¢ 
0.047 | —0.170 | —0.289 -~0.917¢ 
—1.489 | -1.401 | —1.392 
-2.4f 1.127 f 
-o8f | -1.4137 
2 0.356 | 2.199 1.367 O4e | —2.000¢ 
0.230 | 0.347 
—1.939 | -1.424 
0.6/ -28f | oss7f 
-1857f 
3 1.014 | 0646 | 2.032 0.857 ¢ 
—3.096 | —0.032 
—1.710 
1.2 | -2.571¢ 
4 1.225 | 0.853 
—3.222 
5 1.354 


two atomic coefficients. Dividing by this common 
factor we obtain the elements of a dimensionless 
determinant the roots of which are needed to 
describe the quadrupole-quadrupole interaction 
between atoms. 

In general let h, h* be a pair of numbers, 
integers or half-integers, corresponding to L, L*; 
J, L*;or J, J*. Let 6 stand for A, T, or Q as the 
case may be. The range of the remaining 
parameter or is essentially the same in 
all three cases. A comparison of (7a), (15) and 
(20) shows that for the same numerical values 
h, h* the elements of W are the same aside from 
the atomic coefficients. If h* is numerically 
different from h, we see from (10), (16) and (21) 
that on dividing out the atomic coefficients the 
dimensionless determinants are the same. For 
unlike atoms with h*=h there is one further 
simple reduction of the determinant possible. In 
this case (L*=L, L*=J, or J*=J) we can form 
linear combinations ¥(e) and ¥(f) which are 


formally the same as for like atoms. It can be 
shown that 7 has no cross terms between these 
two groups and that the determinant is the same 
as obtained from the expressions for like atoms, 
that is from (12) or (23) on dividing out the 
coefficients. Thus we have the simplifying result 
that the dimensionless determinant depends only 
on the numerical values of the pair of numbers 
h, h*. We shall designate the roots of a sub- 
determinant by A(#h*5). How many roots there 
are for a given 6 depends in general on how many 
times the 6th representation is contained in the 
reducible representation obtained the 
transformation under two-dimensional rotations 
and reflections of all the possible simple product 
functions of the atomic states. 

Roots of the dimensionless determinant for 
some values of h and h* are given in Tables I 
and II. 
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The atomic coefficients A(L) and C(J) are 
defined by the relations 


A(L)(L2M ,0|L2LM;) 
(24) 
C(J)(J2M0| J2JM) =(JM | F(0)| JM). 


To obtain explicit expressions, we transform to 
the zero-order representation scheme charac- 
terized by the one-electron quantum numbers 
a=nim,m,. If we let A stand for the collection of 


TABLE II. A(hh*5) for half-integral 5. 


=1.5 2 4.5 
h*¥ =1 1.5 1 
0.5 2.746 2.864 1.977 
0 1.551 — 0.399 
— 1.049 —0.377 — 1.345 
— 1.168 

1.5 0.176 0.916 1.916 
— 2.722 — 1.195 — 0.206 
— 1.872 — 1.476 
2.5 0.848 0.619 1.883 
—2.770 0.125 
—1.774 
3.5 1.434 2.054 
0.518 
— 2.339 
4.5 0.934 
— 3.270 
5.5 1.401 


| | 
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such quantum numbers in a zero-order state of 
the atom, we have for A(L) the relation 


A(L)(Z2M,0| L2LM1)=> (ySLMsM_/A) 
A.B 


X(A| (25) 


But only A =B terms occur or are different from 
zero. For (A| FB) vanishes if A differs from B 
by more than one individual electron quantum 
number set and has the value +(a fd) if all sets 
but a#b are the same. But even this term is 
ruled out since if all but one of the m, are the 
same, they must all be the same because of the 
condition 2m,= M, for both A and B. We have 
therefore for (25) the expression 


(ySLMsM1| (a f\a)}. (26) 
But 
1(/+1) 
X (12m | 121m (27) 
where 
= “drPR(nl), (28) 


R(nl) /r is the radial part of the wave function for 
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an electron in an nl state. We note that since the 
coefficient 
3m ?—1(1l+1) 


(0 | 121m ,) = (29) 


when summed over the electrons in complete 
shells gives zero, contributions to A(L) will come 
only from incomplete shells. All atoms whose 
normal states are not S states have one incom- 
plete shell having angular momentum. We can 
write 


A(L) = 


(30) 


The second factor in each case is the average 
value of 7* for electrons in that shell. We have 


1(1+1) 
(27 —1)(2/+3) 


Y 
L2LM1) 


BL) =(-)| 


Likewise if we start in the SLJ.W scheme we 
obtain 


1(l+1) ] 
(2/—1)(2/+3) 


TABLE III. Atoms with incomplete p-shells. A star after an element means the spin-orbit splitting of the lowest energy term 
is less than 0.036 ev. Abbreviations (s.c) and (s.c.f) are used to distinguish numbers obtained using shielding constant and self- 


consistent field wave functions, respectively. 


GrounD LEVEL rm IN BOHR RADII (r?),y(np) IN (BOHR RADII)? 

ELEMENT np B(L) D(J) (s.c) (s.c.f) (s.c) (s.c.f) 

B* 2p — 0.6325 0 1.54 1.80 4.4 9.12 

Al* 3p 2.57 10 

Ga 4p 2.74 1 1 

In 5p 3.20 14 

Tl 6p 3.53 17 

c* 2p? 3Py +0.6325 0 1.23 1.26 2.8 4.882 

Ge 4p? ‘ 2.42 8.5 

Sn 5p? ae ae ‘ 2.83 1 1 

Pb 6p? 3.12” 14 

ap, — 0.6325 — 0.3743 0.88 0.85 1.4 2.440 

Se 4p 1.97 5.6 

Te 5p! “ “ 2.30 75 

F 2p spe, +0.6325 +0.4472 0.77 0.73 1.1 1.82 

Cl 3p “ “ 1.48 3.4 

Br 4p > 1.80 4.7 

I > 1 1 6.2 


the 


31) 
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The transformation coefficients to the zero-order 
scheme are usually simplest for extreme values 
of MsM, or M. They can often be obtained for 
extreme values by inspection. Otherwise they are 
readily obtainable with the help of the angular 
momentum operators as used by Gray and 
Wills.* 


8 


We come now to the problem of obtaining the 
numerical values of (7*)y. In Tables III and IV 
we have listed most of the elements for which 
these considerations are of interest. The second 
column in the tables gives the nature of the 
lowest level of the atoms. We restrict our 
discussion to the lowest term or level. 

The average value of r? has been estimated for 
each atomic valence / electron by using the 
shielding constants and effective principal quan- 
tum numbers of Slater.* Such a method gives 
unreliable values, however, especially for light 
atoms. This is because the major contribution to 
the average value comes in the expression (28) 
from values of r out beyond the maximum of the 
density function. In this range it is well known 
that simple functions with shielding constants 
are poor approximations. The error is presumably 
less for an electron with a large principal quantum 
number, 7, since the peak of the density function 
is sharper and the value obtained more nearly 
the square of the radius at the maximum, r,,. The 
magnitude of the error for small m is shown by 
the values which we have calculated for B, C, O 
and F using wave functions obtained from the 
self-consistent field.” The numbers so obtained 
are approximately twice those calculated from 
shielding constants. 

The self-consistent field probably gives wave 
functions which are too large at large distances. 
Electrons are kept away from the neighborhood 
of the nucleus because of the static nature of 
the electron field which the method uses. This 
tendency to favor a more diffuse distribution is 


§ Gray and Wills, Phys. Rev. 38, 248 (1931). 

*J. C. Slater, Phys. Rev. 36, 57 (1930). 

1 Self-consistent field data have been calculated from 
functions given by the following: B, Brown, Bartlett and 
Dunn, Phys. Rev. 44, 296 (1933); C, C. C. Torrance, 
Phys. Rev. 46, 388 (1934); O, Hartree and Black, Proc. 
Roy. Soc. 139, 311 (1933); F, F~, and Ne, with exchange 
hie p electrons, F. W. Brown, Phys. Rev. 44, 214 


TABLE IV. Atoms with incomplete d shells. Numbers in 
parentheses are positions of maximum radial density of the 
two electrons in the s shell lying outside the incomplete d shell, 
which electrons are most effective in determining the radius 
of the atoms. 


Grounp LEVEL rm IN (r2)ay(nd) 
Bour (Bour 
ELe- Rabu Ravi)? 
MENT nd B(L) DJ) (s.c) (s.c) 
Sc* 3d4s? 2Daye | —0.5345 | —0.4472 | 3.00(4.55) 14 
La 5d6s? 5§.33(5.88) 40 
Ti 3d245? 3F —0.2213 | —0.1833 | 2.47(4.35) 9.5 
Zr 4d25s2 ” 3.75(5.08) 20 
Hf 5d*6s? 4.38(5.60) 27 
3d3452 | +0.2213 | +0.1533 | 2.09(4.15) 68 
WwW 5d45s2 +0.5345 0 3.23(5.11) 15 
Fe 3d6452 —0.5345 | —0.4006| 1.44(3.65) 3.2 
Co 3d74s2 {Foye | —0.2213 | —0.1935 | 1.32(3.51) 2.7 
Ni 3d*45? +0.2213 | +-0.2003) 1.21(3.38) 2.3 
Ir ?2Dsj2 | +0.5345 | +0.4782 | 2.58 94 


more clearly the case when no account is taken of 
exchange. The value of (7*)y is a very sensitive 
function of such effects. For the chlorine negative 
ion the value of (7*)4(3p) has been calculated with 
exchange to be!! 5.137a,? (a9 = Bohr radius). The 
closely related diamagnetic susceptibility is 
calculated to be —30.4-10~*. Without exchange 
the diamagnetic susceptibility is calculated to be 
—41.3-10-® which corresponds to a value of 
of about 7.449”. The self-consistent field 
functions for oxygen and carbon were obtained 
without exchange and hence the values for 
(r*) obtained for these atoms are probably 
much too large, judging from the effect of 
exchange on Cl-. The importance of exchange in 
neutral atoms may be slightly less however. The 
function for F (also Ne and F~, Table V) was 
calculated with exchange between the 2) elec- 
trons. Boron has but one 2 electron. The values 
for these atoms are probably quite good.” 
Additional information is obtainable from the 
diamagnetic susceptibility of closely related 
atoms. The diamagnetic susceptibility of a gas of 
atoms in singlet states is given by Van Vleck’s 
expression, —0.79-10-®2(r’),, where the sum is 
over all the electrons of an atom. In Table V we 
have some experimental values for Ne and A 
obtained by Havens and in addition some values 


" Hartree and Hartree, Proc. Roy. Soc. 156, 60 (1936); 
Cl~ with exchange. The value of 7.4 for (r? )y(3p) without 
exchange is only approximate since the numbers for inner 
electrons without exchange are not given. 

”® Unfortunately some doubt as to the accuracy of the 
functions is cast by Brown's result for the electron affinity 
of fluorine. He obtained a negative number for the electron 
affinity when effects due to the change in the function on 
removal of the electron were included. 
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for the negative halogen ions obtained by 
Brindley and Hoare from salt crystals using the 
assumption of additivity of the susceptibilities of 
the two ions of the salts.'"* The magnetic sus- 
ceptibility of Ne obtained from the self-con- 
sistent field with exchange between the 2p 
electrons is only five percent greater than the 
experimental value. If we calculate back from the 
experimental values using Slater values for the 
inner electrons, which contribute only a small 
fraction of the sum, we obtain the numbers for 
(r*) listed in the fifth column (the numbers in 
parentheses have been obtained from self-con- 
sistent field values for the inner electrons). We 
see that for Ne the difference is about ten 
percent. However with ions the agreement is not 
so good. The experimental values are much 
smaller than the values from the self-consistent 
field. This result is not surprising since each ion 
is in the disturbing field of its neighbors. The 
crystals are not completely ionic; Ewing and 
Seitz have shown by the application of the 
Hartree approximation to LiH and LiF crystals 
that there is considerable charge in the neighbor- 
hood of the metallic ion." The picture of the 
valence electrons localized in a region about the 
electron-negative atoms is not a very accurate 
one. 

We have a lower bound for (r?),, for chlorine in 
the number 2.88 for argon and an upper bound in 
the number 5.137 for the chlorine negative ion. 


9 


If one of the three cases discussed in Sections 
3, 4, and 5 is applicable, the quadrupole- 
quadrupole interaction energy between two 


TABLE V. Data for some closely related atoms. 


r?)ay(np) IN (BOHR 
MAGNETIC SUSCEPTIBILITY 108 ( Rapit)? 
FROM 
Expt. (s.c) (s.c.f) Expt. (s.c) | (s.c.f) 
Ne — 7.651+0.008 | — 5.58} — 8.08} 1.20(1.23) | 0.88 | 1.32 
A —19.23 +0.20 —18.5 2.88 2.77 
F- — 94 — 8.13] —17.1 | 1.48(1.50) | 1.28 | 3.12 
cl- —24.2 —25.2 | —30.4 | 3.64(3.82) | 3.81 | 5.137 
Br- —34.5 —39.3 4.38 5.14 
—50.6 —58.6 5.58 6.84 


3G. G. Havens, Phys. Rev. 43, 992 (1933). Brindley and 
Hoare, Proc. Roy. Soc. 152, 342 (1935). 
4 Ewing and Seitz, Phys. Rev. 50, 760 (1936). 
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atoms is given by one of the expressions 


(e?/ R®)A(L)A*(L*)M(LL*A), 
(e?/ 


The results would be of more significance if 
we knew at what distances other forces begin 
to be of equal strength. It is not easy to say 
however. At large distances the Heitler-London 
method gives the repulsive forces to a good 
approximation. From Sugiura’s formula'® we 
obtain for hydrogen at a separation distance of 
4a, for the lowest singlet and triplet states 
—0.27 ev and 0.20 ev, respectively. The induced 
polarization or van der Waals attraction is 
slightly smaller at this distance. Using the 
coefficients of Pauling and Beach,’ we obtain for 
the first three terms (terms in the inverse sixth, 
eighth and tenth powers of R) the numbers 
of the expression — (0.0429+0.0513+0.0293) ev 
= —0.124 ev.'® For purposes of comparison let us 
find the distance at which the Heitler-London 
repulsion between two helium atoms is 0.2 ev. 
Slater has calculated this energy'’ and has 
obtained the expression 484 exp (—2.43R/ao) ev. 
For 0.2 ev we find R=3.2a). At R=4r,, 42.2ao 
the repulsive energy is 2.3 ev (7,, is the position of 
the maximum of the radial density function). 
Using Margenau’s coefficients we find for the 
van der Waals energy of two helium atoms 
— (0.039+0.030+0.011) ev =—0.080 ev at 
R=3.2a). The Heitler-London energy, being 
essentially a penetration energy, builds up and 
falls off faster than an inverse power energy. As 
pointed out in the preceding section, for heavy 
atoms the electron distribution function dies off 
more rapidly beyond the last maximum than for 
light atoms. If for comparison we use a distance 
of twice the sum of the atomic radii, that is twice 
the sum of rm, the exchange forces should be 
relatively smaller at such distances for heavy 
atoms than for light. For atoms with incomplete 
d shells the situation is complicated by the 
presence of an outside s shell. 


(33) 


1 Y, Sugiura, Zeits. f. Physik 45, 484 (1927); see H. 
Bethe, Handbuch der Physik, Vol. 24, p. 537. . 

16 In hydrogen the van der Waals energies are more im- 
portant than exchange energies for R greater than about 
7ao. For R=7a the three van der Waals terms are in the 
proportion 1, 1/2, 1/8. The approximations made in cal- 
culating the coefficients are poor for smaller distances. 

17]. C. Slater, Phys. Rev. 32, 349 (1928). 
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TaBLeE VI. Energies of attraction at distances of twice the 
sum of the atomic radii. The atomic radius is defined as the 
position of the last maximum in the electron density functicn. 
The spin-orbit coupling has been neglected and the values for 
(r? wy taken from the last column of Table III. The last four 
molecules have in addition attractive A states at this distance 
with energies one-sixth the energies of the = states and the 
same multiplicities (neglecting other forces). The known states 
listed are those which have been observed and which on dis- 
sociation go into normal states of the atoms. 2r, is for the 
equilibrium position of the lowest observed state. 


KNOWN STATES |27¢ 108 cm) R108cm|_ E’ev STATES 

Be 3.80 —0.112 

BO} 2.4 2.80 —0.139|} 

Ce | 2.62 2.66 —0.190) 
CF 2.10 —0.231 “11, 

Oz | 'Zg*, u* 24 1.80 —0.341 My, ly 
| 154 | -0.406) 

BC 3.23 —0.202| 

BF 2.67 —0.195 

OF 1.67 —0.549; 


If we have two atoms in P states and neglect 
the atomic multiplet splitting, we see from the 
first column of Table I that the states of the 
system fall into four energy levels due to the 
quadrupole-quadrupole interaction of the two 
atoms. In Table VI are given the energies for the 
lowest (attractive) levels at distances of twice 
the sum of the atomic radii for the molecular 
pairs obtainable from B, C, O and F. Consider Be 
and F,. The atoms have spin 3, which give 
S=0, 1. Also c=1. From Table I we see that the 
lowest states (A= —2.4) are f states. This means 
that 'II and ‘II are g and u, respectively. Similarly 
we find that for the next level, which has zero 
energy, the states are 'Z,~, 
s11,. The third level has the states 'A,, *A, and an 
energy —} that given in Table VI, that is it is 
repulsive. The fourth level has the states 'D,*, 
8y,,+ and an energy —1.5 that of the lowest level. 

For Cz and Oz since the spins are both 1, 
S=0, 1, 2. The lowest level has the states 'II,, 


‘II,, *Il, for which the energy is given in Table 
VI. The states 'Z.~, *Z,-, *Z.~, 'Z,*, *Z*, *Z,*, 
'II,, *II,, *Il, have zero energy. The states 'A,, 
are repulsive, the energy being — } that 
of the lowest level, and the states 'X,*, 
are also repulsive, having energy —1.5 that of 
the lowest level. For BO and CF since the 
atomic spins are }, 1, S=}, 3. Also c= —1. AA* 
is positive. The lowest level has the states *II, II. 
The states *=+, ‘=+, ‘Il have zero 
energy ; the states *A, ‘A have energy — } that of 
the lowest level; the states ?2~, ‘Y~ have energy 
—1.5 that of the lowest level. 

For BC and OF the atomic spins are }, 1 and 
therefore S=}, 3. We see that c= —1 and that 
AA* is negative. The lowest level is therefore 
given by \=3.6. Its states are ?2~, *L~. The next 
energy is also attractive, } that of the lowest 
level, and has the states *A, ‘A. The states 
211, ‘1 have zero energy; the 
states *II, ‘II are repulsive and have an energy 
— 3 that of the lowest energy. For CO the spins 
are 1 and S=0, 1, 2. C=1 and AA’* is negative. 
The lowest level has the states The 
next has {th the energy and the states 'A, *A, °A. 
have zero energy. The states 'II, *II, ‘II are 
repulsive with an energy —% that of the lowest 
level. The level scheme of BF is the same as CO 
except that both spins are } and therefore there 
are no molecular quintets. 

The writer takes pleasure in acknowledging his 
indebtedness to Professor Wigner, who brought 
the problem to his attention, for advice and 
stimulating discussion. He wishes to thank the 
Wisconsin Alumni Research Foundation for its 
support of this work. 
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The Process of Diffusion in a Centrifugal Field of Force 
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The general equation for the settling of particles and molecules in a liquid suspension or ideal 
solution in a centrifuge has been solved for the case of a sector-shaped cell extending tothe center 
of the centrifuge. Curves are given which show the concentration at all points in the cell for 


different values of the time. 


F a liquid solution be placed in a gravitational 
field of force, the distribution of concentration 
changes from a uniform distribution initially to 
an exponential one after an infinite time. The 
differential equation which represents the con- 
centration at a point in the solution at any time 
t has been completely solved by Mason and 
Weaver! for the case of a uniform gravitational 
field. It is the purpose of the present paper to 
consider the differential equation which applies 
to an ideal dilute solution placed in a centrifugal 
field of force. This equation has been derived by 
Lamm? and an approximate solution of it has 
been obtained by Faxen.* But the solution given 
by Faxen cannot be used conveniently for nu- 
merical computation and consequently has been 
of little use for comparing theory and experi- 
ment. 

Consider a dilute solution in a sector-shaped 
cell in a centrifuge which is rotating with a 
constant angular velocity w. In Fig. 1 OAB isa 
cross section of the cell in which the solution is 
confined. The distances of two cylindrical sur- 
faces from the axis of rotation which passes 
through O, are r and r+dr. Let the concentration 
at the surface r be c, where c is understood to 
mean the number of dissolved molecules per cm* 
of solution. The quantity of dissolved substance 
crossing the surface at r in an outward direction 
in unit time is 

Q=(w*rsc — D(dcdr))r6h, 


where h is the depth of the cell. In this ex- 
pression s, which is the velocity of settling of 
the dissolved substance in unit gravitational 
field, is known as the sedimentation constant. 


! Mason and Weaver, Phys. Rev. 23, 412 (1924). 


? Lamm, Arkiv Mat., Astro. Fysik 21B, No. 2 (1929). 
3 Faxen, Arkiv. Mat., Astro. Fysik 21B, No. 2 (1929). 


D is the coefficient of diffusion which will be 
defined precisely below. The net accumulation 
of substance per second in the layer of thickness 
dr is —(0Q dr)dr, which is also represented by 
(dc dt)réhdr. Thus the following relation holds 


0 
( wrsc — = 


r rar 


The solution of (1) will give the concentration c 
as a function of r and ¢. 

The diffusion constant D can be expressed in 
terms of other quantities as follows. In the final 
steady state the quantity of dissolved substance 
moving outward under the action of the cen- 
trifugal field is equal to the quantity diffusing 
inward under the action of the concentration 
gradient. Consider a single particle of the sub- 
stance in solution. If a gravitational force X acts 
on it, it will acquire a terminal drift velocity »v, 
since it is moving through a viscous medium, 
which is given by Stoke’s law as v= Xf, where f 
is the frictional constant. In consequence of this 
velocity there will flow per second through an 
area r0h a quantity of dissolved substance 
ver6h or 


(2) 


For the steady state it can be shown‘ that the 
concentration at any point 7 is given by 


c=c)' exp [M(1— Vp)#*r’/2RT ], (3) 


where Co’ is the concentration at r=0, MV is the 
molecular weight of the dissolved substance, 
V its partial specific volume, p the density of 
the pure solvent, R the gas constant and T the 


4 Svedberg, Colloidal Chemistry. 
746 


[IE 53 


ill be 
lation 
kness 
by 
Ids 


(1) 
ion 


ed in 
final 
tance 
cen- 
using 
ation 
sub- 
acts 
ity 2, 
lium, 
vere f 
f this 
h an 
ance 


(2) 


t the 


(3) 


s the 
ince, 
y of 
the 


DIFFUSION IN A CENTRIFUGAL FIELD 747 


absolute temperature. From (3) one gets 


dc Nm(1—Vp)e*r 
— (4) 
or RT 


where N is Avogadro’s number and m the mass 
of a single molecule. But the quantity m(1— Ip) 
may be regarded as the effective mass of the 
molecule (i.e., the actual mass corrected for 
the effect of buoyancy) so that 


m(1— (5) 
Therefore, from (4) and (5) 
dc dr=NXc RT. (6) 


The number of molecules which, as a result of 
the field, cross in one second a surface at r where 
the concentration is ¢ is, from (2) and (6) 


RT dc 
—- —roh. 
Nf or 


This is also equal to the number Dr6éh(dc dr) of 
molecules crossing the same surface in the 
opposite direction as a result of the concentration 
gradient. Thus we have 


D=RT/Nf. (7) 


But m(1— =fv from Stoke's law and 
v=w’rs from the definition of s. Therefore 


sf=m(1— Vp). 
Substitution of this value of f in (7) gives 
D=RTs/M(1—Vp). (8) 


We return to Eq. (1) and try a solution of 
the form 
c(r, 


where R is a function of r alone and T is a func- 
tion of ¢ alone. On substitution this gives 


1 
(0 r| = 
rR dr T dt 


In this equation the variables are separated and 
the two equations 


1/T(dT ‘dt)=—8 (9) 


are obtained, where 8 is a constant. From Eq. (9) 
T=e-*', 
Eq. (10) may be written in the form 
@R dr+(1 r—sr)(dR dr)+oR=0, (11) 
where D, o=(8—2w’*s) D. 


Assume R(r) to be a solution of Eq. (11). Then 
a particular solution of (1) is 


c(r, 1) = (12) 


As yet there are no restrictions on 8. However, 
from the nature of the problem certain boundary 
conditions must be satisfied by the function 
c(r, t) which determine definite values of 8. 

The boundary conditions are as follows: (i) the 
net flow through the surface r=/ is zero; (ii) 
the function c(r, t) is an even function in r; 
(iii) when ¢=0 the function c(r, t) represents an 
arbitrary initial distribution of concentration 
which is assumed to be known. The first con- 
dition may be expressed mathematically by 
saying that 


dc for r=l. (13) 


It is now necessary to consider the solution of 
Eq. (11). Change the independent variable from 
r to x by the transformation 


x= 6r?/2. 


The equation which results is 


1 dM a 
(--1 —-—-M=0, (14) 
dx? x dx 
A 
OK—K] 
L 
B 


Fic. 1. Sector-shaped cell. 
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where M is the function of x which corresponds 
to R as a function of r and where 


B=2w*s(1—a). 


Thus the solution (12) may be written in the 


form 


c(r, t) = M(6r?/2)e-)*, (15) 


where 2st. 
Investigation of (14) for a series solution of 
the form 


a,x” 


n=0 
yields the single solution which follows: 
a(a+1) 
1° ‘ 22 


(16) 


ax 
M(x) =1+—+ 


This solution is acceptable in that it is absolutely 
and uniformly convergent for all values of a and 
x and satisfies the second boundary condition 
automatically by being a series in 7’. Let this 
solution be denoted by u(x). A second solution 
v(x) may be found by the substitution® of 


M=u(x)v(x) 


in Eq. (14). On carrying through the substitu- 
tion, the complete solution is found to be 


a fl 2 
Au(x) +5] log 
17 \a 1 


a(a+1) sl 1 22 a(a+1)(a+2) 
12.23 a atl 1 2 12-22. 3? 


1 1 1 22 3 
x(-+— 
a atl at+2 12 3 
where A and B are constants. 

This complete solution cannot be used because 
the factor log x becomes infinite at the origin: 
this requires that B be zero. Thus we are justified 
in accepting (16) as the required solution in 
this case. 

If the expression for c(r, t) given by Eq. (15) 
be substituted in the boundary condition (13), 
the following relation must hold 


(d/dr) M(6r?/2)} —érM(6r?/2)=0 for r=l 


’ Whittaker and Watson, Modern Analysis. 
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or 


(d/dx){M(x)}—M(x)=0 for x=6?/2=b. (17) 


On substituting the series expression for (x) 
in (17) and letting x=) we arrive at the following 
equation 


2a 
(a—1))1+ b+ 
12-2? 


3a(a+1) 
12.22. 32 


=(a—1)f(a)=0. (18) 


Thus one value which a may assume is a=1. 
The other values of a are those which satisfy the 
equation 


f(a) =0, 


which has an infinite number of real, negative 
roots. 

Thus the complete solution of the problem is 
given by 


c(x, An Man(x)e"—*, (19) 


n=0 

where apo, a1, @2, --- denote all the values of a 
which satisfy Eq. (18); Meo(x), Mai(x), Ma2(x) 
-++ denote the values of M(x) which correspond 
to these different roots; @i1, d2, are con- 
stant coefficients which are determined by the 
initial state of the solution in the cell. 

It is now necessary to evaluate the coefficients 
of Eq. (19). For t=0 we have 


Multiply each side of Eq. (20) by o,dx and 
integrate from x=0 to x=0. This gives 


(20) 


b 
f c(x, O)o,dx =do f OnMa,dx 
0 0 
b 


+a f onMa\dx+- 
0 


Each term of the right-hand side of this equation 
is an integral of the form 


I=a, f on Madde. 


If possible, the function ¢, is to be chosen such 
that J=0, for n¥s. If J=0 except when s=n 


(17) 


M(x) 
owing 


(18) 


a=1. 
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then the coefficients ao, a;, --- are given by 


am f ct, / OnMandx (21) 
0 0 


Since Ma, is a solution of Eq. (14) we may 
write 


as dMas d* Ma; 
af =“ (1 —z) +x lax. 
0 


0 a P 


x dx? 


Integration by parts of the right-hand side of 
the above equation yields 


b dMa; 
Os XOn 


0 dx 


b 


d 
— Ma,—(xon) +o,(1 Ma, 
dx 


J On x x. 


If the quantity which multiplies a, under the 
sign of integration be equated to a,¢, it is 
possible to get rid of the integral on the right- 
hand side of Eq. (22), and thus to have the 


integral 
f aM. 
0 


given by functions which are not under any 
integration sign. Thus setting 


(d?/dx*)(xon) (23) 
in Eq. (22) we get 


dMa, 
f onMadx= XOn 
0 


dx 


d 
— Ma;—(xon)+¢,(1 . (24) 
dx 


Suppose the factor ¢, to satisfy the differential 
Eq. (23) in the same manner as the function 
Ma, satisfies the Eq. (14). A very simple relation 
between o, and Ma, is discovered if in Eq. (23) 
we suppose 


Tn 


As a result of this substitution we get the 


equation for s, as 
d*s, dx?+(1/x—1)(ds, dx)—(a, x)s,=0, 


which is the equation satisfied by Me,. Thus 
we obtain as the value of o,, 


o, *Ma,(x). 


Introducing this value of ¢, into Eq. (24) gives 


b 


Man 


1 dMa, 
Mu.) 
As — An dx 


dMan 
Ma( Mon) i (25) 
dx 


But at the point b 
(d Ma,/dx) Ma, = (d Ma,,/dx) Man = 0, 


since a, and a, are values of a which satisfy 
Eq. (17). Thus for the allowed values of a we 


have 
b 


e-* Ma,Madx=0, n#s. 
0 


If s=n the above expression gives for the integral 
an indeterminate form, which will be evaluated 
later. 

From Eq. (21) it is now evident that the 
coefficient a, is given by 


b b 
c(x, O)e-* / J e-*{ May (x) 


Thus the complete solution of Eq. (1) satis- 
fying the boundary conditions is 


f c(x, O)e~*Ma,(x)dx 
c(x, ~ May (x) 
e-*| Max(x)\2dx 
| Ta. (26) 
n=0 


The symbol 7a, is introduced for brevity: its 
definition is obvious. Usually the inifial dis- 
tribution c(x, 0) is a constant ¢o at all points in 
the cell. 
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In the above summation suppose a= 1 as one 
value of a determined by the boundary con- 
ditions. It is also the largest value of a; all the 
others are negative as an inspection of Eq. (18) 
will easily show. Then we can write Eq. (26) in 
the following form, singling out the term for 
which a=1. 


f coe 7M, (x)dx 


0 


de, M(x) +¥. Tan. 
n=l 


b 
f e~*{ }?dx 


0 


But /,(x) and therefore 


coe 7M (x)dx 


M,(x) 
0 
Thus 
bco 
c(x, T,. (27) 
e’—1 n=1 


In this final form 
M(1— 


7 


2 2RT 


M(1—Vo)e* 
RT 


and b= 


It is to be noted that the complete solution con- 
sists of a series of terms, the first term of which 
is the only one independent of the time, and 
which represents the state of the solution when ¢ 
is infinite. The other terms depend on the time 
exponentially and become less as time passes. 

To complete the solution it is necessary to 
evaluate the integrals 


b b 


I,= f Ma,(x)dx and e~*{ Man(x) }*dx. 


Consider first J;. In J;, Ma,(x) is a solution of 
Eq. (14). If Ma,(x)=e7Sa,, Sa, is found to 
satisfy the equfation 


xSan!’ + (1+x)San’—(an—1)San=0, (28) 


ARCHIBALD 


where the superscripts denote differentiation 


with respect to x. 
=> f Seats: 


Thus 
From (28) 
(an—1) f Saude = f 


= x( Say’ + San) ff Sond 


Therefore f Sand x = (x/atn) ( San’ + San) 


b 


and) e-* Ma,dx = 
0 


since (d/dx)Ma,(x)=Ma,(x) when x=). Now 
consider Js. Here it is convenient to consider 
Ma,(x) as a particular case of a more general 
class of functions which are defined as follows, ® 


a(a+1) 

x 
1-2-y¥(y+1) 
According to this definition the function /,(x) 
which has been used all through the above 
discussion is equivalent to \/(a, 1, x). Two simple 


formulae’ will be needed which may be verified 
easily from Eq. (29). These formulae are 


(d/dx)M(a, y, x)=(a/y)M(a+1, y+1, x) (30) 
and 


aM(a+1, y+1, x) 
=(a—y)M(a, y+1, x) +yM(a, y, x). (31) 


To evaluate J. we refer back to Eq. (25). The 
expression 


(29) 


M(a, 7,2) =14+—2+ 
1-y 


dMa,/dx— May, 
which occurs in this equation may be written as 
(d/dx)M(a,, 1, x) — M(as, 1, x). (32) 


By the use of (30) and (31) it can be shown that 
(32) is also equal to 
(a,—1)M(a,, 2, x). 


Thus (25) may be written as 


(a,—1) M(an, 1, 6) M(as, 2, 6) —(a@n—1) 1, M(en, 2, 5)} 


e-*M (an, 1, x) M(a,, 1, x)dx= 


As — Ay 


® Bateman, Partial Differential Equations of Mathematical Physics. 


7 Webb and Airey, Phil. Mag. 36, 129 (1919). 
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--- 720.5 


0.8 
x 


01 02 03 04 05 


Fic. 2. Values of coneentration change ¢/¢o (¢o=in- 
itial constant concentration) as a function of position 
X (=*sr?/2D) for various values of r( = 2w*st) when =0.5 
(=e*sl?/2D). 


When s=1 this gives as the value of the integral 
the indeterminate form 0 0. However the ex- 
pression may be evaluated by differentiating 
both numerator and denominator with respect 
to a, and then letting a,=a,. This gives (since 
M (an, 2, b) =0) 


b 
h=f e~*{ Man(x)}*dx 
0 
= be (an — 1) Man(b)(d/dan)[ M(an, 2, 5) ]. 


Introducing the values obtained for J; and J. 
into Eq. (27) we get the solution in a form which 
is free of definite integrals. 


bco 
c(x, = e* 
+o ~¢ 


An(Qn— 1)(d/dan)[ M(an, 2, b)] 


Recursion formulae as well as a short table of 
values for the functions used have been given by 
Webb and Airey’ and are of considerable service 
for the purpose of computation. 

The curves given in Figs. 2 and 3 represent 
the values of c/co at all points in the cell for 
different values of the parameter r(=2w°s?). 


— 


14 


Ts 1.0 


0.6 x 


02 04 06 08 10 


Fic. 3. Values of concentration change ¢/¢o (¢o=in- 
itial constant concentration) as a function of position 
X (=w*sr?/2D) for various values of = 2w*s/) when b= 1.0 
(=w*si?/2D). 


Fig. 2 is for the case in which b=0.5, and Fig. 3 
is for b=1.0. As an example of the sort of infor- 
mation which the curves yield, consider the 
case of carbon-monoxide hemoglobin for which 
Svedberg® gives the following data: molecular 
weight = 68,000, partial specific volume =0.749 
cm*/g. If a dilute aqueous solution of this sub- 
stance were placed in a cell, such as has been 
discussed above, for which /=0.5 cm and rotated 
at a speed of 385 r.p.s., the value which 6} 
assumes is 0.5 approximately. Fig. 2 shows that 
for this case the equilibrium state has been 
reached to within one percent when r=0.5. 
That is when 2w*s!=0.5. For hemoglobin s=5.0 
X10-" approximately. Thus r=0.5 corresponds 
to a time of 23.4 hours. 

The solution of the above equation for the 
case of a sector-shaped cell placed near the 
periphery of the rotor would be of great interest 
from an experimental point of view. The present 
discussion, although of interest in itself, has been 
given primarily as an introduction to the solution 
for the more general case which has been worked 
out and will appear in a separate article at a 
later date. 


8 Svedberg, article in Colloid Chemistry, Theoretical and 
A pplied, edited by Alexander, Vol. 1, page 851. 


752 


The curves derived from this equation are of 
value for molecular weight determinations. One 
method‘ of determining molecular weights by 
the use of the ultracentrifuge depends upon the 
final equilibrium distribution of concentration 
having been established. As the solution of the 
equation shows, this final state is approached 
asymptotically, it is important to know how 
long a time must pass before it is approximated 
sufficiently accurately. The discussion given 
above of the case of hemoglobin shows that the 
solution of the equation readily yields this 
information if an approximate value of the 
molecular weight can be obtained. 

Further, a comparison of the theoretical and 
experimental distribution curves would be useful 
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for determining whether or not the liquid in the 
cell was free from stirring. If the liquid was 
found to be convection free (as in the case of 
the vacuum type ultracentrifuge® and Svedberg's 
ultracentrifuge’), the equation might make it 
possible to shorten the time of centrifuging by 
rendering it unnecessary to wait for the final 
equilibrium state. A knowledge of the sedimenta- 
tion constant s should allow the molecular 
weight to be determined from the concentration 
curve for any time ¢. 

The writer wishes to express his appreciation 
for the encouragement given him by Professors 
Beams, Snoddy, McShane and Hoxton in this 
work. 


® Beams, J. App. Phys. 8, 795 (1937). 


MAY 1, 1938 


PHYSICAL REVIEW 


VOLUME 53 
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Ryerson Laboratory, University of Chicago, Chicago, Illinois 
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The accidental n-fold coincidences (m >2), are of various 
kinds. The first kind involves m events that accidentally 
occur nearly simultaneously. The second kind involves two 
events that are causally simultaneous, and n—2 events 
that are accidentally nearly simultaneous; etc. For theoret- 
ical purposes, it is sufficient to treat only accidentals of the 
first kind, since the theory of the other kinds is very similar. 
It is also possible, for theoretical purposes, to replace each 
of the various coincidence circuits in actual use by the 
same schematic circuit. This consists of three parts: the 
counting, relay, and recording circuits. For n-fold coinci- 
dences, the relay circuit will contain » relays (e.g. type 57 
tubes). It is shown that the time constants (7;) of these 
relays determine the accidental counting rates. The time 


1. ANALYSIS OF THE PROBLEM, AND 
SCHEMATIZATION OF THE CIRCUIT 


HE increasing precision to which counter 

measurements are being carried makes it 
necessary to consider the statistical theory of 
accidental coincidences in greater detail than 
has hitherto been done. When n counters are 
connected in such a way that only n-fold 
coincidences are recorded, the recording rate is 
the sum of various contributions. Every coinci- 
dence consists of m events that occur simul- 


constants of the counting and recording circuits are not 
important in this connection. It is shown experimentally 
that each of the relays ordinarily has a different time con- 
stant, 7;, even though it is constructed of commercially 
identical parts. By the use of these results, it is shown that 
when JN; particles per second traverse the ith counter, the 
rate at which accidentals of the first kind occur is given by 


A, 2, n= BC, 
Noro: 
++ +1/tn. 
The rates of occurrence of the other kinds of accidentals 
are readily deduced from this result by slight modifications. 
This formula is verified by experiment. 


taneously, or nearly so. If all m events are 
causally related, they form a true coincidence. 
If no two are causally related, they form an 
accidental of the first kind. If two are causally 
related to each other, but no other causal 
relations exist among the events, they form an 
accidental of the second kind. The terminology 
may be extended because there is one kind of 
accidental for every partition of the number n 
into a sum of integers. 

Fortunately, the theories of the various kinds 
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of accidentals are very similar, so that it is 
sufficient if the theory of the first kind is con- 
sidered in detail. The theories of the other kinds 
require only that trivial complications be intro- 
duced. Consequently, the unmodified term 
“accidental”’ will be used to refer to accidentals 
of the first kind. | 

There is general agreement on the formula for 
twofold accidentals: if the rates of the two 
counters are NV; and N2 and the resolving time is 
7, the rate at which accidentals occur is 


(1) 
For threefold accidentals, the rate is 
A 193 = aN,N2N37’, (2) 


but various values of a have been used. The 
correct value is a=3, as was shown by Schiff,! 
but a=4 and a=6 have been used. Not much 
work has required higher values of n. 

The validity of these equations depends on 
assumptions concerning the action of the counters 
and circuit. Although there are several circuits 
in use, it appears that all can be schematically 
represented? by Fig. 1. When one of the counters, 
C, is actuated, it furnishes a pulse of current to 
the corresponding relay. The relay thereupon 
closes and the duration of the contact is r. The 
recorder is actuated only if all relays make 
contact simultaneously. 

There is some question as to what will occur if 
the same counter is actuated twice within a 
time of length r. It is also likely that different 
actual circuits will behave differently, and that o, 
the duration of the pulse of current delivered by 
the counter, will be important in this connection. 
Schiff (p. 94) apparently assumes that or in 


order to avoid this question. This appears to be’ 


a very dubious assumption; in the following 
derivation, the question will be evaded by 
supposing that the counting rates are so low that 
the probability that one counter is actuated 
twice within a time 7 is negligible. This makes the 
efficiency of the relay circuit 100 percent. 

The duration of the pulse of current delivered 
by the relay circuit to the recorder will not be r, 


1L. I. Schiff, Phys. Rev. 50, 88 (1936). 
_ ? An alternative schematization would consist of relays 
in parallel, the recorder being actuated only when all 
relays are open simultaneously. The theory of such a 
circuit is identical with that of Fig. 1. 


but will vary from 0 to 7, depending on the 
timing of the events involved in the accidental. 
(For true coincidences, this duration will always 
be r, unless there should be delays between the 
occurrence of the events in the counter and the 
closing of the relays.) The simplest assumption to 
make concerning the action of the recorder, is 
that all those pulses delivered by the relay, whose 
duration is greater than p, will be recorded. 
Pulses whose duration is less than p are not 
recorded. This limit, p, must be distinguished 
from the resolving time p’ of the recorder: if the 
relay delivers two pulses whose separation in 
time is less than p’, only one of these will be 
recorded. The influence of p’ will again be ignored 
by assuming a sufficiently low counting rate. It 
will be seen that the influence of p is the same as 
if r were diminished to r—p. 

It is doubtful if any actual circuit corresponds 
exactly to the scheme of Fig. 1. The experimental 
investigation described below indicates that the 
idealization is not too great, provided only that 
each relay be supposed to have a different 
contact time, 7;. This is found to be the case 
even when commercially identical parts are used 
in the construction of the relays. 


2. CALCULATION OF THE ACCIDENTAL RATES: 


Let the relays be numbered i=1, 2, ---, and 
let N; be the independent rate of relay i. That 
is, N; is the probable number of events per 
second that close relay i, and are not causally 
related to some event that closes some other relay. 

A coincidence is an interval of time during 
which all relays are closed. Such a concatenation 
of events is diagrammed in Fig. 2. The horizontal 
arrows, 7;, extend from the instant at which 
relay 7 closes until the instant at which it opens 
again. The lengths, 7;, of the arrows are thus the 
contact times of the various relays. The coinci- 
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RELAY 2 
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Fic. 1. Schematic circuit for recording coincidences. 
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Fic, 2. Diagram illustrating a fourfold coincidence. 


dence begins at ¢,, the instant at which the last 
relay (No. 3 in the figure) closes, and ends at ?,, 
the instant at which the first relay (No. 2 in the 
figure) opens. 

The problem is to calculate P(t, ¢’)dt, the 
probability that m independent events fall 
together in such a way as to constitute a coinci- 
dence with ¢<¢,<t+dt and ¢t’ <t,. As a prelimi- 
nary, it is convenient to calculate Q(t, ¢’), the 
probability that the accidental coincidence occurs 
so that ¢4,<¢ and ¢’<t,. The relation between P 
and Q is then clearly 


P(t, t')=9Q(t, t’)/dt. (3) 


If it is permissible to ignore the probability 
that the same relay is actuated twice during its 
contact time, then Q(t, t’)=0 unless ¢’ —t > 1, for 
all relays 7. The probability that relay 7 will close 
during the interval is N,(r;—t'+2), so 
that 

Q(t, t’) = +2). (4) 


Because of the assumption concerning the 
action of the recorder, the rate at which acci- 
dentals are recorded will be P(t, t+). From (4) it 
is obvious that this is the same as P(t, ¢), provided 
that 7; is replaced by r;—p. Hence, if the 7; are 
determined empirically from the recorded counts, 
it is permissible to set p=0. Therefore the 
accidental rate of n-fold coincidences is 


Pit, t) 


= BC, (5) 
where 
B=WN,7: Nore: 
C= 1/ro+ 1/r.. (6) 
For twofold circuits, this reduces to 
(1’) 


and for threefold, to 


A123= Ni N2N3(tit72+ 727347371). (2’) 


If all the 7; are equal, these reduce to (1) and (2), 
with a=3. 

As has been said, these formulae apply only to 
accidentals of the first kind. When threefold 
coincidences are recorded, accidentals of the 
second kind will also occur. The problem is then 
to calculate the total rate, N, at which coinci- 
dences are recorded, given that 
(a) the independent rate of relay 7 is Ni; 

(6) on the average, relays 7 and j are actuated N;; times 
per second by two events that are causally related to 
each other, but are not causally related to an event 


that actuates counter k; 
(c) on the average, true triple coincidences occur N23 times 


per second. 


Let 7;; be the smaller of 7; and 7;; then the 
probable number of accidentals of the second 
kind arising from causally related events in 
counters 1 and 2, and an unrelated event in 
counter 3 is easily seen to be 


A 12, 3= Ni2N3(t12 +753). (7) 
The total rate of recording will thus be 
N= Ni23+A 12, 3+A23, 1+A 31; 2+A 23. (8) 


If all the +; are equal, this reduces essentially to 
Schiff’s Eq. (14). 


3. Discussion OF ACTUAL CIRCUITS 


In order to show that the foregoing assump- 
tions regarding the time constants apply to all 
circuits in use, the various types of circuits will be 
analyzed. It is most convenient to consider all 
these circuits (such as the Neher-Harper, the 
direct resistance coupled, and others) as com- 
posed of three separate circuits: the. counting 
circuit, the relay circuit, and the recording 
circuit. 

A typical onefold circuit is shown in Fig. 3. 
It consists of a Geiger-Miiller tube, an applied 
voltage on the tube, and some arrangement for 
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Fic. 3. Single-fold counting circuit. 
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Fic. 4a. Older relay and recording circuit for coincidences. 


quenching the glow in the tube when a discharge 
occurs. With an activation of the Geiger-Miiller 
tube, a surge of current is produced in the 
counting circuit, resulting in the delivery of a 
pulse of duration o to the relay circuit. This is 
true regardless of the type of circuit used for 
the quenching of this glow. o depends on 
the characteristics and operating conditions of 
the Geiger-Miiller tube, and the constants of the 
remainder of the counting circuit. ¢ is probably 
not only dependent on the constants of the 
circuit, and the characteristics of the Geiger- 
Miiller tube, but also on the voltages, the 
counting rate of the tube, etc. 

The relay circuit, consisting of an amplifier, 


‘acts as a switching device which sets off the 


recording circuit. The plate current is normally 
constant but is reduced nearly to zero by the 
pulse from the counting circuit. The effective 
time during which the current is zero is called r. 
It is easily seen that the duration of the pulse in 
the relay circuit must be at least as great as the 
effective duration of the pulse in the counting 
circuit. (See Part 1.) Since the value of 7 depends 
on the value of ¢, 7 will likewise be affected by the 
factors mentioned above.* 

The recording circuit is controlled by the 
relay circuit. The pulse produced in the latter 
trips the recording mechanism, which in this 
case consists of a thyratron tube and a me- 
chanical count recorder, R. 

Two types of relay circuits have been in use 
for obtaining coincidence records. The first type 
is shown in Fig. 4a where each type 57 tube 
corresponds to one of the relays of Fig. 1. The 


3 In the course of the experiment + was found to vary 
somewhat when the individual counting rate was varied. 
The dependence of r on such factors is being studied at the 
present time. 
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Fic. 4b. Schematization of the same. 


pulse in the recording circuit is due to the 
individual pulses in the relay circuits and is 
equal to their sum. This may be represented 
schematically by an arrangement of switches and 
resistances, as shown in Fig. 4b, where the count 
recorder is so adjusted that it registers only the 
pulses occurring when all relay switches are 
closed simultaneously. Pulses of lesser magnitude 
do not affect the count recorder. This is in 
agreement with the schematization used above. 
The type of circuit shown in Fig. 4a has been 
superceded by the one shown in Fig. 5a. In this 
arrangement the current in resistance Rs; is 
the sum of the currents in all the relays. If the 
current in some of the relays is reduced to zero, 
that in R; is not appreciably changed‘ because 
the remaining relays supply more current than 
before. This is due to the fact that there is a 
slight increase of voltage on the plates of the 
remaining type 57 amplifying tubes and the 
fact that the effective resistance of the tubes is 
not constant. This type of coincidence circuit is 
shown schematically in Fig. 5b. The counting 
mechanism is now so arranged that a count will 
be registered only when the current in Rs; is 
greatly reduced. Hence a coincidence will occur 
only when in an interval of time all the relay 
switches are opened. This too is in agreement 
with the schematic assumption made above. 


4. EXPERIMENTAL VERIFICATION 


The apparatus used in the experimental 
verification of formula (5) consisted of a threefold 
Geiger-Miiller tube counting system as shown in 


‘The following experiment was performed with a 
fourfold counting system: the current in R; was measured 
with a milliammeter, and the currents in the relays were 
successively reduced to zero by putting a negative bias on 
the control grids of the type 57 tubes. It was found that 
the current in R; did not vary by more than 3 percent, 
until all relays were biased, when it dropped to zero. 


| 
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Fic. 5a. Newer relay and recording circuit for coincidences. 


Fig. 3 and Fig. 5a. Each of the three branches 
of this counting system was heavily shielded 
from electrostatic disturbances. 

The Geiger-Miiller tubes were made by a 
special technique developed in this laboratory.’ 
They have a threshold voltage of about 1100, and 
a plateau width of about 800 volts and are filled 
with hydrogen. Their metal cylinder is 9/16 inch 
in diameter and 153 inches in length. Their 
efficiency is high, and the single counting rate, 
produced by cosmic rays and local radiation is 
about 150 per minute when the tubes are not 
shielded, and is decreased to about 60 per minute 
when they are heavily shielded with lead. 

To verify the formula for accidental counts, it 
was necessary to have an experimental determi- 
nation of A123. To obtain this with the ordinary 
circuit constants would be difficult, since with the 
ordinary value of 7 only one accidental in a few 
hours would be expected. 

Furthermore, the cosmic rays would make the 
terms N23, A 12, 3, etc., in (8) much greater than 
A123. It was therefore necessary to increase A 123 


TABLE I. Typical set of data. The A’s and N’s are in 
counts/min., the r’s in minutes. 


N,=182+1 Ai2=16.1+0.2 7, =12.0+0.3 «107 
N2=435+2 Ao3=21.0+0.2 r2= 8320.2 
N;=481+2 A3,=12.0+0.1 t3= 1.7+0.02X10- 


5J. B. Hoag, Electron and ‘Nuclear Physics (D. Van 


Nostrand Company, 1938), p. 432. 
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Fic. 5b. Schematization of the same. 


by increasing the 7; and N;. This was done by 
using abnormally large resistances and capacities 
in the counting and relay circuits, and by 
artificial activation from radioactive sources. 
The quantities N23, Nie, etc., were also decreased 
by shielding two of the Geiger-Miiller tubes with 
2.5 inches of lead, and by placing all three out of 
line in a horizontal plane, at distances of 12 feet 
from each other. In this way, the total number of 
recorded counts was made equal to A 193. 

By means of suitable switches in the circuits, 
the quantities N,, N2, N3, Ai2, Aes, and As; were 
also measured. From these and (1’), values of 
71, T2, and r3 were calculated. These three 
numbers were not equal; a typical set of data is 
shown in Table I. 

These data were obtained with a relay circuit - 
whose homologous parts consisted of com- 
mercially identical resistors, condensers, and type 
57 tubes. When these parts were interchanged, 
the values of 7; altered. Thus the practical 
importance of the generalization of (1’) is 
demonstrated. 

From the values of 7; obtained above, one 
calculates A y23=0.51+0.03 from (2’). The 
value of A 123 observed under these conditions was 
0.513+0.007. Several series of data similar to 
this were obtained ; in no case was the difference 
between the calculated and observed values of 
Ais3 greater than 6 percent. Eq. (2’) was also 
checked by V. C. Wilson for the apparatus 


described by him in a recent paper.® 


®V. C. Wilson, Phys. Rev. 53, 337 (1938). 
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The Magnetic Structure of Cobalt 
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(Received March 7, 1938) 


The magnetic powder method has been used to investi- 
gate the magnetic structure of cobalt crystals at room tem- 
perature. A stable colloid of gamma-ferric-oxide was de- 
posited on suitable cuts of cobalt crystals by a magnetic 
field applied normal to the surface under observation. 
Smooth surfaces were prepared by metallurgical polishing 
followed by additional electrolytic polishing. (The latter 
method for polishing cobalt is described.) The surfaces 
chiefly studied approximated basal or prism planes. On all 
cuts three related patterns were found corresponding to 
no applied field or a normal field applied outward or inward. 
The two patterns with field were reciprocal in the sense that 
their spaces and deposits were interchanged. On basal planes 
the pattern appeared lace-like possessing considerable 
detail which was not easily modified by increasing the ap- 
plied field. On prism planes the no-field pattern consisted 
of straight lines running parallel to the hexagonal axis. A 
small applied normal field spread the colloid in these lines 


into the alternate spaces between them. Further increase in 
field widened the resulting stripes in a somewhat erratic 
manner. Other interesting details of the patterns are 
described. In the discussion are listed the various energies 
involved in the demagnetization of an ideal crystal, i.e., one 
in which all energies concerned are reversible. The attempts 
of other authors to use these energies for devising magnetic 
structures are examined for their relation to magnetic 
powder patterns. It is concluded that the cobalt patterns 
do not accord with the structure demanded by simple 
theory which neglects irreversible energy changes accom- 
panying demagnetization. The structure which seems best 
suited to explain the various patterns of cobalt is then pro- 
posed and discussed, To make this structure appear reason- 
able, the process by which a crystal goes from the fully 
magnetized to the demagnetized state (as suggested by 
various pattern observations) is described in terms of the 
dendritic growth of regions of reversed magnetization. 


HE old magnetic powder method has re- 
cently been used with considerable success 
to establish new data concerning the magnetic 
structure of ferromagnetic crystals. Bitter! who 
first applied the method for this purpose pub- 
lished typical photomicrographs of patterns 
which he found on smooth (but unpolished) 
ingots of iron, nickel and cobalt. Other workers,”:* 
including the present writer, have improved the 
method and have extended the study of the iron 
and nickel patterns. No further investigation of 
the cobalt patterns has yet been reported. 
Interpretation of the cobait patterns in terms 
of a model should be easily made in view of the 
simplicity of the magnetic anisotropy of hexag- 
onal cobalt which has but one direction of easy 
magnetization—the hexagonal axis. The present 
work shows that the structure is more complicated 
than that predicted by simple theory. It has been 
found possible, however, to propose a model which 
will account for the essential features of the 
various cobalt patterns. 
1F. Bitter, Phys. Rev. 38, 1903 (1931); 41, 507 (1932). 
*For a brief discussion with references see F. Bitter, 
Introduction to Ferromagnetism (McGraw-Hill Book Co., 
New York, 1937), p. 59. 
3 For more recent work see W. C. Elmore, Phys. Rev. 


51, 982 (1937) and T. Soller, Zeits. f. Physik 106, 485 
(1937). 


SPECIMENS AND EXPERIMENTAL DETAILS 


Several cobalt specimens were kindly supplied 
by Professor P. W. Bridgman of Harvard who 
had earlier secured them from K. Honda in 
Japan. They were in the form of cylinders about 
one centimeter in diameter and several centi- 
meters in length, containing a number of large 
crystal grains. Various cuts of selected grains 
were prepared by careful sawing followed by the 
usual metallurgical polishing. Many of the sur- 
faces were then electrolytically polished by the 
method devised by Jacquet‘ for copper and brass. 
The orientation of exposed grains was deter- 
mined by etch reflections. The surfaces chiefly 
studied were either approximately parallel or 
perpendicular to the hexagonal axis, the direction 
of easy magnetization in cobalt at room temper- 
ature. 

To polish cobalt by Jacquet’s method an elec- 
trolyte of orthophosphoric acid (specific gravity 
about 1.35) was used. A satisfactory polish 
resulted with 1.2 volts between a comparatively 
large cobalt cathode and the specimen anode. 
Lower voltages, or agitation at 1.2 volts, pro- 
duced an excellent macroscopic etch useful in 


*P. Jacquet, Comptes rendus 201, 1473 (1935); 202 
403 (1936). 
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determining crystal orientations. Somewhat 
higher voltages were found to produce passivity. 
As with copper and brass, a black film appeared 
on the surface at the start of the polishing. For 
cobalt this film proved to be ferromagnetic, as 
will be made clear by the following observation. 

The electrolytic polishing leaves a surface 
not completely flat. If there are parallel ridges 
these are found to coincide with no-field lines in 
colloid patterns (Fig. 3b). This variation in re- 
moval of surface material was at first inter- 
preted as a ‘‘magneto-chemical”’ effect. However, 
a microscopic examination of the surface shortly 
after the appearance of the black film showed it 
to be distributed as in Fig. 1. The relief polish, 
therefore, was due to an interference of uniform 
electrolytic action by the magnetically held skin 
fragments (cobalt oxide?). To avoid this effect 
the black film was carefully wiped away with 
cotton as soon as it could be seen. Since this 
procedure temporarily destroyed in the electro- 
lyte the concentration gradient necessary for the 
polishing action, but did not result in the forma- 
tion of a new film, it seems probable that the film 
was formed during previous exposure to air and 
not in the electrolyte. 

As in previous experiments, magnetic fields 
were applied by mounting the specimen, usually 
irregular in shape, on the pole of a vertical elec- 
tromagnet attached to the movable stage of a 
metallurgical microscope with vertical axis. The 
applied field was thus approximately normal to 
the surface under observation. It was never 
more than a few thousand gauss. Details con- 
cerning the use of magnetic colloid remained 
essentially as described for earlier experiments.® 


PATTERNS ON BASAL PLANE 


The patterns of colloid on the basal plane of 
cobalt or on cuts inclined at small angles to this 
plane were similar to those previously reported 
by Bitter.! However, considerable improvement 
in the detail of the patterns was found to result 
from the use of a stable colloid instead of coarser 
powder which permanently settles out. Upon 
reversing the applied normal field these patterns 
always exhibited a shifting effect similar to that 
of the “‘maze”’ patterns found on polished iron 


5 For a discussion of technique see reference 2, p. 55. 
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crystals. Fig. 2b illustrates the typical pattern 
found with no applied field. An applied normal 
field changed this pattern to that of Fig. 2a or 
Fig. 2c depending on its direction. Careful com- 
parison of the latter two patterns reveals a 
complete reciprocity, i.e., a detailed interchange 
of deposits with spaces. 

The colloid of the no-field pattern is spread 
out more evenly except in numerous small areas 
which are completely lacking in colloid. The 
position of these areas coincides with the centers 
of the meshes of the applied field patterns. The 
colloid of all cobalt patterns was found to be 
much more securely held by stray fields than 
that of patterns on polished iron or nickel crys- 
tals. Indeed, washing the surface with running 
tap water did not seem to disturb it. This pro- 
cedure was found useful for removing excess 
colloid not yet congregated into the pattern. 

In contrast to the “‘maze”’ patterns of iron 
crystals, none of the cobalt patterns required a 
cold-worked surface for their existence. In one 
experiment a mechanically polished surface was 
electrolytically polished to the extent of removing 
a layer 0.1 mm thick. The pattern on the new 
surface was very nearly identical with the old 
one. Patterns on other surfaces electrolytically 
polished to much greater depths possessed the 
same typical appearance. It was observed that 
the patterns on different grains and sometimes 
on different parts of the same grain varied con- 
siderably in scale. 

Except for the shifting effect just described, 
the basal patterns proved to be very stable in the 
presence of the available normal field. This field, 
due to the large demagnetization factor of the 
specimen, was only sufficient to magnetize it to 
perhaps 25 percent of saturation. The stability 
of the basal patterns differs essentially from that 
of the line patterns on prism planes (see next 
section). The only change in basal patterns 
noticed was a sideward motion of some of the 
meshes near the edge of one grain. Bitter has 
shown (Figs. 7-9 of his second paper') what 
happens to the basal patterns on applying intense 
fields parallel to the basal plane. No attempt has 
yet been made to investigate further these 
changes associated with intense fields. For pur- 
poses of discussion it seems correct to consider 
that the present patterns are produced by an 
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Fic. 1. Distribution of film which was set free from the 
surface (a prism plane) early in the process of electrolyti- 
cally polishing cobalt. No applied field. Magnification 
approx. 215X. 


unchanging, or but slightly changing pattern of 
local fields with or without the superposition of 
an applied normal field. 


PATTERNS ON PRISM PLANES 


The patterns on planes approximately normal 
to the basal plane are illustrated in Figs. 3a, 3b 
and 3c. These patterns, whose general appearance 


Bas 


Jk 


a b c 


Fic. 3, Patterns of magnetic colloid on a cobalt surface 
approximately perpendicular to the basal plane with 
applied normal field (a) outward; (b) zero; (c) inward. 
The lines run parallel to the hexagonal axis. Magnification 
approx. 65X. 


does not depend on the angle between their plane 
and a type {1010} plane, were obtained, respec- 
tively, with plus, zero, and minus applied normal 
fields. The no-field pattern consists of narrow 
lines whose general direction is parallel to the 
hexagonal axis. In Fig. 3b the average spacing 
of the lines is about 60g. On other grains different 
spacings were found varying from a few microns 


Fic. 2. Patterns of magnetic colloid on the basal plane of a cobalt crystal with applied normal field (a) outward; 
(b) zero; (c) inward. The spaces and deposits are interchanged in (a) and (c). Magnification approx. 75 X. 
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Fic. 4. Line patterns on several crystal grains with no 
applied field. The lines are spaced closer together in the 
smaller grains. Lines on adjacent grains tend to join at the 

ain boundaries. The V shaped line deposits near these 

undaries are especially to be noted. Magnification 


approx, 65X. 


up to values slightly exceeding that of those 
shown. A variety of spacings can be observed in 
Fig. 4 which shows line patterns on several 
grains. The smaller grains exhibit the closer 
spacings. 

The two patterns with applied field consist of 
bands of colloid located in alternate colloid-free 
stripes of the no-field pattern. They bear the now 
familiar reciprocal relation to each other. The 
width of these bands proved to be very sensitive 
to the magnitude of the applied field. For the 
patterns shown in Fig. 3 this field was the 
minimum necessary to spread some of the colloid 
into bands, though much of it, as can be seen, 
remained in the former no-field lines. The bands, 
however, are somewhat wider than the inter- 
vening spaces. When the field was gradually 
increased still further the bands continued to 
grow in width at the expense of the vacant 
spaces. At the same time the colloid spread itself 
out more evenly over the bands. The edges of the 
widening bands moved with a continuous motion 
punctuated locally by sudden jerks. It appeared 
as if a certain advancing edge would occasionally 
meet an obstacle which for a time prevented its 
progress. After other neighboring boundaries had 
advanced further the impeded boundary would 
abruptly jump into place, and then continue in 
its uniform motion. Two widening bands oc- 
casionally merged either by joining throughout 
their visible length, or quite often by first touch- 
ing only at a point. In the latter case the apices 


Fic. 5. Set of patterns which illustrate widening of bands 
when the applied normal field is increased. (a) No field; 
(b) small applied field; (c) applied field 5 times as intense 
as in (b). Magnification approx. 190. 


of the two V-shaped spaces receded from each 
other with a continued increase of the applied 
field. 

If at any time the field was diminished, the 
motion of the boundaries reversed. Soon, how- 
ever, discontinuities in motion of the sort just 
described would again be observed. Upon com- 
plete removal of the field the no-field pattern was 
similar in scale but in detail apparently unrelated 
to that first observed. Fig. 5 illustrates the 


Fic. 6. Pattern of magnetic colloid associated with 
structure interpreted as a twin band. This pattern was 
found on the crystal grain which gave Fig. 3. Magnifica- 
tion approx, 215 X. 
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widening of bands associated with increasing 
fieid. 

The transition of a prism pattern to a basal 
pattern was studied at the rounded corner con- 
necting the two planes. On approaching the 
corner V-shaped line deposits appeared which 
multiplied the number of lines at the edge of the 
prism pattern. These lines merged with a 
pattern resembling that of Fig. 2 when viewed 
at an oblique angle, the distortion decreasing 
further until the typical basal pattern appeared. 
The increase in the number of pattern lines was 
very similar to that found at the grain boundaries 
in Fig. 4. A somewhat similar group of V-shaped 
deposits was noticed near certain line markings 
found on many of the grains showing prism 
patterns as illustrated in Fig. 6. In this figure the 
narrow band, possessing a fine magnetic structure 
and crossing the black no-field lines, was also 
revealed by etching. This structure may there- 
fore be identified as the trace of a twin band.® 
The plane of such a twinned lamina will in 
general approach the surface making an acute 
angle with it on one side of the trace, thus ac- 
counting for the unsymmetrical occurrence of 
the V-structures which according to this analysis 
depend for their existence upon large demag- 
netizing fields. 


DISCUSSION 


Before discussing the experimental results the 
basic facts of ferromagnetism will be reviewed 
briefly. In this review the attempts of other 
authors to explain the nature of magnetic 
secondary structure will be considered. According 
to the Weiss-Heisenberg theory,’ a ferromagnetic 
crystal will be at all times spontaneously mag- 
netized to a saturation value which depends on 
the temperature and becomes zero at the Curie 
point. Actual crystals, however, possess little or 
no net magnetic moment unless placed in a mag- 
netic field. To resolve this apparent contradiction 
the crystal has long been considered to be sub- 
divided into numerous small domains each 
saturated, so that the crystal as a whole appears 


° The twinning of cobalt has apparently not been re- 
ported. The present twinned regions are tentatively 
identified as occurring on the {1012} pyramidal planes, as 
found for other close-packed hexagonal crystals. 

7E. C. Stoner, Magnetism and Matter (Methuen and Co., 
London, 1934), Ch. XI. 


2 d, 


Fic. 7. (a) Schematic diagram (according to Landau and 
Lifshitz) showing arrangement of magnetization in plane 
perpendicular to a surface which is perpendicular to direc- 
tion of easy magnetization. (b) Previous diagram so 
modified as to include the presence of an applied normal 
field H,. Neither of these arrangements of magnetization 
demanded by minimum total energy of the crystal gives 
rise to superficial stray fields which can produce magnetic 
powder patterns. 


unmagnetized. This division into domains can be 
most simply dealt with by considering the various 
energies involved in the absence of applied fields. 


(1) An energy density £,, associated with the 
Weiss field (exchange interaction) 


Ew is increased by changes in direction of the 
local magnetization I,. Explicitly 


Ew= —3NI 2+ 


where JN is the Weiss field constant, ‘‘a’’ is the 
lattice separation of nearest neighbors and I, is 
the intensity of magnetization (varying only in 
direction). The first term is well known.’ The 
second term, used by Landau and Lifshitz* in a 
manner to be mentioned presently, may be ob- 
tained readily by averaging the energy of inter- 
action of a particular elementary magnet with 
its nearest neighbors. This interaction is assumed 
for each pair to be of the form —clI,-I,’ where 


in which Ar gives the position of the selected 
neighbor. Here c=3N since E,= —}NJ,? when 
I, does not vary in direction. In what follows Ey 
will be taken to stand for the second term alone. 


8 L. Landau and E, Lifshitz, Physik. Zeits. Sowjetunion 
8, 153 (1935). 
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(2) An energy density E, associated with the 
magnetic anisotropy of the crystal 

This energy may include the effects of strain 
as well as the natural crystal anisotropy. FE for 
unstrained cobalt at room temperature is found 
empirically to be represented fairly well by 
Ey=K' sin? 0+K” sin‘ @ where @ is the angle 
between the hexagonal axis and I,, K’=5.1 10° 
and K’ =2.2X10°.° The behavior of there- 
fore, depends essentially on the quadratic term, 
and the approximation Fy+ 38/,? sin? @ may be 
conveniently used in discussing magnetic struc- 
ture. In the simplified expression 8 may be taken 
to be 2(K’+K”") /72=7.2 so that is correct 
for both @=0 and 6=90°. 


(3) An energy density E,,= —}H-I, where H is 
the macroscopic field due to the entire 
crystal calculated from the potential 


I-de 
r 

Frenkel and Dorfman" seem to have made the 
first approximately correct estimate of the 
diameter dy) of magnetic domains. By a simple 
argument (neglecting Ey) they are able to es- 
timate that E,dr=C,/d and that Jj" 
=(C.d for a space-filling array of disorientated 
domains of diameter d within a crystal of volume 
V. It follows that do=(C,/C2)! for minimum 
energy. They find for iron crystals of moderate 
volume V that dy= 100x. 

Bitter," working with a torque equation equiv- 
alent to minimizing /(E,.+£,)dr, integrated 
over regions of the crystal remote from the sur- 
face, has discussed what he terms the ‘elastic 
deformation of ferromagnetism.”’ The particular 
cases which he studied are unfortunately chosen 
since in all of them V-1,#0. The torques due to 
the resulting internal magnetic fields are suf- 
ficient to render his arrangements of I, physically 
unstable. His attack on the problem, if carried 
far enough, should lead to the more detailed 
conclusions of Landau and Lifshitz* which will 
now be discussed. 

Landau and Lifshitz minimize /(E,,+s)dr 
integrated over the entire crystal. They conclude 


*J. H. Van Vleck, Phys. Rev. 52, 1178 (1937). All 
energy densities are expressed in erg. cm~*. 

10 J. Frenkel and J. Dorfman, Nature 126, 274 (1930). 

uF, Bitter, Sci. Rep. Téhoku Imp. Univ. (Honda 
Volume), 228 (1936); also reference 2, p. 186. 
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that in the demagnetized state the magnetization 
should occur chiefly in plane layers magnetized 
alternately parallel and anti-parallel to the 
direction of easy magnetization. Near a surface 
perpendicular to this direction they propose the 
arrangement of I, shown in Fig. 7a, for which 
So’ Endr=0 since V-I,=0 and J,=0. The 
transition of I, from one layer to the next is 
found to occur gradually, not abruptly. It takes 
place by a twisting of I, about the normal to the 
boundary between layers, thus making V-I,=0. 
Their equations indicate that the thickness of 
the transition layer for cobalt is about 45 lattice 
spacings. The thickness do of the layers is found 
by minimizing the sum of the energy of the 
transition layers (of the form E;=C, d) and the 
anisotropy energy in the regions of triangular 
cross section (of the form E,=Ced). If / is the 
distance between two basal surfaces of a cobalt 
crystal they find that dy=(8/)'(a@/8)* for mini- 
mum energy, where a= }a?N=8.3X10-" and 
8=7.2. This gives dop=15u for 1=0.3 cm, the 
approximate size of the larger crystal grains in 
the specimens here studied. This spacing is less 
than that of the parallel lines of the prims 
patterns. It is more nearly a measure of the 
scale of the finer structure in the basal patterns. 
The minimum energy /9"(E,,.+£e)dr is some- 
what increased if instead of layers, prisms of 
square or triangular cross section are assumed, 
together with the arrangement of I, near plane 
surfaces which makes Eydr=0. Hexagonal 
prisms cannot be arranged so that JS," Endr 
vanishes, hence a model with these units would 
have a still greater minimum energy. 

Landau and Lifshitz have also considered 
crystals with surfaces inclined to the basal plane. 
By suitable arrangements of I, they retain the 
condition that /9”Eydr=0 and calculate for the 
layer model the dependence of do upon the in- 
clinations of the exposed surfaces. If their con- 
dition that V-I,=0 within the crystal be applied 
to the arrangement of Fig. 7a when an applied 
normal field is present, the situation becomes 
that depicted by Fig. 7b. For this arrangement 


V-I,=0, J,=I], cos 6=H,/47, 


The change in scale of the structure indicated by 
the last equation results from the decrease in £, 
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due to: (1) a decrease in volume of the triangular 
prisms seen end-on in Fig. 7; (2) a decrease in 
E, in these prisms resulting from the rotation of 
I, therein. From the equation it appears that d)’ 
differs appreciably from dy only for large J,,. The 
difference in scale is important, however, in con- 
sidering the regrowth of magnetic structure in a 
crystal which has been magnetized to saturation 
by an external field. 

It is a characteristic of the crystal internally 
divided into layers or into square or triangular 
prisms, that /j"Eydr=0. Hence, in the absence 
of an applied field, a crystal possessing one of 
these structures should give no pattern. Even an 
applied normal field should produce no pattern 
if Fig. 7b is correct, for the field above the 
specimen is uniform although the magnetization 
beneath possesses a definite structure. In most 
cases, however, some stray field will be present ; 
the result of an effect neglected in the Landau- 
Lifshitz model. I, in the surface prisms of tri- 
angular section will in general rotate so as to 
diminish (Z,+£y). This requires Ey, >0. This 
effect can be estimated for the arrangement of 
I, in Fig. 7b, as follows: replace the triangular 
prisms by an extensive parallel slab of crystal 
initially magnetized parallel to its surface. By 
minimizing the energy E= Ey+ Ey, cos 6 it 
is found that a normal field //,, will tilt I, by an 
amount given by cos @=//,,/[(42%—8)/,] rather 
than by cos As for cobalt, 
this effect is quite appreciable, hence in the 
event of a layer structure, some deposition of 
colloid should be expected. 

For ideal crystals there seems little doubt that 
the treatment just outlined is essentially correct. 
The discrepancy between the ideal structure and 
the one indicated by the patterns probably 
arises during the process whereby the crystal 
goes from the magnetized to the demagnetized 
state. Irreversible energy losses accompanying 
this process dissipate some of the potential 
energy which would otherwise be available to 
establish the ideal structure. Nevertheless, the 
Landau-Lifshitz treatment is important since it 
justifies the description of magnetic structure in 
terms of geometrical models. Furthermore, it 
suggests arrangements of I, likely to occur, such 
as that of Fig. 7a. For instance, if the layers 
should extend all the way to the surface, the 


resulting surface energy would be 
much greater than that when the triangular 
prisms lying in the surface afford a return path 
for lines of induction within the crystal. With this 
background the patterns can now be analyzed. 

The prism patterns (parallel lines and bands) 
can be interpreted as indicating sections of 
regions magnetized alternately parallel and anti- 
parallel to the hexagonal axis. The no-field 
patterns indicate that the surfaces so far ob- 
served have not been exactly parallel to the axis 
so that alternate strips have 7, plus and minus. 
This accounts satisfactorily for the observed 
patterns. The widening of the bands with in- 
creasing applied field is, of course, an additional 
feature. The component of the applied field 
parallel to the axis must be responsible for this. 
Thus by a suitable choice of field direction near 
the normal it should be possible to control inde- 
pendently the width of the bands and the par- 
ticular set of strips preferred by the colloid. 

Since the field H below the surface arising 
from the potential {I,-de/r is an image of the 
stray field above the surface, the magnetization 
near the surface will be rotated so as to decrease 
the stray field. This decrease amounts to perhaps 
60 percent for surfaces slightly inclined to the 
hexagonal axis. It should not alter the essential 
features of the prism patterns. 

The spacing of the parallel lines requires some 
comment. Its regularity is probably not typical 
of the magnetic structure existing before the 
surface was prepared. For surfaces not too 
oblique to the axis the number of regions should 
not be changed, but the energy will be minimized 
by making the areas of adjacent regions of 
opposite polarity more nearly equal. 

A clear interpretation of the basal patterns is 
more difficult to achieve. As already pointed out 
none of the simple structures of layers, square or 
triangular prisms is adequate. A clue is found in 
the fact that the spacing of the parallel lines is 
about equal to the spacings of the irregular areas 
of light and heavy deposit in the basal patterns. 
It is suggested that the more or less definite 
boundaries between these distinguishable areas 
are traces of the boundaries separating regions 
of reversed magnetization deep in the crystal. 
The finer detail seems closely connected with 
dendrite-like reversed regions extending short 
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distances into the crystal and producing in cross 
section the V-shaped line deposits. This picture 
is probably further complicated by the occurrence 
of dendrites within dendrites, as some of these 
deposits in Fig. 4 suggest. The following analysis, 
suggested by the behavior of the cobalt patterns 
studied so far, is an attempt to make this mag- 
netic structure appear a reasonable one. 

Imagine a disk shaped crystal having basal 
planes as faces and magnetized to saturation by 
a field applied parallel to the disk axis. If H is 
reduced slightly below 47/, regions of reverse 
magnetization must form. Presumably many of 
them occur first at the two basal surfaces and 
then grow inward towards one another, finally 
meeting by pairs to form complete thread-like 
regions extending through the crystal. In their 
early stages of growth these regions resemble 
dendrites of a solid phase growing into the melt. 
The conditions are similar in respect to the 
gradual supply of energy necessary for the change 
occurring and in respect to the stability in space 
of the transformed regions. 

As /7 is still further reduced, more dendrites 
form, those already present continue to grow 
inward and the completed threads (or prisms) 
increase in cross section. The surfaces of the disk 
will on the average be uniformly populated with 
reversed regions. The number of reversed regions 
will become greater as demagnetization proceeds. 
This corresponds to the change in scale of the 
ideal layer structure, previously pointed out. 
The requirement of finer scale at later stages 
necessitate that some of the new dendrites form 
within regions already anti-parallel to H, and 
hence be themselves parallel to H. Photographs 
have been taken of V-structures near the bound- 
aries of prism patterns which, if one judges from 
the deposit on them, were magnetized parallel to 
H, and which grew as H was decreased. Finally, 
when HT becomes zero, many of the dendrites 
started last will remain as such, there being in- 
sufficient energy available for them to grow far 
into the crystal. This is equivalent to stating 
that the residual magnetic field at the boundaries 
of the dendrites is insufficient to cause their 


further advance. The resulting magnetic structure 
near a basal surface is not one of plane parallel 
layers, but one more nearly described as consist- 
ing of hexagonal prisms of about the same size, 
arranged so that the energy /o"(E£.+4:+Ey)dr 
is as small as possible. A short distance from the 
crystal surface many of these prisms pinch out 
leaving the much coarser structure responsible 
for the typical prism patterns. 

Some of the hexagonal prisms will be com- 
pletely surrounded at a basal surface by others 
of opposite magnetization. An arrangement of I, 
corresponding to the layer model in Fig. 7a is 
demanded by minimum energy for such a 
grouping. Hence there will be locally no stray 
field above the surface. The absence of colloid in 
certain areas of the no-field basal patterns, and 
the very small amount appearing above these 
same areas in the normal field patterns may thus 
be accounted for. At other boundaries between 
oppositely magnetized prisms only a part of the 
lines of induction can return within the crystal. 
Indeed, the best possible arrangement of regular 
hexagonal prisms permits the return of but two- 
thirds of the lines of induction in this manner; 
the remainder give rise to the stray field which 
is chiefly responsible for the basal patterns. 

In conclusion it should be pointed out that 
the erratic movement of the parallel lines is 
visual evidence of the Barkhausen effect. This 
effect in cobalt does not seem to consist of the 
sudden reversal of definite domains in the crystal, 
but may be ascribed to the uneven motion of 
boundaries between regions considerably larger 
than those usually assigned to account for 
Barkhausen discontinuities, containing about 
10" atoms. Local imperfections in the crystal are 
no doubt responsible for this effect. They also 
must be largely responsible for energy lost 
during the process of demagnetization. 

The author is greatly indebted to Professor 
F. Bitter of the department of metallurgy for 
providing certain facilities for conducting this 
research, and to Professor L. W. McKeehan of 
Yale for reading and discussing the manuscript. 
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Prompt publication of brief reports of important discoveries in physics may be secured by 
addressing them to this department. Closing dates for this department are, for the first issue of the 
month, the eighteenth of the preceding month, for the second issue, the third of the month. Because of 
the late closing dates for the section no proof can be shown to authors. The Board of Editors does 
not hold itself responsible for the opinions expressed by the correspondents. 


Communications should not in general exceed 600 words in length. 


Radio Isotopes of Nickel 


2.6-hour Ni®* 

The early observations':? of a radioactivity of a few 
hours half-life induced in nickel by slow neutrons has been 
rather definitely identified as due to Ni*® through the work 
of Oeser and Tuck* with slow neutrons and that of Heyn* 
with fast and slow neutrons. It has also been produced by 
the bombardment of copper with neutrons,’: 4 and of zinc 
with neutrons.‘ Chemical identification has been made in 
all these cases. The half-life is quoted as 2.5 and 2.66 hours. 
It is rather surprising, therefore, that this period should 
not have been found by Thornton’ after bombarding nickel 
with 5 Mev deuterons, since its presence is to be expected, 
if due to either Ni® or Ni®. Only a positron emitting isotope 
of copper was found, presumably Cu”, with a half-life of 
3.4 hours; the activities in the nickel and cobalt precipi- 
tates were disregarded, as their intensities were only 0.01 
and 0.005 that of the copper. 

This note is to point out that the deuteron bombardment 
of nickel does produce the nickel activity as well as that 
of the copper (and in addition a cobalt period of some 
months’ half-life on which we shall report later). The 
exposure of nickel for four hours to 100 microamperes of 
5.5 Mev deuterons yields a chemically identified nickel 
isotope of half-life 2.6+0.03 hours; the emitted particles 
are negative electrons which are absorbed to half-value 
by 0.13 gram/cm? Al and which have a range of 0.90+0.1 
gram/cm? Al, indicating a maximum energy of 1.9 Mev, 
by Feather’s relation.* The gamma-rays are absorbed to 
half-value by 10.8 grams/cm? Pb, equivalent to an energy 
of 1.1 Mev, according to the data of Gentner.? Heyn’s 
value for the half-thickness of the electrons from the 
neutron-produced activity is 0.10 gram/cm*. The identi- 
fication with this previously known period seems to be 
complete. 

We cannot give an exact estimate of the intensity ratio 
of the Ni® and Cu" activities, since only part (perhaps 
one-quarter) of the voluminous nickel dimethylglyoxime 
precipitate was followed. The extrapolated initial intensity 
of this fraction of the nickel was 3.4 X 10° times background, 
while that of the entire copper precipitate was 4.0X 10° 
times background (approximately 4 millicuries). 

The production of this 2.6-hour nickel through the re- 
action Ni®(d,p)Ni® may alter the interpretation of the 
excitation function for Cu® produced by Ni®(d,n)Cu® as 
given by Thornton, since a chemical separation was not 
performed on the stack of nickel foils used to determine 
the activation probability. It is difficult for us*to under- 


stand why the gamma-ray from Ni® was not observed by 
Richardson® in his examination of the gamma-rays from a 
deuteron-on-nickel bombardment, unless a chemical sepa- 
ration for copper and nickel was made, which is not so 
stated. 

Our neutron bombardments of nickel have not disclosed 
any nickel periods shorter than the 2.6-hour activity; 
short periods from deuteron irradiations of nickel have 
not been looked for. 


36-hour Ni*™ 

We wish to report a new radioactive isotope of nickel, 
formed as the result of the exposure of iron to several 
microampere hours of bombardment with helium ions 
at 12.6 Mev and also at 16 Mev. The chemically separated 
nickel emits positrons with a half-life of 36+2 hours. 
Absorption measurements on the particles indicate a 
half-value thickness of 0.034 gram,cm? Al and a range of 
0.2540.05 gram/cm? Al, corresponding to a maximum 
energy of 0.67+0.1 Mev by Feather’s relation.* A gamma- 
ray is also evident, probably annihilation radiation. 

The only unstable nickel isotopes that could be formed 
from iron by the Fe(a,n)Ni reaction are and 
The lightest stable nickel is Ni®* (abundance 68 percent), so 
that if the 36-hour period were due to Ni’, it should be pro- 
duced by slow neutrons or deuterons on nickel. We have not 
been able to detect this activity after strong irradiation of 
nickel with deuterons or slow neutrons, so we feel justified 
in ascribing the activity to through Fe®(a,n)Ni°’, 
followed by positron decay to stable Co’. With fast 
neutrons (from lithium plus deuterons) on nickel, we have 
found a weak and somewhat doubtful indication of the 
period, intermixed with some longer lived activities, so 
that possibly the reaction Ni®*(n,2n)Ni*’ has also been 
observed. 

We are glad to acknowledge the continued assistance of 
the Chemical Foundation, the Research Corporation and 
the Josiah Macy, Jr. Foundation. 

J. J. Livincoop 
G. T. SEABORG 


Radiation Laboratory, Department of Physics, (j.J.L.) 
Department of Chemistry, (G. T.S.) 
University of California, 
Berkeley, California, 
April 14, 1938. 
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Drift of Adsorbed TH on W Filaments Heated with d.c. 


A study (soon to be reported) of the surface structure 
produced on tungsten lamp filaments during long heating 
with direct current has suggested that W ions drift over the 
wire surface toward the negative terminal. Such migration 
can be demonstrated directly, for the case of Th adsorbed 
on W, in a simple electron-optical tube. 

A loop of 1 mil thoriated wire 6 mm long is mounted on 
30-mil leads at the center of a spherical bulb 3 inches in 
diameter coated inside with willemite. The anode, a wire 
ring near the mouth of the bulb, is made a few thousand 
volts positive, and the emitted electrons form on the screen 
a magnified orthogonal image of the loop; the screen is kept 
near anode potential by secondary emission. A 10-mil Ni 
wire, welded to one filament lead and projecting between 
the ends of the loop, casts an electrical shadow which 
permits distinguishing the emission from the two ends, and 
also prevents evaporation of Th from one end across to the 
other. With the anode floating, the filament is first flashed 
and then is heated, from a 2 v d.c. supply, for 20 or 30 
minutes at a low activating temperature. The electron 
image, observed briefly at a low temperature where Th is 
immobile, is then brightest opposite the negative end of the 
loop, with a sharp cutoff toward the lead and a gradual 
decrease toward the center. The positive side of the screen 
is dark. With 60-cycle a.c. heating, the activation is always 
symmetrical about the center of the wire. These observa- 
tions have been repeated in a similar tube, with a loop of 2 
mil thoriated wire. 

The processes involved in activation of a short filament 
are obviously manifold, but the net result, that Th 
accumulates preferentially at the negative end, strongly 
indicates a drift of ions in the electric field. Ionization in the 
space is ruled out by the low voltage of the heating supply, 
the shadow-caster stops evaporation from one end to the 
other, and thermionic emission from either end is at the 
worst limited by space charge to a value too low to give 
appreciable cooling. 

This experimental arrangement is not adapted for 
yielding values for either the mobility or the activation 
energy for migration under the influence of a field. Deposit 
of Th on an isothermal length of smooth undoped filament 
(preferably a single crystal), by evaporation from an 
external source, at present seems a promising technique for 
future study of the effect. 


Research 

General Electric Company, 

Schenectady, New York, 
April 6, 1938. 


R. P. JoHnson 


On the Form of the 8-Spectrum of Ra E in the Vicinity of 
the Upper Limit and the Mass of the Neutrino 


The spectrum of the 8-particles emitted by Ra E has 
often been studied and at the present time! it is established 
that the form of the spectrum over a wide range may be 
accounted for satisfactorily by the formula of Uhlenbeck 
and Konopinski. However the spectrum in the immediate 
neighborhood of the upper limit has been studied in much 
less detail. Nevertheless the results of Lyman’s? measure- 
ments and those described in this note indicate that the 
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form of the spectrum in this region departs considerably 
from that predicted by theory. The study of the energy 
distribution of the electrons in the neighborhood of the 
spectral limit is for a number of reasons accompanied by 
considerable difficulty, chief of which is the low intensity 
of the radiation in this part of the spectrum: The choice of 
apparatus and the sensitivity of the method of registering 
the electrons are therefore of prime importance. In the 
experiments described below, the method of electron focus- 
ing in a uniform magnetic field has been used. The electrons 
were registered by the method of coincident discharges in 
Geiger-Miiller counters. It was necessary to guard against 
electrons scattered from the walls of the apparatus and 
from the diaphragm entering the counters. This was 
attained by separating out a narrow beam of electrons 
(@=10°) and by a special system of traps. 

By placing the first counter with its narrow slit at the 
focus it was possible to eliminate almost completely the 
entry of scattered electrons. The source was a thin nickel 
strip on which 0.3—1.0 mC of Ra E were deposited. The 
spectrum was studied in detail from E=1030 kv to 
En=1180 kv. The results are represented in Fig. 1(I) 
and it will be seen that the energy intervals were very 
narrow. The random error of the individual points on the 
curve were very small and the measurements of the 
magnetic field constitute the main error. A noticeable 
feature in Fig. 1 is the sharp change in the curve at an 
energy of about 1120 kv after which the number of elec- 
trons decreases practically linearly. This part of the curve 
is reproduced on a larger scale in the same figure for the 
sake of clarity (V). The results of repeated experiments are 
indicated by different signs. The insignificant number of 
electrons observed at E>E, is undoubtedly due to 
scattered electrons and indicates that the experimental 
conditions were satisfactory. 

The spectra computed according to the formula: 
P~W-(W?—1))-(Wo— W)? 

— (u/m)*}'-f(Z, W). 


which follows from the theory of Konopinski and Uhlen- 
beck are given in Fig. 1. Curve II refers to the case ».=0 
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and Wo=3.69moc? and curve III to the case n=0.3mo and 
Wo=3.69moc*. It is clear that curve III agrees excellently 
with the experimental results. It may be noted that the 
remainder of the experimental spectrum also fits the Uhlen- 
beck-Konopinski theory satisfactorily when « =0.3mo. The 
curve IV of Fig. 1 has been computed according to Fermi’s 
formula, from the experimental value of the spectral limit 
(Woexp =3.31moc? while » was taken equal to zero. This 
curve clearly does not agree with that found experimentally. 
A. I. ALICHANIAN 

A. I. ALICHANOW 


Oy 
B. S. DZELEPow 


Leningrad, U.S 
March 16, 1938. 


1A. Alichanian and A. Zavelsky, Comptes rendus d. Acad. d. 
L'URSS. 17, 467 (1937). 
2 Lyman, Phys. Rev. 51, 1 (1937). 


The Shape of the 8-Spectrum of Th C and the Mass 
of the Neutrino 


In the preceding letter of A. I. Alichanian, A. I. 
Alichanow and B. S. Dzelepow it was shown that the 
intensity curve of RaE cuts the axis of abscissa at a 
certain angle. 

It was of importance therefore to investigate whether 
the shape of the Ra E spectrum is typical for all radio- 
active elements. The end of the spectrum of Th C was 
accordingly studied by the magnetic focusing method with 
registration by means of coincidences. The apparatus was 
the same as that used for Ra E.' A narrow beam of elec- 
trons with an angle of spread ¢=10° was focused onto the 
slit of the first counter (d~0.8 mm). A considerable 
number of diaphragms and traps prevented scattered 
electrons from entering the counting system. An active 
deposit of Th(B+C+C”) (1.0-1.5 mC) was deposited on 
an aluminum strip 0.3 mm wide and 5 mm long. 

The end of the spectrum is represented in Fig. 1(I). 
There is a distinct break in the curve at E=2130 kv, 
similar to that observed in the spectrum of Ra E, after 
which the curve falls off practically linearly to the very 
limit. The small number of coincidences observed for 
H>Hm is due partly to Compton electrons liberated 
from the walls of the vessel by y-rays, and partly to 


+ exp 
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scattered electrons and random coincidences. It is clear, 
however, that this background does not depend on the 
strength of the magnetic field. The remainder of the Th C 
spectrum down to E=1275 kv is given in the same figure 
(V). By constructing the Konopinski-Uhlenbeck plot from 
this spectrum (assuming the mass of the neutrino »=0) 
we found the extrapolated limit Wo=6.18 with the help 
of which the K—U spectrum for 4=0 was constructed 
(curve II). Curve III was computed according to Fermi’s 
theory with n=0 from the experimentally determined 
spectral limit for the total energy of disintegration. Finally 
curve IV has been computed according to the Konopinski- 
Uhlenbeck theory, for 4 =0.8»10 from the formula 


P~W(W?—1))-(Wo— W)8-[(Wo— W)?— (u/m)?})-f(Z, W) 


where Wo=(Wo)exp+0.8moc?. It is clear that none of the 
theoretical curves cited, with the exception of curve IV, 
agrees with the experimental spectrum. The satisfactory 
agreement between the latter and the Konopinski-Uhlen- 
beck curve for 1 =0.8mp was not confined to the immediate 
vicinity of the limit but extended to the rest of the spec- 
trum. A similar analysis in the case of the Ra E spectrum 
showed that agreement between the theoretical and experi- 
mental curves is obtained when un =0.3mpo. 

The following conclusions may be drawn from these 
results. The existing theory of Konopinski and Uhlenbeck 
describes the experimental results correctly when »+0. 
u however is not a constant for all elements and varies 
with W’o. 

It was found that the most suitable value of « could be 
deduced from the equation: 

(Wodextr— (Woexp =p 
which gives 0.8moc? for Th C and 0.3moc* for Ra E. The 
fact, however, that uw varies from element to element 
means that no simple significance can be attached to this 
parameter. 

Preliminary measurements of the end of the spectrum 
of Ra C showed that in this case also the shape is similar to 
that in the vicinity of the spectral limit of Ra E and Th C. 

The three maxima observed in the figure at FE; = 2526 kv, 
E,=2588 kv and E;=2610 kv are due to internal con- 
version of the y-rays of Th C” (hv=2614 kv) in the K, L 
and M levels of Th Pb. The sharpness of the lines bears 
witness to the high resolving power of the apparatus and 
the absence of disturbing factors. The relative intensities 
of the lines 5: 1:0.2 give the relative probabilities of 
conversion of the y-line 2614 kv in the K, L and M levels. 
The small maximum at E=2070 kv corresponds to the 
internal conversion of the y-line hv=2160 kv found by 
Alichanow and W. DZelepow? in the positron spectrum of 
Th (C+C”). The present measurements.yield a value of 
a=1.77 X10-* for the coefficient of internal conversion of 
the y-line 2614 kv whereas according to Taylor and Mott® 
the theoretical value is a=1.5 107%. 


Physico-Techaical 
Leningrad, U.S.S.R., 


A. I. ALICHANIAN 


March 16, 1938. S. J. NikITIN 
a A. Aiheaten, II Union Conference on Nuclear Physics, Moscow, 
pt. 27, 
2A. and W. P. Dizelepow, Physik. Zeits. Sowjetunion 
in print). 


3H. M. Taylor and N. F. Mott, Proc. Roy. Soc. 138, 665 (1932). 
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Note on Atmospheric Absorption Caused by the Rotational 
Water Band 


In a recent article! the writer computed the infra-red 
absorption of the atmosphere which is caused by the 
rotational water band. The work was based upon the 
analysis of the rotational spectrum given by Randall, 
Dennison, Ginsburg and Weber.? The average absorption 
for an interval containing many individual lines is pro- 
portional to the square root of the optical mass, a known 
result for the absorption of a line spectrum when the 
distance between the lines is large compared with the line 
width. 

The black-body curve for T=300° extends from about 
4u towards longer waves and has its maximum around 10u. 
Thermal radiation of the atmosphere is mainly due to 
the vibrational water band at 6u and to the rotational 
band. It appears from the above-mentioned calculations 
that the highest rotational lines are too weak to cause 
appreciable absorption and with an average water content 
of the atmosphere there should be practically complete 
transparency between about 84 and 22u. This result is 
not in agreement with the measurements of the heat 
radiation of the sky which give much higher values of 
the total radiation than one would expect on this basis. 
The discrepancy can be removed by taking into account 
the effect which is produced by the very intense rotational 
lines at a large distance from the line center. Although 
these lines are situated at the far infra-red side of the band 
(50u), they produce absorption between 10 and 20u. The 
absorption coefficient due to the line wings is readily 


seen to be 
k =Liva/e(v— vo)’, 
vo 


where the summation extends over all lines of the band. 
The io are the line intensities and a is the half-width for 
which, as previously,! a value of 0.25 cm was assumed. 
Under atmospheric conditions line broadening is due to 
collisions and therefore @ is proportional to the air pressure. 
The above value refers to standard pressure. Numerical 
calculations have been carried out for T=220° (tempera- 
ture of the stratosphere) and T=300°. They yield the 
results given in Table I, if the optical thickness of water 


TABLE I. Values of k at 200 and 300° Kelvin. 


v(em~*) | 400 450 500 600 700 800 900 1000 


1.43 0.96 0.54 0.35 0.23 0.18 0.14 
1.36 0.90 0.46 0.30 0.21 O15 0.12 


2.61 
3.50 


Roo 


vapor is expressed in grams of water per cm?. For »>500 
cm~ the following formula gives a good approximation 


k = (78,000 — 2402) /(»— 200)? 


where ¢ is the temperature on the centigrade scale. It 
must be emphasized that this absorption is exponential, 
J=Joe-**; it is therefore far more sensitive towards 
changes in the water content of the atmosphere than the 
absorption within the band itself which varies only with 
the square root of the water mass. This is in agreement 
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with the observations of the thermal radiation of the sky 
which shows a very considerable dependence upon hu- 
midity. Under laboratory conditions the continuous 
absorption is too weak to be observed in moist air: it 
does however show up in the absorption of steam.’ The 
absorption is proportional to the half-width of the lines 
and since the latter is not precisely known, the relative 
variation of the above figures is more significant than 
their absolute values. 

A more detailed account of this work will be given later. 

W. M. Evsasser 


California Institute of Technology. 
Pasadena, California, 
April 15, 1938. 


1W. M. Elsasser, Astrophys. J. April, 1938. 

2H. M. Randall, D. M. Dennison, N. Ginsburg and R. L. Weber, 
Phys. Rev. 52, 160 (1937). 

3G. Hettner, Ann. d. Physik 55, 476 (1918). 


Recalescence in Uranium 


During the summer of 1937 the thermionic work func- 
tion and the emissivity of well-outgassed uranium were 
investigated. The results will be published at a later date. 

In the course of these experiments a striking double 
recalescence was observed, resembling the well-known 
phenomenon in iron at the 6-y transition just above 900°C. 
It is thought that there may be corresponding changes in 
the crystal structure of uranium at temperatures in the 
vicinity of 600 and 800°C. An attempt is being made to 
investigate the matter further by means of x-rays. 

W. L. 
R. Wricut 
H. B. WaHLIN 
Department of Physics, 
University of Wisconsin, 


Madison, Wisconsin, 
April 12, 1938. 


Cosmic Rays; Barometric Effect, Variations of Second Kind 
and Disturbances Produced by the Earth’s 
Magnetic Field 


Simultaneous records of the zenithal intensity of cosmic 
rays have been made during a period of nearly a year with 
large Geiger-Miiller counters in two separated double co- 
incidence arrangements. The new method of mutual control 
proved very reliable and rendered even semi-hourly values 
trustworthy as may be seen from the Figs. 1-3. In these 
the cosmic-ray intensity J reduced to 750 mm Hg together 
with the variations of the atmospheric pressure, the hori- 
zontal X (N-—S) and vertical Z components of the earth’s 
magnetic field are plotted separately for each arrangement 
for the month of January, 1938. I am very indebted to 
Professor Bartels and his staff for these magnetic data 
from the Adolf-Schmidt Observatorium Niemegk. 

In the first half of January, 1938 the records run nearly 
undisturbed both for J, X and Z contrary to the second 
half of the month, during which three great magnetic 
storms, each of different behavior, occur. At once the 
cosmic-ray intensity diminishes simultaneously and almost 
in the same manner as _X, and inversely as Z, but distinctly 
different each time as the magnetic storms themselves 
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Fic. 1. Cosmic-ray intensity J, horizontal X (N-S) and vertical Z components of the earth's magnetic field. 


vary. The first storm (Fig. 2) of January 16 exhibits a long 
perturbation of Stérmer’s toroidal stream and, in conse- 
quence of the large distance, a greater influence on cosmic- 
ray intensity; whereas the shorter storm (Fig. 3) with 
the extended aurora of January 25, occurring in lower 
regions, affected the rays not so much. Finally the second 
storm of January 22 brings only a relative small diminution 
of J; X, and Z changed at first inversely. It is to be hoped 
that closer investigation now in progress will reveal further 
details, also for the smaller variations of the earth’s mag- 
netic field. At any rate magnetic disturbances almost 
always seem to reduce cosmic-ray intensity contrary to 
the hitherto prevailing opinion. 

For the so-called “variations of the second kind” of 
cosmic-ray intensity, especially the very fast decreases 
which are followed by more or less slow increases, are 
characteristic features. The magnetic influences here re- 
ported perform in essentially the same manner so that 
these variations of the second kind are really the conse- 
quences of magnetic disturbances. 

The magnetic variations and their influence on cosmic- 
ray intensity reveals the variability of the barometric 
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Fic. 2. Cosmic-ray intensity J, horizontal X (N-—S) and vertical Z 
components of the earth's magnetic field. 
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Fic. 3. Cosmic-ray intensity J, horizontal X (N-—S) and vertical Z 
components of the earth's magnetic field. 
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Fic. 4. Daily mean values of atmospheric pressure, cosmic-ray 
intensities and magnetic field values. 


effect. Fig. 4 of the daily means both for intensity J (un- 
reduced) and atmospheric pressure shows a perfect recip- 
rocal behavior during the magnetic undisturbed period of 
January, 1938 contrary to the disturbed second half. If 
however the daily means of J are reduced with a constant 
barometric coefficient, derived from a long series, the in- 
tensity remains constant during magnetic undisturbed 
periods but changes at once when magnetic variations 
occur. 

The real barometric effect is therefore a pure absorption 
effect with a nearly constant coefficient and has nothing 
to do with the observed variations of the usually so-called 
barometric effect, which result both from magnetic and 
absorption phenomena perhaps also from other but less 
important influences. However these absorption effects are 
somewhat complicated, since they originate from the dif- 
ferent energies of the rays, as will be shown elsewhere. 

W. 

Institut fiir H6henstrahlenforschung 

der Universitat Berlin, 


Berlin, Germany, 
April 4, 1938. 
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Artificial Radioactivity Produced by Protons 


By the use of a beam of about one microampere of 
protons at an energy of 3.8 Mev, accelerated by the 
Princeton cyclotron,! we have found that Ni, Cu, Mg and 
Ag are strongly radioactive after bombardment. The 
following identifications are suggested for the radioelements 
formed. 

Ni: Four periods are found in the decay of Ni bombarded 
with protons. The half-lives are 80+2 seconds, 7.9+0.6 
minutes, 3.4+0.1 hours, and 12.8+0.3 hours. The relative 
initial intensities in a thick target, corrected to infinite 
bombardment time, are about 3 :1:31:35, respectively. The 
three shorter periods all emit positrons, while the 12.8-hour 
period has been found to emit both electrons and positrons. 
The last activity is then almost certainly to be identified? 
with Cu", formed here by the reaction Ni™(p, n)Cu®™. The 
3.4-hour period can also be identified with some certainty*® 
as Cu", formed from the stable isotope Ni®™ by the (/, 7) 
reaction. Two radioelements of approximately these 
periods were reported as formed by proton bombardment of 
Ni in the preliminary data of DuBridge and his collabora- 
tors,‘ but were not mentioned in their recent paper.® 

While chemical identification of the two short periods 
has not been carried out, the fact that positrons are 
emitted makes it seem very likely that the radioelements 
are Cu isotopes, formed by the (p, ”) reaction. Our decay 
curves give no evidence for a half-life of 10 minutes, which 
is that® of Cu®, so that one of the short lived radioelements 
observed probably arises from Ni®, and the other from 
Ni®*. We may tentatively assign the 80-second period to 
Cu*®* and the 7.9-minute period to Cu®, on the doubtful 
basis of relating the activity ratio of these two radioelements 
with the abundance ratio of Ni®* to Ni®. 

Cu: Copper bombarded with protons displays a moder- 
ately strong activity in which positrons are emitted. The 
half-life is 38+1 minutes. This radioelement may be identi- 
fied’ as Zn**, formed here by the reaction Cu®*(p, n)Zn®, 

Mg: Measurement of the activity of this element was 
not begun soon enough after bombardment to detect the 
presence of Al?*, which has a 7-second half-life’ and could be 
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Fic. 1. Absorption in Al of the radiation from silver which has been 
bombarded with protons (6.67-hour period). The mass range of the 
electrons is seen to be about 70 mg/cm?, while the mass absorption of 
the gamma-radiation is that of the Ka radiation of Ag. 
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formed here by the (p, 2) reaction. The decay curves of 
activated Mg showed a fairly strong period of 20+2 
minutes, together with a much weaker activity of some 
hours half-life. The short period radioelement emits 
positrons and has been chemically identified as an Al 
isotope. It is presumably either Al* or Al*®. 

Ag: The decay curves in this element show a weak 
activity of a few minutes half-life, which has not been 
investigated, and a very strong activity of half-life 
6.67 +0.06 hours. Chemical identification of the long lived 


‘ radioelement shows that it is an isotope of Cd, presumably 


either Cd! or Cd!. The particles emitted are negative 
electrons, as has been shown by deflecting them in a 
magnetic field, and their energies, as measured by ab- 
sorption in Al, are not greater than about 250 Kev. A 
very much weaker radiation with greater penetrating 
power (see Fig. 1) has a mass absorption coefficient in Al 
which is that of the Ka radiation of Ag. There is some 
indication of a harder radiation still which is of such small 
intensity that we have not been able to investigate its 
absorption with the Lauritsen electroscope used in these 
measurements. 

Since neither In! nor In"? is a stable isotope, we 
tentatively regard the negative electrons emitted by this 
radioelement not as disintegration electrons, but as 
photoelectrons ejected by a gamma-radiation which 
follows extranuclear electron capture by the unstable Cd 
isotope. If the absorption measurements are correct in 
indicating that the K radiation of silver is also emitted, the 
production of this radiation is presumably accomplished in 
two ways; by the filling of the K shell in Ag after the 
capture by the unstable Cd nucleus of one of its own K 
electrons, and by the replacement of electrons which have 
been removed from the atom by the internal conversion 
process postulated to explain the negative electron emission 
from this radioelement. The harder radiation suggested by 
the absorption curve may be the unconverted fraction of 
the nuclear gamma-radiation postulated above. Rough 
considerations of intensity lead to the conclusion that the 
internal conversion coefficient for this gamma-radiation 
must be about 90 percent. The whole situation in this 
radioelement is similar to the case of Ga®’ discussed by 
Alvarez.°® 

We are glad to be able to thank Mr. R. C. Newton, of the 
department of chemistry, for suggestions concerning the 
chemical separations, and Mr. John E. Walter for assist- 
ance in making the measurements. 

Louis N. RIDENOUR 
L. A. DELsAsso 

M. G. WHITE 
RuBBY SHERR 


Palmer Physical Laboratory, 
Princeton University, 
Princeton, New Jersey, 
April 15, 1938. 
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Concentration of Radioactive Copper by Use of a High 
Speed Rotating Cathode* 


Radioactive copper (period 12.8 hours), formed from zinc 
by neutron bombardment, has previously been concen- 
trated in three ways: addition of inactive copper ion, 
followed by precipitation of CuS;! addition of inactive 
copper ion followed by deposition on a zinc plate;! and 
partial solution of zinc dust in acid, followed by electro- 
chemical displacement on a lead plate.? This communica- 
tion describes a fourth method, which may find wide 
general use in problems of concentrating minute quantities 
of such isotopes. The high speed rotating cathode,’:‘ it 
has been found, will separate a considerable proportion 
of the radioactive copper from solution, without the 
addition of its nonradioactive isotope as a ‘‘carrier.”’ It is 
thus possible to obtain extremely thin layers of radio- 
active copper on a variety of surfaces, including copper 
itself. 

The procedure consists of irradiating overnight 90 grams 
of zinc dust with a radon-beryllium source (300 to 600 
millicuries) of neutrons. The zinc is then treated with 
concentrated hydrochloric acid until all but five grams 
have dissolved, which undissolved part carries practically 
all the activity.2 This is next dissolved in hydrochloric 
acid ; 30 to 40 cc of concentrated sulphuric acid are added, 
and the mixture evaporated until fumes of SO; are given 
off. It is then diluted to 400 cc and filtered. This solution 
is then electrolyzed with the use of a rapidly rotating disk 
cathode of copper, platinum or other metal, of 3 cm 
diameter, 0.25 mm thickness. With a platinum anode, 
under these conditions, satisfactory operation was ob- 
tained with currents of 10 to 18 milliamperes and voltages 
ranging from 2.1 to 2.7 volts. The optimum speed of 
rotation of the cathode was about 4500 r.p.m. 

Activities of various disks, sulphide precipitates, etc., 
were measured with a thin silvered-wall Geiger counter 
or with a pressure ion chamber electrometer system. 

The following preliminary results are reported: 

(a) The amount of radioactive copper ion deposited is 
approximately a linear function of the duration of the 
electrolysis in the initial period up to 30 minutes. 

(b) The effect of the stirring rate on the removal of 
activity is such that there is a maximum between 2300 
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and 5000 r.p.m. The speeds tested ranged from 2300 to 
12,000 r.p.m. 

(c) The efficiency of cathode materials in the removal 
of radioactive copper decreased in the order: copper, 
molybdenum, platinum, tungsten, tantalum. 

(d) It is estimated that a copper disk, running for 45 
minutes, at 4500 r.p.m., 11 milliamperes and 2.1 volts, 
removed at least 50 percent of the activity in solution. 
The estimate is based upon a comparison of the activity 
of an aliquot of solution before electrolysis as CuS, and 
of the activity of the solution after electrolysis as CuS. 

(e) The high speed rotating cathode gives very sharp 
separations. When a current of 10 to 11 milliamperes was 
used, it was found that the appearance of the disk after 
electrolysis was unchanged, and retained its mirrored 
finish, with no evidence of the deposition of other elements. 
At higher currents the appearance of the disks after 
electrolysis indicated that some impurities had plated out 
as well. They then presented a somewhat darkened 
appearance, although still retaining the mirrored surface 
which had been imparted to them by buffing. 

This method permits a much higher concentration ratio 
of active to inactive isotope to be achieved than under 
ordinary conditions, since it eliminates the need for the 
addition of the inactive isotope to active solutions. It 
should be applicable to a large number of elements. It can 
facilitate the study of various diffusion and exchange 
phenomena and any researches which require high con- 
centration and purity of the radioactive indicator. 
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